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INTRODUCTin^ 



^^Igh School mathematics teachers have long been aware that their students 
HBshould know not only something about the development of pure mathematics 
Wmbut also something about its applications. Several years ago NASA, recog- 
iizing the appeal of aerospace activities, initiated and supported the development 
of curriculum supplements for several high school courses. Because attainments in 
aerospaee would not be possible without mathematics, it was most appropriate 
that Q iupplementary publication deaUng with space activities be prepared for 
teachers of mathematics. 

The first mathematics curriculum supplement. Space Mathemaiics, A Resource far 
Teachers, was published in 1972. One of the most popular and oft-requested of the 
supplements, the book has been unavailable for several years. This volume up- 
dates the earlier work. We hope that a new generation of students will become 
more interested in mathematics as the result of seeing some of its significant appli- 
cations in recent and current space projects. Working problems such as those in 
this book should enhance both the mathematical knowledge and skills of students 
and their appreciation and understanding of aerospace technology and achieve- 
ments. 

NASA's Technical Monitor for this project was Muriel M, Thorne, Educational 
Programs Officer, under the general direction of WilHam D. Nixon, Chief of Edu- 
cation Services. NASA. 

National Aeronautics and Space Administration 
Washington, D.C. 
September 1985 
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PREFACE 



■ n 1972, a collaction of mathematicai problems related to space science entitled 

■ Space Mathematics, A Resource for Teachers was published by the Educational 
H Programs Division of the National Aeronauiics and Space Administration 
(NASA). As an early user of that publication, I can say that it has been both a 
pleasure and a challenge for me to undertake the revision of that volume of en- 
richment materials, especially in the light of another twelve years of activity in 
space exploration. This interval has been a period of much progress in both the 
science and the technology associated with the space program, and it has offered a 
wealth of new material in which to find applications of high school mathematics. 

The basic format of the original publication has been retained, as well as many of 
the classical problems and those which complemented the new material. In deveU 
oping the examples and problems presented here, we have aimed at preserving the 
authenticity and significance of the original setting while keeping the level of math» 
ematics within the secondary school curriculum. The problems have been grouped 
into chapters according to the predominant mathematical topic. Within each chap- 
ter we have attempted, as far as possible, to group problems involving similar 
themes. There is a wide range of sophistk*ation required to solve the various prob- 
lems. Since this is a resource book for teachers, we have assumed that the reader 
will be interested not only in problems that can be brought directly into the class- 
room, but also in those that, although beyond the current level of their students, 
will increase the teacher's own awareness of some of the interesting applications of 
mathematics in the space program. 

Perhaps the most valuable potential of a collection such as this lies in its ability to 
convey a sense of how secondary school mathematics is actually used by practicing 
scientists and engineers. Attitudes and approaches may thereby be fostered, on 
the part of teachers, that can help students to be more insightful users of the 
mathematics they learn. The present school mathematics curriculum, for example, 
gives no hint that many real-world problems do not have analytic solutions in 
closed form but may nevertheless be satisfactorily "solved'' by using carefully cho- 
sen approximations or the numerical methods made possible by modern computers. 

In this connection, we stress that in order to use numerical analysis correctly or to 
make good approximations, it is necessary to know something of the theoretical 
background of the subject and to understand the concepts of precision and accu- 
racy md the use of significant digits. Also, methods that reveal meaningful aspects 
of a procedure are preferable to purely algorithmic prescriptions; the perhaps un- 
familiar ^'factor unit'' method of unit conversion presented in Chapter 2 is actually 
quite commonly used in science and engineering. It not only removes all uncer- 
tainty about whether to multiply or divide by a conversion factor but also is far 
more likely to contribute to an understanding of the underlying concepts than, for 
example, the more usual metric system algorithm expressed in terms of "moving" 
the decimal point. 
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CHAPTER GNE 



MATHEMATICAL ASPECTS OF 
SOME RECENT NASA MINIONS 




^peachers of mathematics, like most adults in today's world, can hardly fail to be 

■ aware of the rapid development of space science, We realiza that the spec- 

■ tacular achievements of the space program have depended heavily on mathe- 
matics—mathematics that is generally complex, advanced, and well beyond the 
level of most secondary school curricula. Even though this perception is valid, 
there are many significant aspects of space science that can be understood using 
only high school mathematics. 

The exploration of space naturally uses the tools and techniques of astronomy. 
Astronomy in turn is gaining much new information as a result of sending sci- 
entific probes and satellites beyond Earth's atmosphere. Because astronomy has 
stimulated the growth of many of the concepts and methods of mathematics, the 
high school teacher will find here much that is familiar. However, in some in- 
stances the way mathematics is used to solve reaUlife problems is rather different 
from methods emphasized in school courses. 

In this opening chapter, we shall examine several recent achievements of the 
National Aeronautics and Space Administration and identify mathematical ideas 
and questions that may be of interest to high school teachers and students. When 
appropriate, we will refer to a problem illustrating some aspect of the subject 
and worked elsewhere in the book. 



The Space Shuttle 

The Space Shuttle (Fig. 1.1) is a true aerospace vehicle— it takes off like a rocket, 
operates in orbit as a spacecraft, and lands like an airplane. Tb do this takes a 
complex configuration of three main elements^ the Orbiter, a delta-winged 
spacecraft=aircraft, about the length of a twin-jet commercial airliner but much 
bulkier; a dirigible-like external tank, the only expendable element, secured to the 
Orbiter's belly and containing two million liters of propellant (Chapter 4, Prob^ 
lem 5); a pair of reusable solid rocket boosters, each longer and thir*:er than a 
railway tank car and attached to the sides of the external tank. 



Each Space Shuttle is meaat to be just one element in a total transportation 
system linking Earth with spac^ . In addition to providing for continued scientific 
investigations by transporting such systems as the Spacelab and the Large Space 
Telescope, recently renamed the Edwin P. Hubble Space Telescope, into orbit 
(Chapter 3, Problem 4), the Space Shuttles are also expected to carry the build- 
ing blocks for large solar-power space stations or huge antenna-bearing structures 
for improved communication systems (Chapter 4, Problems 9 and 10), Structures 
that would be too fragile to stand up under their own weight on Earth will be 
folded up in the Shuttle's cargo bay and assume their final shapes in the micro- 
gravity environment of space. The Shuttle will also be capable of carrying a work 
force of seven people and returning them home after the completion of their work. 
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Mathematical Aspaets of Some Recent NASA Missions 



One of the most basic mathematical problams raised by the launching and control^ 
ling of a Shuttle or any other spacecf aft is that of describing its motion. This 
problem requires the ability to specify the position of the spacecraft's center of 
mass and its attitude (ofientation) and to describe changes in both during flight. 
The specification of position and attitude can be accomplished by setting up suit- 
able coordinate systems (Chapter 7, Problem 10), Instruments to determine a 
spacecraft's attitude are most effectively referenced to a spacecraft-based coordi- 
nate system, whereas ground control is best accomplished in terms of an Earth- 
based system. This duaUbased system necessitates transformations between 
coordinate systems (Chapter 7, Problem 1, and Chapter 8, Problem 2). 
Describing a change of position and attitude requires an understanding of the 
measurement of time (Chapter 2, Problem 11). It is interesting to note here that 
our definition of a day on our rotating Earth must be redefined for a Space Shuttle 
Orbiter crew. For them the Sun might rise again and again every hour and a half! 

The Planetary Probes 

The launch of the two Voyager spacecraft in the summer of 1977 climaxed a series 
of fruitful missions of planetary exploration including the Manner, Viking, and 
Pioneer series of probes to Mercury, Venus, Mars, Jupiter, and Saturn. All these 
missions sent back new information about the structure and composition of these 
planets and their associated moons. We focus in this book on some of the results 
of Voyager J and Voyager 2, These probes, which benefited from more highly 
developed instrumentation and computer capability than their predecessors, 
approached closer to Jupiter (Chapter 7, Problem 11) and Saturn than previous 
flights did. Stunning pictures resulted, showing the unanticipated presence of 
active volcanoes on Jupiter's moon lo (Chapter 10, Problem 6) and the fine 
structure of Saturn's rings. 

Among the mathematical problems that arose in these missions were the 
following. 

1. Transmitting spacecraft observations back to Earth (Chapter 5, Problems 2 
and 3, and Chapter 8, Problem 1). 

2. Determining the time of transmission of spacecraft observations (Chapter 3, 
Problem 5). 

3. Calculating the rotation period for planets such as Saturn, which is not solid 
and has no outstanding observable features like Jupiter's Great Red Spot (Chapter 
2, Problem 13). 
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Satellites 

NASA began its formal existence in 1958 and by the end of 1979 had iuecessfully 
launched more than 300 large and small satellites with missions as diverse as 
observing Earth's weather (Synchronous Meteorohgical Satellite [SMS] series) and 
resources {Landsat series), providing comrnunication links for television signals 
(Applications Technology Satellite [ATS] series), and measuring solar radiation out- 
side Earth's atmosphere (Orbiting Solar Observatory [OSO] series). 

The design of these satellites and their experiments and the analysis of the data 
gathered involve a variety of mathematical questions. We shall consider some of 
the following examples. 

1. The connection between the conic sections and the law of gravitation 
(See Appendix). 

2. For elliptic orbits, the connection between the orbit parameters and the period 
of revolution (Chapter 9, Problem 11) and the determination of the exact posi- 
tion of a satellite in its orbit at a specified time (Chapter 9, Problems 19 and 20). 

3. The geometry necessary to correct for distortions arising when flat pictures 
are made of a curved Earth (Chapter 7, Problems 7 and 9, and Chapter 10, 
Problem 2). 

4* The need for logarithms to understand how radiation is absorbed by Earth's 
atmosphere (Chapter 6, Problem 3)* 

The mathematical analysis of the reflective properties of the conic sections 
needed to design an X-ray telescope (Chapter 9, Problems 21 and 22). 

6, The judicious use of approximation (Chapter 3, Problem 8; Chapter 4, 
Problems 6 and 8; Chapter 7, Problem 6; Chapter 9, Problem 22). 
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MEASUREMENT 
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A photograph of the planet Jupiter made 
from Images obtained by Voyager 1 on Feb- 
ruary 5, 1979, showing the Great Red Spot 
and three of Jupiter's four largest satellites: 
lo (In front of Jupiter), Europa (brightly lit, to 
the right), and Callisto (barely visible at the 
bottom left). 
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Chapter Two 



I pace science is based on a mathematical description of the universe. This 
. inathematical description is in turn based on defining physical quantities 
clearly and precisely so that all observers can agree on any measurement of 
these quantities. Every measurement has two parts: a number and a unit. In 
mathematics, wa tend to focus on the numbers and assume that the units are taken 
care of; but in scientific work, units receive careful attention through a proce- 
dure known as dimensional analysis, which is illustrated in the first problem. 

Among the physical quantities used to descfib^^ the universe, some are considered 
fundamental quantities whereas others are derived quantities, comparable to the 
designation of definitions and undefined terms in a mathematical system. Although 
it does not really matter which particular quantities are the ones designated as 
fundamental, the most common are length, mass, and time. In scientific work the 
two major systems of unitf for these quantities are the mks (meter-kilogram- 
second) and the cgs (centimeter gram-second). Every measurement is a compari- 
son with the standards that are universally accepted as definitions of these funda- 
mental units. In astronomy and space science, where large distances are common, 
the meter and even the kilometer are too Small to be convenient; in Problems 5, 
9, and 10 of this chapter, we show how more suitable units for length are defined. 

Dimensional analysis (manipulation of units according to the rules of algebra) is 
the procedure used to ensure consistency In the definition and use of units. For 
©xample, since force is, by definition, the product of mass and acceleration, mea- 
sured respectively in kg and m/s^ in the mks system, the unit of force in this 
system must be equivalent to kg - m/s^. A new term, the newton, was created to 
describe the unit of force: 1 newton = 1 kg - m/s^. 

PROBLEM 1, Newton's law of gravitation, one of the most important ideas in space science, states 
that the force of gravitational attraction between two bodies of masses M, and Mi 
is proportional to the product of the two masses and inversely proportional to the 
squap of the distance R separating the two masses. If G is the constant of pro- 
portionality, called the tdniversal gravitational consiant, this law can be stated in 

symbols as F ^ — - What must be the unit for G in the mks system? 

Solution: Using dimensional analysis, we equate the known units in accordance with the 
relationship above without worrying about the numbers, then solve algebraically 
to get the unknown unit. This gives 

newton ^ (unit for C) (kg) (kg) (m)^^ 

or (unit for G) ^ (newton) (m)^ (kg)^^ 

- (kgmsec-2j(^32(^gj- 

= m^- kg°* sec-. 
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Co mputation and Mea surement 



PROBLEM 2. 



Solutioni 



PROBLEM 3. 



Solutioni 



We know that in a circle of radius r, if an arc of length 3' subtends an angle 6 and 6 is 
measured in radians, then s ^ rd. Show that the radian is essentially dimension- 
less (i,e., an angle of ^/4 radians is just the real number 7r/4), 

Since r and s are both lengths, in the mks system they will both be measured in 
meters. From s ^ rO we have 6 = f — ^ , Since the units cancel, & is 



dimensionless. 

Scientific theories and technological development both require accurate mea- 
surements. Since every measurement is an approximation, an important aspect of 
scientific and technical work is the analysis of experimental error and the control 
of the propagation of error when computations are made using measured quan- 
tities. The use of computers to solve complex problems by numerical methods 
has made error analysis even more important because computers approximate real 
numbers using finite decimals. Moreover, computers represent numbers inter- 
nally using a floating point binary representation. Even though it is not really 
necessary to understand the binary numeration system to work with computers, 
such knowleuae is essential to the analysis and control of error propagation in 
computational work. The next problem considers the floating point binary repre- 
sentation of our familiar numbers. 

The binary (i,e,, base two) representation of a number U'jes only two digits, 0 and 1 , 
Whereas in base ten the actual value of a digit is the product of its nominal value 
and the appropriate power of 10 according to the position of the digit with respect 
to the decimal point, in base two the value of a digit is the product of its nom- 
inal value and the appropriate power of 2, So, for example, the binary number 
10011 has the value that we represent in base ten as 1 x 2^ + 0 x 2^^ + 0 x 2^ + 
1 X 2^ + 1 X 2-\ or 16 + 2 + 1 ^ 19; the binary number 10.011 is the same as the 
decimal number 1 x 2* + 0 x 2" + 0 X 2^* + 1 x 2^^ + 1 X 2^\ or2 + 0,25 + 
0,125 ^ 2,375, 

a* Determine the binary representations of the decimal numbers 625, 6.25^ and 
0.0625. 

625 can be written as the sum 



512 + 64 + 32 + 16 + 1 ^ 2'' + 2*^ + 2'' + 2"* + 2", 



so 625e,n ^ 



1001110001 



6,25 can be written as the sum 



4 + 2 + 



= — 2^ + 2^ + 2^', 



4 




110.01 



0.0625 - 



625 



16 ^ 



10000 



so0.0625e,.n = 0,0001, 

15 
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Chapter Two 



Show that it is impossible to represent the decimal fraction 0,2 exactly in a finite 
binary code. 

Solution: 0.2 ^ 2/10 ^ 1/5. To express this in binary notation, we must write 1/5 as a sum of 
unit fractions, each having some power of 2 as denominator. Since 1/2^^ is the 
largest such fraction smaller than 1/5. we begin by finding the difference: 

I _ ±^ 8-5 ^ 3 
5 2'^ 5 ' 2^^ ~ 5 ' 2^ 

Now the largest unit fraction less than 3/(5 - 2') with a power of 2 as denominator 
is 1/2^, so we next find the difference: 

3 _ X ^ 6-5 ^ 1 
5-2^ 2^ 5 ' 2^ 5 * 2^ 

This means that 1/5 - 1/2' + 1/2' + (1/2-^) (1/5). Since the fraction 1/5 has re- 
curred, multiplied by 1/2', we see that the first four digits we have found to the 
right of the binar y '^de cimal point/' 0,001 1 , will repeat continuously. In other 
words, 0.2tg„ ^ 0,001 1^^^. (A quicker but less intuitive approach to finding this 
representation is to express 1/5 as the binary fraction 1/101 and then divide 101 
into 1, using binary arithmetic.) 

The reader can use the method of part (b) to show that the decimal fractions 0, 1 . 
0.3, 0.4, 0,6, 0.7, 0.8, 0,9 also have infinitely repeating binary representations. 

c. Almost all computers use a floating point binary representation for numbers. 
In this system, every number is expressed in the form 0,djd2 . . , dn x 2"", where 
dj = 1, d| ^ 0 or 1 for i - 2, 3, , . . , n, and m is an integer. For example, the 
floating point representations of the numbers in part (a) would be 

625 ^ 0.1001110001 X 2'" 

6.25 = 0,11001 X 2' 

0.0625 ^ 0.1 X 2^'^ 

Different computers have differing capabilities both with respect to the length (n) 
of the string of O's and 1 's that can be stored for any single number and with 
respect to the exponent m that can be stored. The limits available for n and m 
determine the largest and smallest number a computer can represent and also 
the size of the errors that must result when a number with an infinitely repeating 
representation must be stored with only a finite string length available. 

If a certain computer can store only an eight-digit string (n - 8), then the repre- 
sentation for the decimal fraction 0,2 will be stored as 0. 1 1001 100 x 2^-. What 
number is this, and what is the difference between this number and 0.2? 
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Co mputiition an d Measurement 



Solutioni 



8 16 



128 



256 



32+16 + 2+1 
256 



51 

256 




256 - 255 

' 5 • 25'6"" 



1280 



1 



We now state two definitions used in error analysis. These definitions can be 
applied to both measurement errors and the errors that arise because of the way in 
which numbers are represented in computers. It is probably worth noting in this 
context that the term measurement error as used here does not imply that the 
measurement has been carelessly made but rather refers to the fact that every 
measuring instrument is limited in accuracy and can never provide more th/in an 
estimate of a true value. 

Let be the true value of a specified quantity, and let X be the value of this 
quantity as measured or as represented in the computer. Then: 



Observe that absolute error has the same units as the quantity under considera- 
tion, whereas relative error (usually reported as a percent) is dimensionless. 

The relative error is considered to be the indicator of how good a measurement 
or any other approximation is. For example, a measurement of 2.5 mm with a 

possible absolute error of 0.05 mm has a relative error of ^ . or 2 percent. 

whereas a measurement of 1250 km. with a (much larger) possible absolute error 

of 5 km has a much smaller relative error of y-^- , or 0.4 percent. Awareness of 

the appropriate tolerance for relative error is a vital ingredient of scientific work. 



absolute error in X 



Xi - X 



relative error in X 



Xt - X 
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Chapter Two 



PROBLEM 4. What are the absolute and relative errors if a computer that has an eight^bit binary 
digit string represents 0,2 as 0, 1 1001 100 x 2^' ? 

Solutloni From Problem 3c, the absolute error is 1/1280. or about 0.0008. The relative 
error ^ 0,0008/0.2 ^ 0,004, or 0.4 percent. 

The use of significant figures is helpful in error analysis. The number of signifi- 
cant figures is defined as the number of digits that can be assumed to be correct, 
starting at the left with the first nonzero digit, and proceeding to the right. By 
this definition, 10.62, 0,05713, and 4.600 all have four significant figures. A num- 
ber such as 4300 is ambiguous. This ambiguity may be resolved by using scientific 
notation, since we may write ihe number as 4.3 x 1Q\ 4,30 x 10\ or 4,300 x 10^ 
according to whether the number has two, three, or four significant 
figures, respectively. 

When approximate numbers are added or subtracted, it can be shown that the 
absolute error in the sum or difference could be as large as the sum of the abso- 
lute errors of the individual numbers. When approximate numbers are multiplied 
or divided, it can be shown that the relative error of the result could be as large 
as the sum of the relative errors of the individual numbers. This means that for 
sums and differences of approximate numbers, the number of decimal places 
considered significant can never be greater than the number of decimal places in 
the least precise addend. For products and quotients, the number of significant 
figures can never be more than the smallest number of significant figures in the 
individual factors. Wherever appropriate, numerical results will be given in 
accordance with these guidelines, 

PROBLEM 5, Earth's orbit around the Sun is elliptical, but in many cases it is sufficiently accurate 
to approximate the orbit with a circle of radius equal to the mean Earth-Sun dis- 
tance of 1,49598 X 10^ km. This distance is called the Astronomical Unit (AU). 
Listed in the chart that follows are actual Earth-Sun distances, given to five sig- 
nificant digits, on the first day of each month of a representative year, (The Ameri- 
can Ephemeris lists daily distances and the actuaj times for these distances to 
seven significant digits,) 



Date 


Distance (x 10** km) 


1 January 


L4710 


1 February 


1.4741 


1 March 


1,4823 


1 April 


1.4949 


1 May 


L5073 


1 June 


L5169 


1 July 


1,5208 


1 August 


1.5183 


1 September 


1,5097 


1 October 


L4977 


1 November 


1,4848 


1 December 


1,4751 



a. To how many significant digits is It reasonable to approximate the Earth-Sun 
distance as though the orbit were circular? 
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Computation and Meaiurement 



Solutions To two significant digits, each of the distances in the table can be given as 
L5 X 10^ km. 



Solutloni 



PROBLEM 6, 



Solution: 



What are the largest possible absolute and relative Ui.ui. m using the Astro- 
nomical Unit as the Earth-Sun distance in a computation instead of one of the 
distances from the table? 

(1,49598 - 1,4710) x 10^ ^ 0,0240 x 10^ km (smallest table value) 

(L49598 - L5208) x 10^ ^ =0.0248 x lO^km (largest table value) 

absolute error s 0,0248 x 10^ km 

0,0248 



relative error s 



L49598 



0,0166, or 1,7 percent. 



The procedure of dimensional analysis^ described earlier, is easily adapted and 
commonly used in science and iechnology for the task of unit conversion. Recall 
that in dimensional analysis the units are manipulated in accordance with the 
rules of algebra. 

Suppose we wish to change a length of 623 cm to meters. The adaptation of dimen- 
sional analysis for unit coversion involves multiplication by b factor chosen 
according to the following simple principles: the factor unit is a fraction with a 
value of 1, whose numerator is expressed in terms of ihe unit we wish to have and 
whose denominator is expressed in terms of the unit we wish to change. Since 
100 cm ^ 1 m, in order to change 623 cm to m, we perform the multiplication 



623 cm 



1 m 



, "canceling" the cm in numerator and denominator to get 



1 100 cm 
623/100 m, or 6,23 m. 

More complex conversions can be done using multiplication by several factor 
units and those readers wishing to convert between British and metric units can 
also use this method. For example, the speed of light, 3,00 x 10^ km/sec, can be 
found in miles per hour: 

3.00 X 10^ km Imile 60 sec 60 min ^ ^ ■ 
— - — — — X ^ ^ r- — X — =r- X = 6,71 X 10^ miles per hour. 

1 sec 1,61 km 1 mm 1 hour 

The deep space probe Piorieer 10 took 21 months to get from Mars to Jupiter, a 
distance of 998 million kilometers. Use the factor unit technique to find its aver- 
age speed in kilometers per hour during that period. 



average speed 



distance 
time 

998 X 10^ km 



^ 12 months ^ 1 day 
21 months 365 days 24 hours 

« 6.5 X 10^ km/h, or about 65 000 km/h 
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PROBLEM 7. 



Solution: 



a. Recall that the Astronomica! Unit (Earth^Sun distance) discussed in Problem 5 Is 
L^^6 X 10' km, to four significant figures. Find the Earth^Sun distance in miles 
to three significant figures, 

1 AU - 1,496 X 10^ km 



L496 X 10^ km 1 mile 

^ = ~ K 



L61 km 



- 9,29 X 10^ miles (almost 93 million miles) 

b. The chart that follows gi.es the mean distance in kilometers of each planet in 
the solar system from the Sun, Hxpress these distances in AU, using a suitable 
number of significant digits. 



Plane t 


Disiance 


Mercury 


0,579 


Venus 


L08 


Mars 


2.27 


Jupiter 


7,7S 


Saturn 


14.3 


Uranus 


28.7 


Neptune 


45.0 


Pluto 


59.1 



(km X IQ^) 



Solution: Smca eacn distance m the table has three significant digits, and the factor unit 
1 AU 

has an exact number in the numerator and six significant digits 



1,49598 X 10^ km 

in the denominator, the distances in AU can be given to three significant digits. 
Multiplying by the factor unit shown gives the following distances in AU: 



Mercury 


0.387 


Venus 


0.722 


Mars 


L52 


Jupiter 


5.20 


Saturn 


9,56 


Uranus 


19.2 


Neptune 


30.1 


Pluto 


39.5 



PROBLEM 8. The Solar Maximum Mission (SMM) satellite orbits Earth at a height of 560 km In 
many computations, the Earth-Sun distance of 1.5 x 10' km is used to approxi- 
mate the distance of SMM from the Sun. What is the maximum relative error of 
this approximation? 
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Solution: The distance of SMM from the Sun is contained within the range (Earth-Sun 
distance ± (Earth diameter ^ 560 km)), or (Earth-Sun distance ± 6930 km). 
From Problem 5, if 1.5 x 10^ km is used as the Earth-Sun distance, the absolute 
error s2,48 x 10^ km. 6930 km is much smaller than this error, so that the abso^ 
lute and relative errors incurred in using L5 x 10'^ km as the SMM-Sun distance 
are the same as those of part (b) of Problem 5. If greater accuracy is required, it 
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PROBLEM 9. 



will be necessary to use daily Ephemeris values such as those in the listing given in 
Problem 7. In this event, it is stUl true that since 6930 km < 7 x 10^ km, the 
relative error in approximating the SMM-Sun distance with the Ea:th-Sun 



distance < 



7 y W 
L5 X 10^ 



= 4.67 X 10^\ or about 0.005 percent. 



The Astronomical Unit (AU), although useful for measuring distances within the 
solar system, is too small to be convenient for distances to stars. We shall therefore 
consider two other units of length used by astronomers. The first is called 
the lighhyear. 

The light-year is the distance travelr-d by light during one Earth year. To three sig- 
nificant digits, the speed of light is 3,00 k 10^ km/s. Find the length of the light- 
year in km and in AU. 



Soiutioni 1 Earth year = 365.25 days ^ 365,25 x 24 x 60 x 60 seconds. In one year, light 
travels 3.00 x 10- x 365.25 x 24 x 60 x 60 km ^ 9,47 x 10^= km. To express^this 

1 A(J 



distance in AU, 1 light-year — 9,47 x 10^^ km x 



1.50 X 10^ km 



= 6.31 X 10- AU. 



EKLC 



The parsec is the astronomical unit of distance that relates to observational mea- 
surements. In order to define this unit, we must consider the fact that when we 
observe the heavens, we have no direct perception of depth or distance, A useful 
model developed to portray the heavens is the celestial sphere. In this model. 
Earth is surrounded by an imaginary sphere with infinite radius, A coordinate 
system, similar to latitude and longitude, is imposed on the celestial sphere by 
projecting Earth^s rotation axis on the sphere to identify the celestial north pole 
(CNP) and celestial south pole (CSP) as shown in Fig. 2.1. Since the radius of the 
celestial sphere is infinite, all parallel lines point to the same spot on the sphere, 
and so every line parallel to Earth*s rotation axis also points to the celestial north 
and south poles. 

The extension of Earth's equatorial plane intersects the celestial sphere in a 
great circle called the celestial equator. Now a system of small circles of declin- 
ation (5), comparable to latitude ciicles on Earth, is imagined on the celestial 
sphere, and a system of great circles called right ascension (a) circles, comparable 
to longitude, passing through the two poles, completes the coordinate system 
(Fig, 2,2), 

Celestial norih po!ej 

8 ^ +90^ 

Fig* 2.1 




CSP 



2 ^^'^^^'^ 




1 equator; 5=0^ 



Vernal equinox 



23 



Chapter Two 



Every star or celestial object can now have its position identified bv the ordered 
pair (a, 6). Because Earth rotates with respect to the celestial sphere, the time of 
observation must also be known in order to use the coordinate system. Diffefences 
in th^ positions of two objects on the celestial sphere are expressed in terms of 
the angle subtended at Earth by the arc joining these points. 

As Earth revolves around the Sun, very distant stars show no discernible changes 
in position, but closer stars will show apparent motion with respect to the celes- 
tial sphere when viewed from different points in Earth's orbit, as shown in Fig. 2,3. 
This apparent motion is called parallactic motion, and the change in position is 
called the paraliax angle. In this context, 1 parsec is defined as the distance at 
which the radius of Earth's orbit subtends an angle measuring 1 arc-second 
(see Fig, 2.4). 




Fig. 2.3 



PROBLEM 10. a. Find the length to three significant digits of 1 parsec in terms of AU, km 
and light-years. 

Solution: If d is in radians, we have arc length ^ r^, where r is the distance expressed in the 
same units as the arc length. In this case, arc length * radius of Earth's 
orbit ^ 1 AU. 



e - 1 second - ^ degree - IMMSI ^ ^ ^^^'^"^ 
3600 ^ 3600 180 degrees 

Since we want three significant digits in our answer, let us use tt = 3,142 in 
this computation: 

' = MmWm = 4.85 X 10- 



(We have omitted rad, since the radian is really dimensionless. (Problem 2, 
this Chapter)) 
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Then 



arc length 



1 AU 



2.06 X 10' AU. 



In terms of km, r ^ 2.06 x 10^' x 1.50 x 10^ km ^ 3.09 x 10'^^ km. Since 1 light- 
year ^ 6.31 X 10^ AU, 1 parsec ^ 2.06 x 10' AU x ^ I'ght-year ^ ^ . 

^ " 6JlxlO-*AU -^'-^'"gni 

years, 

b. In general. If p is the parallax of a star and d its distance from Earth, then 
d (in parsecs) - ^ .g^^nds of arc) ^ P^^^ of our nearest star, a Centauri, 
is 0.75 seconds, and the parallax of Sirius, one of the brightest stars in the north- 
ern sky, is 0,38 seconds. Find the distances to these stars in parsecs and in km. 

Solutfon: For c£ Centauri, 

d — Q parsecs — 1,3 parsecs 
^ 1.3 X 3,09 X 10^' km 



^ 4.1 X 10*- km. 

For Sirius, 

d — parsecs ^ 2.6 parsecs 

- 2,6 X 3.09 X IQi^^ km 

- 8,1 X 10^^ km. 

The accurate measufement of time has been one of the most challenging prob- 
lems in human history. We now tend to take for granted the civil time-keeping 
system in general use. This system has evolved over many centuries and from, 
time to time has been substantially revised. The original definitions of day, month, 
and year depended, fespectively, on observations of the periodic motions of the 
Earth, Moon, and Sun with respect to the celestial background as observed from 
Earth. Since all these motions have fluctuations in their periods, it is not possible 
to define a completely regular unit on which to base an accurate time measurement 
in terms of the day, month, or year. 



The first time-keeping instrument that did not depend on celestial observation 
was the pendulum clock. It did, however, depend on the Earth's gravity, which 
varies v^ith geographic location and the positions of the Sun and Moon, The 

second, origmally defined as 34 x 60 x 60 ^ ^^^^ recently has been 
redefined in terms of the microwave emissions of certain atoms (e,g,, quartz 
crystals). This new definition provides a uniform, standard with which to measure 
intervals of time. 
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PROBLEM 11. a. One of the serious problems of the old nonu.nifofm time units was the accumu- 
laiion of error. It might seem that an accuracy of 1 second a day (a possible rela- 
tive error of 24 x 60 x 60 ' ^^out 1 x 10 would be sufficient for most tech^ 
nical or scientific purposes. Show that an error of 1 second a day could result in 
an error of L 1 x 10"* km in the position of the Earth in its orbit after only 1 yesr. 
(Assume Earth's orbit is circular with a radius of 1.5 x 10^ km,) 

Solutions In one year the error could be 365 seconds. Earth moves through an angle of 

365 X 24I 60 X 60 '^^°"^' 24 x 60 x 60 ^" 

onds. The length of arc that subtends this angle in a circle of radius 1,5 x 10^ km is 

^ = " ^ 24 X 60 X 60 

b. The tropical year is defined as the time difference between successive vernal 
equinoxes— in other words, the time it takes Earth to complete one revolution 
around the Sun. This time does not have a simple relationship to Earth's rotation 
period (the day). In fact, it turns out that to the nearest second one tropical year is 
365 days, 5 hours, 48 minutes. 46 seconds, Sh:?w that the current system of 
adding an extra day to each calendar year that is a multiple of 4 but not a multiple 
of 100 (leap years) serves to give each calendar year an integral number of days 
and also keeps the seasons constant with respect to the calendar. 

Solutions If a calendar year has 365 days, the excess time in a tropical year is 5 h 48 m 46 s, not 
quite 1/4 day. Multiplying this excess by 4, 4 x (5 h 48 m 46 s) - 23 h 15 m 4 s/ 
almost 1 day. If we add 1 extra day each 4 years, we will create a deficit of 

24 h - (23 h 15 m 4 s) = 44 m 56 s for each leap year. In each 100 years, there are 

25 years that are multiples of 4; however, after 24 leap years, the deficit will accu- 
mulate to 24 X (44 m 56 s) ^ 17 h 58 m 24 s, almost 3/4 day. This will almost 
balance the excess accumulation for the remaining 4 years of the century, so that 
years that are multiples of 100 should not be leap years. It is clear that further 
juggling will be necessary, since things never balance exactly, 

c. For some computations in astronomy and space science, it is necessary to have 
an absolute time that is a continuous count of the number of time units from 
some arbitrary reference. The univarsally accepted standard is the Julian Day Cal- 
endar, a continuous count of the number of days since 12:00 noon on 1 January 
4713 B.C, This curious starting date was actually chosen in A,D, 1582 by consid- 
ering the cycle that is the least common multiple of the 28-year solar cycle (the 
interval required for all dates to recur on the same day of the week), the 19-year 
lunar cycle (the interval containing an integral number of lunar months), and the 
15-year indiction (the tax period introduced by the Roman emperor Constantine 
in A.D. 313). The year 4713 B.C. was the most recent date prior to 1582 when 
these cycles coincided, and it had the added advantage of predating the ecclesiasti- 
cally approved date of Creation, 4 October 4004 B.C. How long is the Julian day 
cycle, and when is the next year when all three of the cycles used in its creation 
will coincide? 



Solutions The least common multiple of 28, 19, and 15 is their product, since these numbers 
have no prime factors in common. 28 x 19 x 15 ^ 7980, so the next year the 
cycles coincide will be (-4713) + 7980, or A.D, 3267. 
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d. A clever computer algorithrn for converting calendar dates to Julian days was 
developed using FORTRAN integer arithmetic (H, F Fliegel and T. C, Van 
Flandarn, *'A Machine Algorithm for Processing Calendar Dates,'' Communica- 
tions of the ACM 11 [1968]: 657). In FORTRAN integer arithmetic, multiplica- 
tion and division are performed left to right in the order of occurrence, and the 
absolute value of each result is truncated to the next lower integer value after 
each operation, so that both 2/12 and -2/12 become 0, If / is the year, / the nu- 
meric order value of the month, and K the day of the month, then the 
algorithm is 

JD ^ " 32075 + 1461 * (/ + 4800 + (J-14)/12)/4 

+ 367 ^ {J -2-{J-U)l 12^:2)1 12 = 3 * ( (/ + 4900 + (/"14)/12)/100)/4. 

The calendar date 25 December 1981 is JD 2 444 964. Use a hand calculator and 
this algorithm to find the Julian dates of the launch of Explorer 1 (the first U.S. 
satellite placed into orbit), Greenwich Mean Time 1 February 1958 (Eastern 
Standard Time January 31, 1958), and the launch of the seventh Space Shuttle on 
18 June 1983 (carrying the first American female astronaut, Sally Ride, into orbit). 

Solution: For 1 February 1958, / - 1958,/ - 2, i^' ^ 1. 

JD ^ 1 - 32075 + 1461*(1958+4800+(2=14)/12)/4 

+ 367K2=2-(2-14)/12n2)/12 - 3*( (1958+4900+(2-14)/12)/l00)/4 

^ 1 32075 + 1461*6757/4 + 367*(in2)/12 3*(6S57/100)/4 

= 1 = 32075 + 2 467 994 + 367 -51^ 2 436 236 

For 18 June 1983, / ^ 1983, / ^ 6, ^ ^ 18, 

JD ^ 18 32075 + 1461*(1983+4800+(6--14)/12)/4 

+ 367*(6-2=(6-14)/12n2)/12 - 3*( (1983+4900+(6"14)/12)/100/4 

^ 18 32075 + 1461*6783/4 + 367*4/12 - 3*68/4 



^ 18 32075 + 2 477 490 + 122 - 51 2 445 504, 



A large number of satellites require ground processing of spacecraft sensor data 
to determine the spacecraft attitude (i.e., the spacecraft's orientation). Examples 
of sensors used are Sun sensors, Earth sensors, and star sensors. These sensors 
provide information, usually a measured angle, concerning the spacecraft pointing 
relative to a celestial body (e.g.. Sun, Earth, or star). 

Telemetry signals from these sensors are converted on the spacecraft to digital 
counts and transmitted to ground stations. The digital count representation of a 
sensor output can be easily converted to meaningful measurements and units on 
the ground. However, telemetry signals are frequently subject to random inter- 
ference, or "noise." To understand the meaning of noise, one has only to tune 
into a weak channel on a television set; the "snow" that is seen is a visual repre- 
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PROBLEM 12. 



sentation of noise present in an electronic signal. Noise consists of fandom sig- 
nals supefimposed on valid electronic signals ffom any electronic device. An 
in-depth understanding of the cause, effect, and reduction of noise is not neces- 
sary in the context of this problem. Ho^vever, it should be understood that noise 
can sufficiently corrupt any electronic signal to the extent that making use of, 
and properly interpreting, the true signal can be difficult. 

This problem applies to spacecraft instruments and sensors. A number of methods 
have been developed to smooth data and remove the effecis of noise. In the next 
pfoblem, we examine one such method, called the running average. 

Given an ordered set of numbers, Xj, j = 1 , 2, . _ M, a smoothed set of numbers can 
be found by averaging each number with the n preceding and the n following 
numbers. Symbolically, 



1 



^ 2n + 1 



2 X, 



where n is typically a small whole number (n s 5). 
For example, for the data in Table 2.1, if — 2, then 

XI - J (Xs + Xe + X, + Xg + X9) 



Solutions 



^ I (11 + 14 + 18 + 19 + 16) ^ 15.6. 

a. Compute the smoothed values of the data in Table 2.1 for n = 1 and n ^ 2. 
We present a computer program in BASIC, along with the run for this task. 



Table 2.1 



Spacecrslt Sensor Data 


Sample 


Unsmoothed 


Sample 


Unsmoothed 


No.(j) 


Value (Xj) 


No. (j) 


Value (X,) 


1 


2 


11 


20 


2 


7 


12 


20 


3 


10 


13 


18 


4 


6 


14 


19 


5 


11 


15 


20 


6 


14 


16 


20 


7 


18 


17 


17 


8 


19 


18 


19 


9 


16 


19 


18 


10 


17 


20 


16 
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3i& 


^ PRINT "lO" i l 


HTAB. 


10 1 


PRINT 




X(J) 






m 


END 
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18 


17*7 
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b. If the data points are plotted and joined by line segments, we get a graph 
demonstrating data fluctuations. Normally, wEiat is of interest is the underlying 
smooth curve (hence the term smooihini) for t . liis process. Compare the graphs o:of 
the ixnsmoothcd data and the smooth eddata feir « - 2, 

Solution! Fig. 2.5 shows the pJot of the unsmooiheci daca (solid line) and the data smooth^med 
With ^ = 2 (broken line). 
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c- This smoothing technique cannot be used blindly. It ispossible to disguise ^the 
true nature of the data by smoothing. Modify the prograniof part (a) to smoo ^ th 
the data in Table 2.2 and compare the graphs of the unsmoothed data and the 
smoothed data for some value of for instance, 5, How has the smoothing te^ch- 
nique disguised the true nature of the data? 

Solution* The following program incorporates a subroutine to smooth the data of Tabled. 2 for 
values of n from 1 through 5, The graph of the unsmoothed data and the 
smoothed values for n = 5 are displayed in Fig, 2,6, We see that the original Mata 
had a sinusoidal form, with the noise appearing as some slight departures frci^^ 
the smooth curve. The smoothed data are still sinusoidaljbut the amplitude h 
been dfastically reduced. It is evident, if the program with these data is run, tfciat 
each increase in n reduces the amplitude more than the previous n. 
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1110 Yt'^s V.IOO * <YCJV'V 0.0O5) lYT 
:; . J)^ W ; Y%>/ 100' ^ ^ 

;ii2o - print; Jix< j) ,Y(j) . 
1130 ne5<t; j "' . - : 

1140 FOR I - CM - N i ) -TO M 

1150 PRINT like I) ' > ' 

llfO NEXT .1 

1170 RETURN ; " 



Counts 



36.000 



^loi of imoQthgd (■ ) and unsmoothad data dGfived from Table 2, 



N ^ 5 for smoothed data 



0,000 



Sfsmplo Number 



100.000 



ERIC 



Fig. 2.6 

We and this chaptar with a problem that shows how a classical mathematical modal 
is modified so that it can be used to determine the period of rotation of a planet. 
The modification uses some of the principles of scientific accuracy discussed ear- 
lier in this chapter and illustrates the use of successive iterations to refine a com- 
putation, easily done by computer. 

The procedure we develop has been used to determine the rotation period of 
Saturn more accurately than earlier estimates by using observations of variations 
in the planet's radio emissions made by Voyager L Since this planet has neither a 
solid surface nor any distinctive atmospheric features comparable to Jupiter's 
Great Red Spot, what is computed is the period of rotation of th 2 magnetic field of 
the planet. Because of the complex nature of the radio emission data, we illus= 
trate the method by computing the rotation period of Jupiter rather than Saturn. 
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We begin with a question that is essentially the same as the familiar **How much 
time elapses between successive alignments of the hands of the clock?"' but that 
sets the stage for the actual problem we wish to solve, 

PROBLEM 13- ^ Op Jupiter rotates on its axis once every 9,92 hours, and its moon lo revolves around 
^ Jupiter once every 42,5 hours. What is the length of time between consecutive 
passages of lo over a particular spot on Jupiter? 

Solutions r Let Ri and Ri be the angular rotation and revolution rates for Jupiter 
^ and lo respectively. 



1 hen 



360 
9,92 



36,3 degrees/hour 



^^nd 



R, - 



360 

42J 



8,47 degrees/hour. 



X3n Fig. 2,7, lo moves from A to B while the point S on Jupiter makes a complete 
iCTevoIution and then goes on to S' to be under lo again. So we must find the time T 
^^uch that RjT = 360 = RiT. Using the values above for and Ru we get 
&the following: 

36,3 r - 360 - 8.47 T 
36,3 r - 8.47 r- 360 



(27.8) 360 

^_ 360 
27,8 



12,9 hours 



Solution: 
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'^Xo see how this classic problem might be altered, suppose we don't know Jupiter's 
Flotation period (or that we don't know it very accurately). As it approached 
J^upiter, the Voyager was able to make observations of the times at which Jupiter's 
C3reat Red Spot appeared in the center of the disc as viewed from Voyager. We 
^^^^ant to use these observations to determine Jupiter's rotation period, 

ho. Voyagefs trajectory as it approached Jupiter is illustrated in Fig. 2.8. Modify 
il he results of part (a) to find Jupiter's period if the Red Spot is observed to be in 
tl he center at time — 2 h 25 min ± 1 min, when Voyagef s distance from Jupiter is 

^ 7.70 X 10- km, and again at time 6 = 16 h 24 min ± 1 min, when Voyagers 
^Histance from Jupiter is Di = 4,72 >< 10^ km and Voyager has moved through an 
fUAngle a = 147.2^ with respect to Jupiter's center between these two observations, 

Ivf Voyagers trajectory were circular and if Jupiter's period is P, then the analysis of 
l»art (a) can be applied with = 360/P and RiT - a. We now have (360/ P)T - 
^S60 — a, which transforms into T ^ P ^ aP/360. However, the trajectory is not 
c circular, so we rnust take into account the different lengths of time it takes for 
Imght to travel from Jupiter to Voyager. Since light travels at a speed c (c - 3,00 x 



1 km/s), the corrected equation is T 



^ P + 



360^ 
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Fii^ 2.7 Pig, II 



For the data provided, we note that 

* 1 second. 

Since the uncertaincy (possible error) is one minute in both time measure ments, 
we see that in this ease there is no point in muking this correction. The eq iiation 

T ^ P + ~ ^ 1 + may be solvedforP: 

T 



1+ ^ 



360 



Substituting the dati 



^2 - h 



- (16.40 - 2,42)Mr 13,98 hours. 

360 360 

- 1.409, so 
13.98 ^ 



1.409 

9,92 hours 
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c. The observations of the radio emissions from Saturn were much more er ratic, 
A series of measurements tu ■ ^ ^ ^ h of the times of certain characteristic peak 
emissions was recorded, but it was not known how many complete rotations had 

actually occurred between each consecutive pair (Ji^u A)- The correction 

is not negligible in this case, and other corrections not discussed here also apply. 
Under these circumstances the equation relating T, P, and a becomes impossible 
to solve directly. However, the iterative approach illustrated next can be used, 
and by applying it to the case of part (b), we can show that the same solution is 
obtained. We let tjQ — ti + Since P is Jupher's period, ^an be considered a 
"corrected time'' in the sense that if Voyager*s revolution period around Jupiter 
matched Jupiter's rotation period, we would have s ^ 0, P = F, and tjc ^ h- 
Substituting T ^ ^ t\, and P ^ hQ ^ h for the first term on the right in the equa- 
^ b\ aP ^ aP 

aP 

Solving for ^e, we get ^e ^ £2 — Since this equation still contains P, the 

quantity we are seeking, we begin by making an initial guess (^,) at its value; this 
is used together with the known values of t2 and a to find a first estimate of fjc? 
which we may call r2ci' Next we let = ' = fj and repeat the evaluation of to 
get fjes; then Pj ^ faea " ^i- This proce:;s is rcr sated until the difference Pn = is 
less than 1 min — 0,02 h, the possible er r of our time observation. This pro- 
cedure is easily done by computer (and in this case may even be done using a 
hand calculator). Write a computer program to perform this iteration. 

Solution: We display a program written in BASIC. 



3LIST 

10 - JUPITER PERIOD ITERATION 

, lO/.^^'REW/J^wJOUNDlD TQ 3 DEC I HAL PL A 

. 130 PtO)' » T2 - Tl 

135 PklUt "M-ii HTA5 Ss PRINT "T 

^ . ' in >."-.ji ■ pj M),s ..^„>,r :^ .-J ■ ,7 - 

|gS^Y.m^>A;1000: » TECHS 2Y'<nJ ^' ¥% / 
' iTd'.iPRiNT^rffVHTAB 'Bs^-PRINT YtHI 

' lao .; v-if ^ jjABii.-c p cr > . p c# — ;i ) i < 

190 -;Ni'XT,,'M''^':x^;-;-4\i .-^ 
2Q0^*>RlNr'' ".AFTER " |M "ITIRATIQN 
. / Si JUPITER'S PERIOD IS FOUND 
' '-TO 8E ** fK(M) 

V ZIO END 



3RUN 

M 

1 

2 

e 

7 - 
8 



TO<M). 

isvoea- 
12. a5^ 

12*295 
12,36 " 
12.33a h 
12.3aS 







^ AFTER 8 ITERSTXONSf J\li!>trBt:'S'X}S'&iKf^: 
'■ PERIOD IB FOUND TO BE a';BZt^,;ffi^X'';<-W^- 
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Artistes concept of Pioneer 10 leaving the 
solar system, June 13, 1983, 
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Chapter Tliree 



Algebra is the language of quantitative science. As such, its methods and tech- 
niques can be found in most of fhe examples in this volume. The problems 
iselected for this chapter are those that do not also draw heavily on other 
mathematical areas. Several use the distance, time, and rate relationship, and 
others use direct and inverse variation. Some approximation techniques that are 
frequently used to solve otherwise intractable problems are also included. 

During 1982, the planets Jupiter and Saturn were in conjunction, so that they 
appeared very close to each other in the night sky. In the problem that follows, we 
see how frequently such an event happens, 

PROBLEM 1, The planets Earth, Jupiter, Saturn, and Uranus revolve around the Sun approxi- 
mately once every 1, 12, 30, and 84 years respectively. 

a* How often will Jupiter and Saturn appear close to each other in the night sky 
as seen from Earth? 

Solution! The time requirad must be a multiple of both 12 years and 30 years. This event will 
recur at intervals of the least common multiple of 12 and 30, or 60 years. 

b. How often will Jupiter, Saturn, and Uranus all appear in the same area in the 
night sky as seen from Earth? 

Solution: Now we need the least common multiple of 12, 30, and 84, which is 420 years. 

In addition to the electromagnetic radiation that we know as heat and lights the 
Sun continuously sends out charged particles known as the solar plasma (see also 
Chapter 10, Problem 1). From time to time, there is a strong burst of highly ener- 
getic particles called a solar flare from a small source in the Sun's atmosphere, 

PROBLEM 2. When a solar flare occurs on the Sun, it can send out a blast wave that travels 
through interplanetary space at a speed of 3 x 10^ km/h, 

a* How long would it take for such a solar flare blast wave to get to Earth where 
it could be detected by a satellite in orbit? (Recall from the preceding chapter 
that for a satellite close enough to Earth, we can use the Earth-Sun distance of 
1,5 X IQ- km for the satellite-Sun distance,) 



Solution: Since distance ~ speed x time, 

1,5 X lO^km - (3 X 10^ km/h) (time). 
L5 X 10^ 



Then lime 



3 X 10' 
or 50 h (about two days) 



h ^ 0,5 X 102 h. 
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b. When such a solar flare is detected, it is interesting to study the source. Since 
the Sun is rotating, we must determine how far the source has turned between the 
emission and the detection of a solar flare. Because the Sun is a dense gas rather 
than a solid body, it does not have a uniform rotation rate; on the average the Sun 
makes one complete revolution in 25,4 days. How many degrees would it rotate 
(on the average) during the time the blast wave traveled to the orbiting satellite? 



Solution: Since the Sun rotates 360^ in 25.4 days, it rotates 



360 
25.4 



7day, so 



the solar rotation rate — 14,2 ^/day 



14,2^ 



1 day 
24 hour 



PROBLEM 3. 



Solution: 



= 0.59 Vh. 

In 50 hours, the Sun rotates 50 h x = 29.5^. 

A scientific capsule was carried aloft by a rocket and released at the peak of the 
rocket's trajectory. The rocket had an average vertical speed of 920 km an hour 
from liftoff to release of the capsule. The capsule made a controlled descent with 
an average vertical speed of 390 km an hour and landed 68 minutes after the 
rocket was launched. Find the maximum height reached by the rocket. 

Let r be the time of ascent in hours. Then^K ^ the time of descent, and since the 

ou 

distances of ascent and descent are equal, 

920/^ 390 (1 = ^ 

= 442 - 390 / 
1310/ = 442 



/ - 0,34 h. 



The maximum height ^ distance of ascent = (920 km/h) (0.34 h) ^ 310 km, 

Italian scientist Galileo Galilei's introduction of the telescope for studying the 
heavens brought about a revolutionary change in astronomy. It is expected that a 
comparable leap in our ability to examine the universe will take place when the 
Hubble Space Telescope is launched into orbit by the Space Shuttle in 1986. 
Because the space telescope will be above Earth's atmosphere, it will be able to 
see much fainter objects than can now be seen by the best Harth^based telescopes. 
As we shall see in the next problem, this means that astronomers will soon be 
able to make observations and to compare them with the cosmological theories 
about the age and formation of the universe. 
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PROBLEM 4. a. The Space Telescope will be able to see stars and galaxies whose brightness is 

only — of the faintest objects now observable using ground-based telescopes. The 

bfightness of a point source such as a star varies inversely as the square of its 
distance from the observer. How much farther into the universe will the space 
telescope be able to see compared to ground-based telescopes? 

Solutions Let da be the distance from Earth of the faintest object visible to a ground-based 
telescope, and let Bq be the brightness of this object. Let d be the distance of an 

object of brightness ^ 5g. Since brightness varies inversely as the square of dis- 
tance, Bg ^~Bnd~Bc ^~.Then~ ^ di - ^ , so ^ SOdt.otd ^ 7A da. 
The Space Telescope will see about seven times farther. 

b* Because of the time it takes for light to travel from, distant stars and galaxies, 
we see them as they were some time ago— the photons that reach us from an object 
that is 1 parsec away were actually emitted 3,26 years ago (see Chapter 2, Prob- 
lem 12). The best ground-based telescopes can see objects about 10^ parsecs from 
our solar system. How long ago were the photons emitted that we now see when 
we observe such an object? 

Solution: Since 10^ parsecs - 3,26 x 10^ light-years, the photons were emitted 3,26 x 10^ 
years ago. 

c. When the Space Telescope begins its observations, how far back in time will it 
see stars and galaxies? 

Solution: Since it will see 7,1 times farther, it will see photons that were emitted 7,1 x 

3.26 X 10^ ^ 2.2 X 10^" years ago, (If, as suggested by cosmological theory, the 
age of the universe is between 10 and 20 billion years, the space telescope should 
enable us to see stars and galaxies in the earliest stages of formation,) 

Pioneer 10 was launched on 3 March 1972, It outlived and outperformed the fond- 
est dreams of its creators. Designed to last at least 21 months, it has continued 
well beyond the accomplishment of its mission. On 25 April 1983, its distance from 
Earth equaled that of Pluto, and the following June it crossed Neptune's orbit 
and left the solar system. (Although Pluto is normally the outermost planet in the 
solar system, it has a highly eccentric orbit and will be closer to the Sun than 
Neptune will be for the next seventeen years.) To add to its record of endurance, 
inost of Pioneer 2 OH instruments are still working, and Earth-based tracking 
stations were still receiving signals bearing information about the behavior of the 
Sun's extended atmosphere as of this writing. 
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PROBLEM 5, a. How long did Pioneer 10"% radio signals, traveling at the speed of light 

(3.00 X 10' km/s) take to reach Earth from the distance of Pluto in April 1983 
(4.58 K 10^ km)? 

Solution: time - ^listance _ 4.58 x 10^ 



speed 3.00 x 10' ' 
^ 1.53 X lO's 
1.53 K 10' 



3600 



h, or 4.25 hours 



b. What was Pioneer 10"% average speed, in km/h, if it traveled about 
4.58 X 10^ km between 3 March 1972 and 25 April 1983? 

Solutions From 3 March 1972 to 3 March 1983 there were 11 years, of which 2 were leap 
years, and from 3 March 1983 to 25 April there were another 53 days. The time 
for Pioneer 10 to travel that distance was (365 x 11) + 2 + 53 = 4070 days or 
4070 X 24 ^ 97 680 hours. 

Average speed - ^ip^ ^ km/h - 4.69 x 10' km/h 

(We note that the average spaed over this period is less than the average speed 
over the 21-month period of Problem 6 in Chapter 2.) 

The time required for an orbiting satellite to make one complete revolution of 
Earth is called its period. The length of the period depends on the location of the 
observer making the measurement. 

Suppose the observer is located far out in space and views the satellite against 
the background of fixed stars. The period measured in this manner is called the 
sWerea/ period of revolution, or the period in relation to the stars. Note that the 
rotation of Earth does not affect the sidereal period. Now suppose that the 
observer is standing on Earth's equator. A satellite is overhead in low Earth 
orbit moving directly eastward. When the satellite has made one complete transit 
of its orbit, it will not yet be overhead for the observer because the rotation of 
Earth will have carried the observer a distance eastward. The satellite must travel 
an additional distance to again be over the observer's head. The observer mea- 
sures the period of the satellite as the time elapsing between successive passes 
directly overhead. This period is called the synodic period of revolution, or the 
period between successive conjunctions, and it takes into account the 
rotation of Earth. 

Spacecraft usually orbit in the same easterly direction as Earth's rotation: this is 
called ELposigrade orbit. All U.S. manned spaceflights have been launched in posi- 
grade orbits to take advantage of the extra velocity given to the spacecraft by 
Earth's rotation. In this case, the synodic period is greater than the sidereal period. 

If the direction of orbiting is westerly, or opposite to Earth's rotation, the orbit is 
said to be retrograde. In this case, an Earth observer would meet the satellite 
before it made one compjete. revolution around Earth, and the synodic period 
would be less than the sic^retil period. 
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PROBLEM 6. In Chapter 9, Problem 7, we show t hat the sidereal period (in seconds) can be 

coniputed by the formula P ^ lirVa^/GM, where a is the average radius of orbit 
from the center of the body about which the satellite is in motion, G is the constant 
of universal gravitation, and M is the mass of the body about which the 
satellite orbits* 



a. Find the sidereal period of the High Energy Asironomy Observatory (HEAO) 
satellites, which have an average altitude above Earth of 430 km. The radius of 
the Earth averages 6370 km, and the value of the product GM for Earth is 

3,99 X lO^^mVsecl 

Solutions The radius of orbit is the sum of the radius of Earth and the average altitude 
of the satellite: 



a ^ 6370 km + 430 km ^ 6.80 x 10^ m. 
Then the sidereal period in seconds is 

p. 2 (3. 142) x^^^ 
- - V 3.99 H 10'^ 



P - (6,284) (6S0)^^ X 10= - 5580 s 
The sidereal period then is 93.0 minutes, or 1.55 hours. 

b. Compute the synodic period of the HEAO satellites, given that their orbits 
are posigrade. 



Solution! In Fig. 3T, let x be the position of the observer (assumed on the equator) when the 
satellite is directly overhead and let >^ be the observer's position one synodic 
period laiei , due to the rotation of Earth. 



O— * O Fig. 3a 
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If we can find the angular distance A. we shall be able to use it to find the synodic 
period. In one synodic period the observer traveled an angular distance >1, and 
the satellite traveled an angular distance 360^ + A, measuring the angular distance 
in degrees. The observer travels 360" in 24 hours, or 1" in 24/360 hours, so it 
takes the observer (24/360) (A) hours to travel the angular distance A. From part 
(a), the sidereal period is L55 hours. It takes the satellite L55/360 hours to 
travel r, and the time that elapses between successive viewings over the observer 

is therefore (^^) (360 + A) hours, 

24 1-51 

24 A ^ (L55)(360) + 1,55 A 
22.45^ ^ 558 

A = 24.9 degrees 

So the synodic period is {^^) (360 + 24.9) - 1,66 hours ^ 99,6 minutes. We 
observe that the synodic period is 6.6 minutes longer than the sidereal period. 

PROBLEM 7. The statement has been made that Newton's derivation of his inverse-square law of 
gravity from Kepler's third law is among the most important calculations ever 
performed in the history of science. Kepler's third law, based on observation 
rather than theory, states that the squares of the periods of any two planets are 
to each other as the cubes of their average distances from the Sun, Derive Newton's 
law from Kepler's law. 

Solution: If we represent the periods of any two planets by t and Tand their distances from the 
Sun by r and R, respectively, then 

or 



Assuming that we know the values of t and and substituting a constant C for 
the quantity p the equation can be reduced to 

T~ ^ CRK 

Thus if we know either Tor ^ for the second planet, we can solve for the unknown 
quantity. In this problem, however, we wish to use this equation to discover a 
new relationship, Newton's law of gravitation. For a body moving in a circular 
path, the acceleration toward the center is 
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Substituting in F ^ 



t = ^ — 



The velocity of the body in the circular orbit is 



So 

r R V T ) 

Because 7^ ^ we find by substitution in tha previous equation that 

= 4trvn 1 

F = - _ - X — 
C ^= 

That is, the force holding a planet in orbit falls off as the square of the distance R 
to the Sun. Newton expressed the vaiue of K and obtained his law of universal 
gravitation: 

c- GMm 

This law applies not only to the attraction between a planet and the Sun but also 
to the attraction between any two bodies, G is the constant of universal gravita- 
tion, M and m are the masses of the two bodies, and r is the distance between 
their centers of mass. 

In solving the next problem, two special techniques are needed. One is a fre- 
quently used approximation based on the fact that (1 + ^) (1 = ^) ^ 1 - x^. 
If ^ is small (for example, suppose ^ ^ 0.01), then is very much smaller (for 
^ = 0,01, ^ 0,0001), and in this case it is well within the liitiits of experimental 

error to use (1 + x) (1 - ^) ^ l^orj^^-^ ^ I - x. The other is the substitution 
of a single variable for the ratio of two other variable names. 



40 



Ajsgbra 



PROBLEM 8, 



Frofii Kepler's third law, we see thcit the fafther a planet is located from the 
Suiiit. ^he longer its period is. Suppose Earth orbits thne Sun in a circle of radius 
r(r^ IS x 10^ km) with a period r(r ^ 1 year), r^han any spacecraft SC (see 
.2) orbiting the Sun in the same plane but at swine greater distance (r + a) 
wlllli^^ve a period larger than T, and if it starts from a point on the extension of the 
Esfth— Sun line (as shown), it will gradually lag farth er and farther behind. 



Spacccra:-., 



0- 



Earth 



Sun 

-o 



FlgJ,2 



Howev^^er, the situation changes if a is sufficiently smmll, because then tha gravity of 
Earlli, ^ in the configuration shown, adds appreciably -1^ that of the Sun, For the 
force hacnolding the spacecraft in orbit to balance th^ cambined pull of the Earth and 
Sun,tfc^a spacacraft must move a bit faster, In fact^ thjere is a particular value for 
^solh^^at the speeding up of the spacecraft is just suffScient to allow it to keep up 
witliE^^arth. If that happf ns, then the spacecraft orbits the Sun in a circle of 
radiui (r + a), but with period riike Earth. What is «lie value of a that allows 
sucha^^ii orbit? 

Solution! Letwis^ , m^, nttE be the masses of the Sun, Earth, and spacecraft, respectively. For 
Earth's^ motion, we have, as in the foregoing problem, 



Of 



Forllie ^ motion of the spacecraft, similar analysis give s 

- ™r^^"^(r + a)^/^r+a^;Til r~j ^ n^..ir + a) (-j ) 
Canceiadng and substituting = (~ j from (1) =: 

anddivE-iding by G • nts' 



(1) 



We 1 



1 r + a _ 1 + (g_/r) 



ms fl* (r + af- 
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Now let « ^ =s and ^ ^ | ^ Notice that z is very small, since o is much less than r. 
After these substitutions, 

fl= (r + af ~~ H ' 
If we divide each denominator by and again use z = ^ , 

Although this equation contains only u and it still has no simple solution, so v/e 
now make two approximations in the second term on the left: 



1 



— ^ - (1 - z)2 - 1 ~ 22 + z2 i 1 ^ 2z: 



(1 + 4 



^ + (1 - 2z) ^ 1 + z, simplifying to 4 ^ 3z, or z^ ^ | 
Now the quantity u - 'h % lO^S so z^ ^ 10^^, or z ^ 10^^ 



Since z ^ -= , s ^ r k 10^^ ^ L-5 k 10^ km. 
r 



The position we have found in this problem is an equilibfium point of the Sun- 
Earth system, A similar analysis can be used to show the existence of another 
equilibrium point on the sunward side of Earth, and in fact there are five such 
equilibrium points in all for any two-body gravitational system. These are called 
Lagrangian points in honor of the mathematician who first proposed thei^ 
ence (Fig* 3,3). It has been suggested that two of the Lagrangian points ol liie 
Earth-Moon system should be considered as possible locations >r future 
space colonies. 



Center of mass of 
SumEifth lysiem 




Earth 



path of Eorth's orbit 
V 



The five Lagrangian points 



Fig. 3,3 
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The ISEE-^3 satellite (third International Sun-Earth Explorer), a joint vantura of 
NASA and the Europaan Space Agency, was launchad in August 1978 and 
placed in a ^'halo orbit'' around the Lagrangian equilibrium point batween Earth 
and tha Sun, In this orbit, it monitored the Sun's emissions that approach 
Earth— without tha interference that would result if the satalJite were actually at 
the aquilibrium point where its radio antenna would have to point directly at the 
Sun. After this mission was successfully completed in 1982, ISEE^3's orbit and 
direction ware ch/ nged to conduct an exploratory survey of Earth's magneto- 
tail. In December 1983, the satellite was redirected toward the comet Giacobini- 
Zinner and renamed International Cometary Explorer (ICE) in keeping with its 
new mission. It reached this comet in September 1985. 

PROBLEM 9, a. If Af is the mass of Earth, then the mass of the Moon is 0.012Af . The radii of Earth 
and the Moon are 6370 and 1740 km, respectively. Use these facts with Newton's 
law of universal gravitation to find the ratio of surface gravity on the Moon to 
surface gravity on Earth, 

Solution: If we place a mass m at the surface of Earth, then the gravitational attraction 
between the mass and Earth is 

p _ OMm 
- ^ (6370)2 ' 

Similarly, the attraction between the Moon and an equal mass m placed on its 
surface is 

G(0,012M)m 
" "(1740)- ' 

The ratio of to ^ is 

^ = Q'012 (6370)= 
E " (1740)2 ^ j 

^ 4.87 y 10^^ 
3.03 X 10^ 

^ 1 

That is, gravity at the surface of the Moon is g as great as gravity at the surface of 
Earth. 

b. If a man weighs 180 pounds on Earth, what would he weigh on the Moon? 
Solution: Weight on the Moon would be as follows: 



2 X 1801b ^ 301b 
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PROBLEM 10, In soma future space stations it is expected that artificial gravity will be created by 
rotation of all or part of the station. Gas jets or other propulsion devices can be 
used to control the rate of rotation of the station. As with the centrifuge, the 
rotation will produce a force against the astronaut that cannot be distinguished 
from gravity. If r is the distance of a point in the station from the center of rota- 
tion, then the velocity of the point for rotations a second is 

As noted in Problem 7, 

e = ^ , or y = Vor. 



Setting the two velocities equal. 



If r is given in meters^ then a is the acceleration in meters per second per second* 
By contfolling the values of r and N, any desired artificial gravity can be produced. 

a. Compute the rotational rate needed if the radius of the station is 30 m and a 
gravity equal to one-half the gravity of Earth is desired. (Use g — 980 cm per 
second per second, or 9,8 m per second per second,) 

Solutfonf a ^ I (9.8 m/s") ^ 4 J 




V07l57 
27r 



^ 0.063 

The rate of rotation must be 0.063 rotation per second or 60 x 0,063 - 3,8 rota- 
tions per minute. 
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b. Compute the needed rotational rate if the radius of the station is 150 m and 
Earth surface gravity is desired. 



Solution- 



= 0.04 

The rata of rotation must be 0.04 rotation per second or 2,4 rotations per minute. 

PROBLEM 11. The force of gravitation with which one body attracts another is inversely propo 
tional to the square of the distance between them. Consequently, the piill of the 
Moon on the oceans is greater on one side of Earth than on the other. This gravita* 
tional imbalance produces tides^ The Sun affects the tides similarly. Because the 
Sun exurts an enormously greater pull on Earth than the Moon doesj one might 
think that the Sun would influence the tides more than the Moon. Just the oppo- 
site is true. How can this be? 

Solutioni Let N be the point on Earth nearest the Moon and let F be the point on Earth 
farthest from the Moon, We shall assume that the tide-raising force of the Moon 
is in some sense measured by the difference in the pull of the Moon on unit masses 
located at N and F (see Fig. 3.4). If r is the distance from the center of the Moon 
to N and if is the diameter of Earth, then the forces acting at N and F are, 

respectively, g^d = — ^, , M being the mass of the Moon and G the uni- 

n (r + D^y 

versal gravitational constant. The difference between these two forces is the tide- 
raising force^ which we shall call TheUj 



2CMD, ( 1 + ^ ) 



Because ^ is very sman, this expression is approximately 



2GMD, 



Earth 

Moon 

O 



Fig. 3.4 ^ 
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Thus we would expect the tidal effect to be inversely proportional to the cube of 
the distance, whereas gravity is inversely proportional to the square of the dis- 
tance. Because the distance from Earth to the Sun is enofmously greater than the 
distance to the Moon, it is not surprising that the Moon provides the dominant 
tide-raising force. Local horizontal components of this force cause the tides to roll 
in and roll out (i.e., the horizontal movement of the water), 

b. From the foregoing, we can compare the tide-raising forces of the Moon and 
the Sun, If we use the subscript m for variables that apply to the Moon and s for 
those that apply to the Sun, the ratio 

The mass and distance of the Moon and Sun are as follows: 
- 73.5 X 10-^kg; M, - L99 x 10%g 
r„ ^ 3.S4 X lO^km; ^ L5 x lO^km 

Compute^. 



^ . Fxm 73.5 K 10=^ X (1.5 X 10^)3 248 

F[, 1.99 X 10^" X (3,84 x l{^^f 113 



So the tidal force exerted by the Moon is more than double that exerted by the Sun 
on the Earth. 
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A hypothetical space station with "electronic 
mall" capabilities, 
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Chapter Four 



Geometry is fundamental to space science, A multitude of activities, from the 
prediction of flight paths to the design of equipment, make use of geo- 
metric analysis. Geometry enters into many of the problems of the preceding 
and subsequent chapters. Most of the pfoblems in this chapter fall into three 
categories: those involving areas and volumes of plane and solid figures, those that 
use similarity, and those that use properties of circles or spheres. The Sun-Earth- 
Moon system happens to exhibit a striking geometric coincidence, which we exam- 
ine in the first problem, 

PROBLEM 1* To an observer on Earth, the Sun and the Moon subtend almost the same angle in 
the sky. The average angle is 0,52 degrees for the Moon and 0,53 degrees for the 
Sun, Depending on the particular location in its elliptic orbit, the Moon*s angle 
ranges between 0.49^ and 0.55^, whereas that of the Sun ranges between 0.52^ 
and 0,54°, This is why the Moon sometimes completely blocks the Sun, producing 
a total solar eclipse, 

a. If the mean lunar and solar distances are respectively 3,8 x 10^ km and 

1 .5 X 10^ km, what is the ratio of the solar diameter to the lunar diameter, and 

what is the ratio of the solar volume to the lunar volume? 




Solution: The geometry of the eclipse is illustfated in Fig. 4,1. Since the angle at 0 is the same 
for both the large and the small triangles and the triangles are isosceles, they 
must be similar. Letting and Rs denote the lunar and solar distances, 
respectively, and Dm and Ds the lunar and solar diameters, we have 



Rm 3,8 X 10' 



If Vm and V$ are the lunar and solar volumes, respectively^ 

Vu (mMDu/ir - [dJ ^390) 5.9X10. 
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b* Since the angle at 0 is very small, we can approximate the lunar or solar 
diameter by the arc length of the circle with radius or where this arc length 
subtends the angle at 0, Use the relation s ^ rd (6 in radians) to determine the 
actual values of Dm and D^. 



Solution! 0,52^ ^ 0.52^ x 



180^ 



0.0091 rad 



Dm ^ RsiiOmn) = 3,8 X 10-^ K 0,0091 ^ 3.5 x 10^ km 
0.53" — 0,0092 rad, by the same conversion just shown 

Ds ^ Rs (0.0092) ^ L5 x 10^ x 0,0092 ^ L4 x 10' km 
(Note: The reader may prefer to avoid the approximation by using the tangent 



function— that is, tan 



lunar radius 



solar radius 



2 lunar distance solar distance 
to two significant digits the result is the same,) 



note, however, that 



PROBLEM 2. 



All the energy to meet needs on Earth, whether the energy is natural or synthetic, 
ultimately comes or has come from the Sun In the form of electromagnetic radi- 
ation. There has been much interest recently in using this radiant source of energy 
directly to supplement or supplant the existing power sources. Further, since our 
Sun is but one of many stars, it is of interest to compare its energy output with that 
of other celestial objects. 

O a measure of the total energy radiated by the Sun received at a unit area of 
the Earth's surface is called the solar consiant (where radiation is summed over all 
wavelengths of the electromagnetic spectrum). 

A radiometer flown on the Solar ^daximum Ndission (SNINI) is able to measure 
accurately the intensity of solar radiation. SMM is a satellite in orbit around 
Earth at low altitude, and its measurements can be used to provide a good estimate 
of the solar constant. 

The radiometer on SMM admits solar radiation through a small aperture whose 
area is 0.50 cm-, and it measures the rate of entrance of this radiation accurately. 
The spacecraft attitude (pointing direction) is controlled so that the entrance aper- 
ture is perpendicular to the line of sight between SMM and the Sun. 

a# Over one observation period, radiation entered the radiometer at the rate of 
0.069 watts. What is the value of the solar constant, S, as determined by this 
observation? (Use an extra significant digit in the answer, since this quantity will 
be used in subsequent calculations.) 



Solution: 



0.069 watts 
0.50 cm' 



^ 0.138 watts/cm^ 
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b. It is generally assumed that the Sun emits radiation uniformly in all directions. 
If this is true, calculate the total rate of energy radiation by the Sun. 

Solutioni Since the radiation energy rate measurement containb only two significant digits, we 
can use the Earth-Sun distance of 1.5 x 10^ km as SMM's distance from the Sun 
(see Chapter 2, Problem 8). If the Sun emits uniformly in all directions, the total 
rate of energy radiation from the Sun is the product of the solar constant and the 
area of the sphere with radius L5 x 10^ km, or 1.5 x 10-^ cm. 

Letting P ^ total rate of energy radiation from the Sun, 
P - (S)(4w2) 

^ (0,138 watts/cm=)(4^)(1.5 x 10^' cm)' 
- 3.9 X 10^^ watts, 

c. The foregoing are typical values. Variations of approximately 0,05 percent have 
been observed at other times. How much do such variations affect S and Pi 

Solutioni M ^ 0.05 x 10^= x 0,138 ^ 6.9 x 10^^^ watts/cm' 

AP ^ 0,05 X 10^2 X 3.9 x 10=^ ^ 2.0 x lO^Svatts 

(Note: These variations occur on a short time scale (day to day) and are thought to 
average to zero over a long time scale. A 0.05-percent systematic variation in 
solar radiation over a time scale of years could produce significant climate changes 
on Earth.) 

d* li 1981, SMM lost pointing accuracy because of a component failure on the 
spacecraft. Suppose that the orientation of the spacecraft changed so that the 
line perpendicular (the normal) to the entrance aperture made an angle of 30^ with 
respect to the Sun-SMM line, rather than being parallel to it. By how much 
would the radiation entering the radiometer be changed? 

Solutions For simplicity, let us assume the aperture is a square, ABCD (see Fig, 4.2), with side 
length a, where gr — 0,50 cm^. Looking at this square edge-on with AD as the 
tilted edge, if DE is parallel to the direction of solar radiation incidence and AE is 
perpendicular to DE, the aperture is effectively a rectangle whose dimensions 
are AB and AE. We label the angles a, y, Sm shown. Since (Z.a, /^p) and 
i^y, ^6) are complementary pairs of angles, and since ^/B = z^y, we have 
^8 ^ d.a ^ 30^, so 4ADE is similar to the standard 30*-60*-90* triangle, and the 

ratios ^^p,^^, and are equal, giving AE ^ ^AD ^ 0.866a. The area of 

the effective aperture is therefore (0,866 a) (a) ^ 0,866 «r. In other words, the 
radiometer will register only 87 percent of what it did before losing 
pointing control, - 

A 
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Fig. 4.2 
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PROBLEM 3. 



Solution: 



(We observe that the result holds for apertures that are not square,) 

As obsefved in part (c) of the foregoing problem, one of the interesting outcomes 
of modern advances in the precision with which it is now possible to make mea^ 
ourements of the solar constant is that this quantity is in fact not really a constant! 

Solar cells convert the energy of sunlight directly into electrical energy. For each 
square centimeter of solar cell in direct overhead sunlight, about O.Ol watt of 
electrical power is available. A solar cell in the shape of a regular hexagon is 
required to deliver 15 watts. Find the minimum length of a side. 

The total area required is 15 watts/O.Ol watt per square centimeter, or 1500 square 
centimeters. The regular hexagon can be partitioned into six congruent equi- 
lateral triangles, each with an area of 1500/6 ^ 250 square centimeters (see 
Fig, 4,3), 




Fig. 4.3 



The area A of any equilateral triangle with side s may be expressed 
A ^ I (base) (altitude) ^ ^ ^ Solving for s. we have 



^ Ua ^ /4(250) cm=^^^ 
^ Vx^= V 1,73 



8 cm" = 24 cm. 



PROBLEM 4, Solar cells are made in various shapes to use most of the lateral area of satellites, A 
certain circular solar cell with radius r will produce 5 watts. Two equivalent solar 
cells are made, one being a square with side s and the other an equilateral triangle 
with side p. Find r in terms of p and also in terms of j. 

Solution: For the solar cells to have equivalent outputs, their areas must be equal. Thus for the 
circle and square, we have 



For the circle and equilateral triangle, we have 



= 0.564 J. 
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= 0.371 p. 
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PROBLEM S, The largest component of the prelaunch Space Shuttle configufation is the external 
tank, which serves as the "gas tank'^'for the Orbiter^ — ^it contains the propellants 
used by the main engines. Approximately 8.5 minutes after launch, when most of 
the propellants have been used, the external tank is jettisoned. It is the only 
major part of the Space Shuttle that is not reused. 

Fig 4.4(a) shows the launch configuration with the back view of the Orbiter , and 
Fig. 4.4(b) shows the side view. Fig, 4, 5(a), (b), and (c) show the liquid hydrogen 
tank, the intertank, and the liquid oxygen tank, respectively. The intertank serves 
as a mechanical connector between the liquid oxygen and liquid hydrogen tanks, 
and contains the upper dome of the liquid hydrogen tank and the lower dome of 
the liquid oxygen tank. 



a. Using the dimensions provided in the diagrams, estimate the volume of each of 
the tanks by dividing the tanks into shapes whose volumes are easy to compute. 




(a) 

Liquid Hydrogen Tank 



Nose Cap' 




Fip. 4,4 



Liquid Oxygen Tank 



Fig. 4.S 
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(a) 



(c) 



Solutions (This is one possible solution.) The liquid hydrogen tank has the shape of a cylinder 
with ellipsoidal caps on each end. Since the formula for the volume of a hemi- 
sphere is better known, let us approximate the domes as hemispheres. Now the 
total length of the tank is given as 29.6 and the diameter as 8.4 m; our 
approximation, then, consists of two hemispheres (or a single sphere) of radius 
4.2 m and a cylinder of radius 4,2 m and length 29,6 = 8,4 — 21,2 m. The resulting 
volume estimate is 



ho 
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= 310 -+^ 1170, or 1480 

The liquid oxygen tank can be approximated by joining a hemisphere of radius 
4.2 m, a cylinder of radius 4.2 m and length about 4 m, and a cone with base radius 
4.2 m and height 8, 1 m. (This should probably underestimate the volume, since 
the tapered section is larger than a cone,) Using this dissection, we find that the 
volume estimate is 



1 . . 1 
3 ' ' 3 

7 



- (4.2)3 + IT (4.2)2 (4) + i TT (4,2)^ (S. 1) 



^ 155 + 222+ 150 
^ 527 m\ 

b* The actual volumes of the hydrogen and oxygen tanks, respectively, to the 
nearest m-^ are 1450 m' and 541 m\ What are the absolute and relative errors of 
the estimates in (a)? (See Chapter 2 for a discussion of these errors.) 

Solution: For the hydrogen tank: 

Absolute error ^ |estimate - true value] ^ |1480 - 1450j - 30 

o i absolute error 30 ^- 

Relative error = ^^^^ .^^^ x 100% - — x 100% - 2. 1 % 

For the oxygen tank: 

Absolute error - |527 - 54l| ^ 14 m^ 

Relative error ^J^^ 100% ^ 2.6% 



c* The outside of the external tank is covered with a multilayered thermal protec- 
tive coating to withstand the extreme temperature variations expected during 
prelaunch, launch, and early flight. Although there are variations in the exact type 
of material and the thickness at various locations on the tank, the average thick- 
ness is 2.5 cm. Estimate the total volume of the insulation material on the tank, 
assuming a uniform thickness of 2.5 cm. 

Solution: A simple way to get such an estimate is to model the external tank as three sections: 
the lowest section is approximately a hemisphere of radius 4.2 m; the middle 
section is an open cylinder of radius 4.2 m and height (47.0 - 4.2 — 8.1) ^ 34.7 m; 
the top part is approximately a cone of base radius 4.2 m and height 8.1m. The 
volume of insulation is then close to the product of the surface area of this figure 
and the thickness 2. 5 cm, or 0 025 m. 

The surface area of an open hemisphere of radius r is 2w^; the lateral area of a 
cyUnder of radius r and length h is lirrh; the lateral area of a cone of radius r and 
slant height s is wj— in this case we know the vertical he ight h rather than the slant 
height, but s, h, and r are related by s^^r^ + hf , or s - VpT^ (see Fie 4 6) 
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The surface area, then, is 



2^ (4.2)2 + 2^ (4,2) (34.7) + tt (4.2)= \/(4,2)2 + (8, l)^ 
^ 111 + 916+ 120 
^ 1147 ml 

The volume, then, is (1147) (.025) ^ 29 mi 





F|g< 4.6 



Fig. 4.7 
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PROBLEM 6. A spacecraft is at P, at an altitude h above Earth's surface, as pictured in Fig, 4.7, 
The distanca to the horizon is and r is the radius of Earth. 

a. Derive an equation for d in terms of r and h . 

Solutloni Because PA is tangent to the circle at A, angle PAO is a right angle. Then 

r- + d^^{r +hf 
^ (r + hf = 

d^Vlrh +/i2, 

b* The satellite Atmospheric Explorer 3 (AE-3) has an elliptic orbit with apogee 
height 4300 km and perigee height 150 km. Find the distance from AE-3 to the 
horizon at apogee and perigee. The radius of Earth , to two significant figures, 
is 6400 km. 



Solution: At apogee: 

d « V2(6400) (4300) + (4300)^ 



At perigee: 



d-V2(6400) (150) + (150)2 



^ 100V5504 + 1849 
^8600 km 



^ 10Vl9200 + 225 
- 1400 km 
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is small comparad to the forniula for d f g und in part (a) can be siin- 
plif ied by dropping the A' term, leaving d ^ V^^, Redo tha calculations of part 
(a) using the limplified formula and compute the relative errors of 
these approximations* 

Solutioni At apogee; 

dopprq,. ^ V2(6400) (4300) ^ 7400 km 



1 8600-7400 ^ 

relative error ^ o^n Fr~ -0^14^ 14% 

At perigee: 

^.ppr... ^ V2(6400)(150) - 1400 km 

This agrees with the previous result, and the relative error is 0. 

d. For what range of values of h is the approximation d = V2r/i accurate to two 
significant digits? 



S-olution^ We need to know the range of values of h that satisfy the following condition: 

V2rA +/i^< V2^ + 0.0lVi^ 



2rh +h'<iim)H2rh)^2rh (1,0201) 
h^KOmrh 



h<OMr 



For r ^ 6400 km, we need h < 256 km. 



From the foregoing, we see that under certain conditions it is possible to substitute 
a simple formula for a complicated one without affecting the results. Great care 
must be taken, of course, to ensure that the conditions needed for such sim- 
plification are in fact satisfied. Another useful result based on two such approxi- 
mations is developed in Problem 8. But first we consider the basic geometry of 
photographic scale. 



55 



61 



PROBLEM 7. In Fig, 4,8, the flight path of an airplane or satellite carrying a camera with its lens at 
C is shown by the arrow. The camera is at a height H above the ground and has 
focai length/, 




a. If PQ is the image on the film of line AB on the ground, find the scale of the 

. PQ . . 

picture, the ratio assuming the picture is taken vertically (PA is perpendic- 

ularto both the film and the ground). 



Solution: Since triangles ABC and PQC are similar. 



AB H' 



h. If/and H are in the same units, the ratio ^ is called a 1-1 scale factor, Deter- 

ti 

mine the 1-1 scale factor for a photograph taken at a height of 30 km with a 
camera having a 150 mm focal length. 



Solutioni 



X 

H 



150 X 10^ 
30 X lO^ 



^ 5 K 10^ 



c. If the photograph of part (a) shows an image of a straight road that measures 
1,25 mm on the film, how long is the actual road? 

Solutions Let L be the actual length of the road in meters. 



image length _ L25 k IQ^^ ^ 
actual length L 



5 X 10^^ 



1 7^ X 1 n^'^ 
^ = gt in" m = 0-25 X 10' m - 250 m 

5 X lu 



d» With current technology, it is possible to make measurements on photographs 
so the nearest micron (10^^ m). What is the smallest actual length whose image 
can be measured on the photograph of part (b)? (This is called the resolution of 
the photograph,) 



Geometry 



Solutioni Let S be the smallest actual length in meters. Then 



10^ 



^5 X 10=% 



so 



5 K lO^^S^ 10^ 



and 



The curvature of Earth can introduce distortion in a photograph. We need some 
new tefminology to discuss the correction for such distortion. The point on Earth 
vertically below the camera is called the nadir and its image on the film is called the 
photograph nadir point. Because the photograph Is flat, the image of any point 
except the nadir will be closer to the photograph nadir point than it would be if 
Earth were also flat. In the next problem, we develop a formula to correct for this. 



PROBLEM 8, The geometry of the photographic correction for Earth's curvature is shown in 

Fig. 4.9. The image of the point P is a distance r from the photograph nadir point 
Q,/is the focal length of the camera, and H is its height when the picture was 
taken. In order to get a ''corrected*' picture, we need to place the image at P% in 
the plane of the tangent to the nadir, N, This means we need to compute Ar so that 
the corrected image is a distance r + Ar from Q. 




Solution: Let h be the vertical displacement of P' with respect to P and let x be the horizontal 
displacement of P' (also of P) with respect to N. We see from the diagram in Fig. 4.9 
that X and h are related. If Tis the foot of the perpendicular from P to NE, 
where E is the center of Earth, then APTE is a right triangle with PE - R, 
FT^x, and TE ^R-h. By the Pythagorean theorem, R^'^x^+ (R - hf, so 
R^^x^ + R^~ 2Rh + ^^ giving ^ 2Rh -h^. Since h is very small compared to 
we shall use the approximation x^ ^ 2Rh. - ^ 
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There are two pairs of similar triangles in the diagram. 
For the inner pair, 

f _H±h 



or 



or 



For the outer pair, 



or 



or 



fx ^(H ^ h)r. 



f 

r + 4r X 



fx - H{r + 4r), 



fx^Hr + HLr, 



Comparing the two expressions for fx, we see that hr — H4r, Then ^ 

ri 

r X 

' Now we need another approximation in order to eliminate x. In the relation 

fx " Hr -¥ hr, since h is small compared to H, we have fx ^ Hr, so a " = , 

r H^r^ Hr^ 
Making this last substitution, we have Ar — ^j^^ ' ^2 = fRf^ ' 

b* Find the correction Ar and the resulting r + Ar for a photographic image taken 
at a height of 92 km with a camera having a focal length of 132 mm if r measures 
65,24 mm. Recall that R - 6400 km. 

Solution: Since H and R are in km and r and /are in mm, if we do no unit conversions, we shall 
be computing Ar in mm. 

. (92) (65.24)= ^ ^ ^ 
r&Ar^ 65,24 + 0, 11 ^ 65,35 mm 
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The final problems in this chapter deal with soma aspects of planning for the future 
construction of such large commercial space structures as the antenna system in 
the illustration on page 47. It is planned that the materials for the antenna system 
will be carried up and the actual construction done in orbit. This frees the con- 
struction of two considerations: (1) the rigidity that would be required for such a 
structure to break away from Earth's gravity and (2) the strength needed to 
survive transportation into orbit intact. It is, of course, desirable to keep to a 
minimum the number of trips needed to transport the components, and consid^ 
arable effort has gone into the development of materials that are strong and light- 
weight and that maintain their properties over a wide range of temperatures. Let 
us see how successful the effort to minimize the number of trips has been. 

PROBLEM 9, The Space Shuttle can carry 29 500 kg of payload into orbit in a cargo bay that is 
basically a cylinder having a length of 18.3 m and a diameter of 4.6 m. The 
structure in the illustration has 91 antennas, each a paraboloidal cap 20 m in 
diameter and 2 m deep. The material for the antennas is a knitted metallic mesh 
weighing 60 glw?. 

i 

The plan for the truss assembly that holds the antennas is shown in Fig. 4, 10. A 
promising material for the columns is graphite-epoxy, which combines excellent 
strength and stiffness with light weight, having a density of 1522 kg/m^. The truss 
assembly shown has 252 copies of the basic repeating element, with each repeat- 
ing element consisting of a tetrahedron having nine complete struts as shown in 
Fig, 4.10(c). The struts themselves are hollow columns 10.4 m long with radius 
3.8 cm and thickness 0.57 mm as shown in Fig. 11, 




Fig, 4^ 1 1 



Fig, 4.12 



a* How many Shuttle trips wouldbe necessary to get the weight of these elema xits 
(the rrietallic mesh for the antennas and the columns for the struts) into orbit? 

Solutlun: We sliall get an approximation for the weight of metallic mesh by treating the 
antennas as though they were eircb of radius 10 m.(In Chapter 10, *'CalcuIus_ 
we sfeall get a more accurate result.) If we have 91 clreles of radius 10 m, the total 
area x^ill be (91) (7r)(10)^ m^, andttie total weight of metallic mesh will be 
(91) C^) (10)^ (60) g ^ L7 K 10^^ k|-The volume of material in each column can toa 
approximated by (2w) (t) (€) ^ 2if(3J x 10^') (0 J7 x 10^^^) (10.4) m\ so the 
total weight ofthe columns is (252)(9)(eJ8) (3. 8)(0j7)(10.4)(10^')(1522) kg = 
4.9 K 10' kg. The total weight of tlie^e niaterials is (1.7+ 4.9) K 10' kg ^ 6.6 x 1 O^^ 

kg. ^V^ " 0,224, or aboiil22 percent of the Shuttle's weight capacity, 

b* We see that the Shuttle can easily carry this weight on a single trip. Now we 
must consider volume: will the niaterials fit in the available space? Assume thac 
the rrxetallic mesh is 7,5 mm thick iind sufficiently flexible to pack into any shap«. 

Solutioni= The cargo bay's cylindrical volunie is given by ^^^L ^ it(2, 3)^(18. 3) m^ ^ 3.0 >^ 10^ 
m^. ^\^e have already found that tho total area of metallic mesh is (91)(1007t) m^ 
^2.9 X IQ-* m% so the total volumeof mesh is (2,9 x 10^)(7.5 x 10^') m'' 
^2.2 X 10^ m\ 

This 1 eaves (3,0 - 2.2) x 10^ ni - =BOm^'for the columns (and all the remaining 
hardN^are needed for assembly, Nvhich we are ignoring here). 

Since the columns are 10.4 m long, they eannot be placed end to end in the 18.3 m 
long oargo bay. We must considerhow to stack them most efficiently, If we 
consider the cross section of the slaek, w^e see that we need to find the most 
efficient way to pack circles in a pliine. It is intuitively reasonable (although the 
proof is far from simple) that themaximiam efficiency is achieved when each circle 
is tangent to six others, as illustrated in Fig. 4, 12. From the diagram, we see that 
each hexagon has sides of length^r, where r is the radius of the circle, and there- 
fore Ixas area 6(1/2) (2r) (V3r) ^6V3r^ Also, each hexagon contains three coin^ 
plete circles whose total area is 3(irH), So the fraction of area occupied by circles is 

-^^--^-0.907. 
6V3r- 2V3 
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PROBLEM 10, 



Since the fraction of space occupied at the boundary will be smaller than this, let 
us assume t&at the columns pack into the cargo bay solhat 12 percent of tlia space 
is empty. Tbriis means that the columns will pack intoocross-sectional area equal 
to 1/0 JS of ^ their total cross-sectional area. We recalllhat there were 252 >s 9 ^ 
2268 colurncris, each having radius 3,8 cm and thereforea cross-sectional area of 
If (0,038)5 nm- = 4,5 x 10^^^ ml The total cross-sectionalarea needed then is 
2268 >< 4,5 >^ 10^^ 



' mr 11,7 ml Since the columns areIO.4 m long, they will 



occupy (IJ,^) (10.4) m'^ L2 x lO^ of space. However, there was only SO m^^of 
space left af-^er calculating the volume of the metalUcmesh, so it will take more 
than one tri^ to handle the volume, even though the might is not a problem. Our 
success in reducing the weight now places the focus ofour attention on volume. 

In order to Sit more columns into a smaller space, tjodesigners realized that they 
should investigate the possibility of tapering the coluniiisand then **nasting'* 
them for traaisportation, like a stack of paper cups. Fii.4,13 illustrates the idea. 
Under this s«-«heme, each column would be made of t^vo tapered half-coIurT^ns, 
with their wSde^ openings joined; half=columns couldta be nested for stowage in 
the cargo ba^. Tapered columns have been developed iind tested for strength. If 
n is the raciimis of the smaller end and r2 the radius of tlielarser end. tests sliowed 

that an optirsium taper ratio is^ ^ 0.41 and that suchntapered column is actu- 
ally strongei— ; it can carry about 30 percent more load Wore buckling than an 
untapered cc^lumn of the same weight. 




(a) Assembled ColMumn 

Fig, 4* 13. Tap€r^i«d Column Concept 



(b) Nestyd HalfCQiumn Elements 



a. If the mea^^n radius is to be 3,8 cm as before, and - ^0.4, find the valuer of r, 
and r^. 



n + 



Solution I We have ^ 



2 ^ ^ 3,8 and r, = 0.4 Clearing the fraciionand substituting, we get 

0,4 r, + ^ 7.6 
1,4 fa ^ 7,6 

r2 - ^ 5,4 cm 



EKLC 



r, = (0.4) (5.4) = 2.2 c 
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b. Fig, 4, 14(a) and (b) display the georrtstry of the tube nesting, where 



and A is the tube-nesting separ^tt^ion. 



^2 = 2 ' r^. € is the length ofii half-coiufni 

Show that 4 ^ ^ ^ and find an expression in terms of i and 4 for the nunu^ber of 

half-columns that %vill fit inioone stack tK^e length of the Space Shuttle car|^^ 
bay. 

Solutions In Fig. 4.14(b), if we insert the horizontM.! line shown and letter sorne key po-cints as 
indicated (Fig, 4,15), AABE' and 4BCP are similaMo AB/BC ^ BE7CP^ , We 
have AB = 4, BC ^ €, CD ^ anc^ we shall approximate BE' by ^ = /, 

"Then, with this approximation and the pr^oportion above, 4 ± — — . 

r. - ri 

From Fig, 4, 14(a), we haveone half-CQju^Tin of length f on the left end, in ^hdich 

flB 3 ^ f ^ 

we nest additional half-columns, where- n = INT ' 



(INTsigriifiQias the 

greater integer which is lessthan or equaZl to the number in the square br^^ka^et.) 
Now the number of half-columns that \vil 1 fit into one stack is N — 1 + // ^ C 1 + 

INT ri8-3 - r 



For the truss assembly of Problem 9^ 4 . eterniine the volume occupied by ttf he 
strut columns if they are madeof half-col ximns as described hare and neste^ for 
stowing in the cargo bay. 

Solution: V/e have € ^ half-column length =(1/3)^10.4) ^ 5.2m 



4 ^ 



i¥ ^ 1 + INT 



11,3^,2 
9 % 10=3 



(Q.57i< iC3^^) (5.2) .^.^ 

* 1 4— INT [145,6]^ 1 + 145 ^ 146. 



We had a total of 2268 columns, or 4536 laalf-columns, so this means there ^Vt^ill be 

[45361 



INT 



146 



]- 



31 stacks, and one additional shorterstack. 



Each stack is 18.3 m long Colthough one s=£ack will be shorter) and has a t^^I^mjs of 
3 A cm, so its volume is ^(0.054)^(18.3) ^ 0,17 m\ The total volume of th^^32 
stacks is a little less than 32 % 0,17 - SA m^. 



— - 5 4 

By the analysis in the last partof Problenm 9, these stacks will take ^ 

€.2 of space in the cargo bay, and now the materials for the truss asseni^l^^^y and 
The antenna ''dishes'' can all be transpoftaed in a single Shuttle trip. 
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Space Shuttle locker tray containing 18 stan- 
dard orbital flight test menu meals for two 
crew members. 
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robability theory and statistical techniques niake important c^n_ tributions to 
the space program. In this chaptef we examine the role of probability in 
■ menu planning, in some aspects of the transmission and coding of spacecraft 
observations, and in the control of equipment reliability. Some elefii^sntary 
examples of the use of statistics are illustrated in the fina! two pfObl^:^rTis = 

PROBLEM 1, The early manned spaceflights revealed much about the body's response to pro- 
longed weightlessness. An interesting and varied foodiupply ^vas tJiiss needed to 
guard against a loss of appetite in the face of what was learned. The F^^od supply 
for the crew of the Space Shuttle is carefully planned tocompenitt^ tor the high 
fenergy requirements (averaging 3000 calories per person per day) of ^workmg in a 
frictionless environment and the body's tendency to lose essential rtii^ierais (such 
as potassium, calcium, and nitrogen) in microgravity. The Space Shumtle food and 
beverage list contains more than a hundred individual items, A typjca.1 day's 
menu might be the following: 

Meai I Meai II __MeQlJlI 

Peaches Frankfurters Shriinp c^^cktail 

Beef patty Turkey tetraizini Beef steBBc 

Scrambled eggs Bread (2) Rice pilaff 

Bran flakes Bananas BrocCOlJau gratin 

Cocoa Almond crunch bar Fruit cocl^tail 

Orange drink Apple drink (2) Butterscb i:ch pudding 

Grape dn 

In general, each meal III contains a main dish, a vegetable^ and two dl^esserts, with 
an appetizer included every other day. The food list contains 10 it^riis classified 
as main dishes (M), 8 vegetable dishes (V), 13 desserts (D), and 3 ^pp-^tizers (A). 
How many different menu combinations are possible in each of th^ first six days 
of flighty assuming no dish is repeated? 

Solution^ The number of choices is tabulated below: 



Day 



M 



DI 



D2 



Hmber of cof^ibinations 



1 


3 


10 


8 


13 


12 


37 4^0 


2 




9 


7 


11 


10 


6 930 


3 


2 


8 


6 


9 


S 


6 9m2 


4 




7 


5 


7 


6 


1 4'^0 


5 


1 


6 


4 


5 


4 


4S0 


6 




5 


3 


3 


2 


»0 
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PROBLEM 2. The electronic telemetry system aboard a spacecraft transmits the data of spacecraft 
motion in the x, y, and z directions. The system consists of thrje motion sensors, 
a signal conditioner, and a transmitter. The probability of failure for each motion 
sensor and for the signal conditioner is 0=0001. The probability of failure for the 
transmitter is 0.001, Assuming that component failures are independent events 
and that the failure of any component will render the telemetry system inopera- 
tive, compute the probability of a spacecraft telemetry success, 

Solutioii: The probability of success is equal to 1 minus the probability of failure. Therefore, 
the probability of success for each sensor and the signal conditioner is 

P ^ 1 - 0,0001 ^ 0.9999. 

Similarly, the probability of success for the transmitter is 

P ^ 1 - 0.001 ^ 0,999. 

The probability of success for the telemetry system is the product of probabilities 
of success for each component; that is, 

P - (0.9999)^(0,999) ^ 0.9986, 

The signals transmitted by a spacecraft telemetry system are in the form of pulses 
imposed on a radio beam, which can be interpreted as binary digits. For exam- 
ple, the signal fragment . . . ^TLTn — ^ , . , will be read as . , , 0101 10 . . , , since the 
presence of a pulse is read as 1 and the absence of a pulse as 0. Each possible 
representation of a 0 or a 1 is called a "bit,'' 

For a variety of reasons, equipment errors can cause a 0 to be transmitted instead 
of a I5 or vice versa. As a result, error-detecting codes have been developed to 
improve data reliability. All such codes are based on transmitting extra bits that 
can be used to determine whether errors are present and even, for the more 
sophisticated codes, where the errors are. Transmitting these extra bits, however, 
means that fewer message-carrying bits can be sent in a given unit of time, and so 
transmission reliability must be traded against transmission efficiency. Probability 
theory plays an important part in weighing the trade-offs. 

PROBLEM 3. a) The telemetry system of a certain spacecraft has a probability of 1 percent of 
transmitting an erroneous bit. One way to increase data reliability would be to 
repeat each message bit three times. For example, . . . 010110 , , , would become 
. , , 000111000111111000 . , . , if no errors occur. If it is decided to interpret any 
of the triplets 000, 001, 010, or 100 as 0 and any of the triplets Oil, 101, 110, or 111 
as 1, find the probability of error in the interpretation of a message bit, 
assuming the transmission of each bit is independent. 

Solutloni A message bit will be interpreted erroneously if two or three errors have occurred in 
the triplet. 

F(2 errors) ^ (0.01)^ (0.99) ^ 0,000297 
P(3 errors) ^ (0.01)^ - 0,000001 



P(2 or 3 errors) ^ 0,000298 ^ 0,0003 

^^^^ 66 
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We see that we can reduce the probability of a transmission error in a single bit 
from 1 percent to 0.03 percent, but at a cost of sending three times as many bits 
as are actually needed for the message. To put it a different way, the desired 
message would be sent one-third as quickly. 

b* More efficient error detection can be done with pariiy coding. In this method, 
a '^parity bit" is inserted after each string of message bits of a predetermined 
length k so that the sum of the (A^ + 1 ) bits is either always even (even parity) or 
always odd (odd parity). For example, if ^ 4 and even parity is used, the message 
110100101001,,, will become 110110010110010, , . On receiving the trans- 
mission signals, an error is detected if the sum of the appropriate five contiguous 
digits is odd. If the probability of error in a single bit is 1 percent, find (i) the 
probability of at least one error in the transmission of four sequential bits, and (ii) 
the probability of an undetected error after using even-parity coding. 

Solution: (i) The probability of an error occurring in at least one of the four bits is 

1 - P(no errors in the 4 bits) ^ 1 = (0,99)^ 

- 1 - 0,9606 - 0.04. or 4 percent. 

(ii) In each set of five bits under parity coding, if 1 . 3. or 5 errors occur, the sum of 
the binary digits will be odd and the error will be detected. If 2 or 4 digits are in 
error, this will go undetected, 

P(2 errors) ^ (0.01 )= (0.99)^' * 0.00097 
P(4 errors) ^ (0.01)^ (,0.99) ^ 5 x 10^^ 

so 

P(undetected error) ^ P(2 or 4 errors) — 0,1 percent 

By inserting a parity bit after each four message bits, we have reduced the trans- 
mission efficiency to 80 percent of its possible ma.ximum but have reduced the 
probability of an undetected transmission error in each four-bit '*word" from 4 
percent to 0.1 percent. However, when we do detect an error, parity coding does 
not tell us which of the bits is erroneous. In Chapter 8, '^Matrix Algebra," we 
shall examine the Hamming Code, which not only detects a transmission error but 
also tells which bit is wrong. 
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PROBLEM 4. An aerospace consulting company is working on the design of a spacecraft system 
composed of three main subsystems. A, B. and C, The reliability, or probability 
of success, of each subsystem after three periods of operation is displayed in the 
following table: 



A 
B 

C 



/ day 
0.9997 
LOOOO 
0,9961 



3.3 months 
0,8985 
0.9386 
0.9960 



8.5 months 
0.6910 
0,7265 
0,9959 



These reliabilities have been rounded to four significant digits. The 1,0000 in the 
first column means that the likelihood of the failure of subsystem B during the 
first day of operation is so remote that more than four significant digits are needed 
to indicate it, 

a. Consider the case of the series system shown in Fig, 5,1, If any one (or more) 
of the subsystems A, B, or C fails, the entire system will fail. If Is the total 
probability of success of the system, find for each of the three time periods. 



A 




B 




C 







Solution: For the first 24 hours. 



ERIC 



For a period of 33 months, 



For a period of 8*5 months. 



- (0.9997) (LOOOO) (0,9961) 
-0.9958. 

n - PaPbPc 

-(0.8985) (0.9386) (0.9960) 
^0.8400. 

H - PaPePc 

^ (0.6910) (0.7265) (0.9959) 

- 0.5000, 
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b. The system shown in Fig. 5,2 will succeed if B succeeds and at least one of A 
and C succeeds. Find the probability of success for this system for the 3.3-month 
time period. 







1 










1 


A 




B 




c 



A 




C 












Fig. sa 



Fig. S.3 
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Solution: The probability of success for the portion of the system containing A and C is 

Pax ^ 1 - P(both A and C fail) 

^ 1 - (0, 1015) (0.0040) - 0.9996. 

Then 

P, ^ PbPax ^ (0-9386) (0.9996) ^ 0,9382. 

c. For more complicated systems, the use of conditional probability is helpful. If 
an event A can be divided into n mutually exclusive subevents Bi, B2, , ^ Bp (n 
finite), then P(A) - F(A|B,) P(Bi) + P(A|B,) PiB.) H- . . . + PiA\B^) P(B,), where 
the notation F(X| Y) designates the conditional probability of X given that Y has 
occurred. 

Consider the system in Fig. 5.3, where the 3.3^month reliabilities of the sub^ 
systems A, B, C are the same as before and the 3,3-month reliabilities of D and 
E are 0,9216 and 0.9542, respectively. Use conditional probability to find the 
reliability (i .e., P,) of this system for the 3,3-month period. 

Solution: This system will succeed if any one of the paths (A,D), (B,D), (B,E), or 
(C,E) succeeds. We can choose B as our focus and assert that 
P, - P(system succeedsjB succeeds) Pg + P(system succeedsjB fails) (1 - Pb). Now 
we evaluate P(system succeedsjB succeeds). If the system succeeds given that B 
succeeds, this means that at least one of D and E would have succeeded, so 

P(system succeeds|B succeeds) - 1 - (1 - Pd) (1 - Pg) 

-1-(1 = Pd-=Pe + PdPh) 

= ft + Pe - ftPfi 

= (0,9216) + (0,9542) - (0.9216) (0.9542) 
^0,9964. 
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Next we evaluate P{system succeeds|B fails). For the system to succeed in view of 
the failure of B means that at least one of the paths (A,D) or (C,E) must have 
succgeded, so 

P(system succeedsjB fails) ^ I = (I - PaPd) (I - PcFe) 

- PaPd + Pcfls = PaPdPcPh 

- (0,8985) (0,9216) + (0,9960) (0,9542) 
- (0.8985) (0,9216) (0,9960) (0,9542) 

^0,9915 
P, ^ (0.9964)Pb + (0.9915) (1 - Pb) 

- (0.9964) (0,9386) + (0,9915) (0,0614) 
^0.996L 

PROBLEM 5. In Problem 4a we saw that the total reliability of the system deteriorates rather 
rapidly in its present stage of design, with less than a 50-percent chance that it 
will operate after 8,5 months. The reliability of subsystem C remains nearly con- 
stant, whereas the greatest decline in reliability takes place in subsystem A, 
which contains a particular part that is expected to wear out rapidly. The consult- 
ing firm, is asked to determine if enough improvement could be made in sub- 
system A to provide a reliability after 8,5 months of 0,7500. Compute the 
improvement needed in subsystem A. 

Solutions Let-* be the factor by which the reliability of subsystem A must be multiplied. Then, 
as before, 



0.7500 - (0.6910x) (0,7265) (0.9959) ^ 0,5000x 
0,7500 . ^f.^ 

^'^o3ooo"i-^°°' 

The reliability of subsystem A must be L500 x 0,6910 ^ 1.037. The increase in 
reliability cannot be obtained by improving subsystem A alone, since the 
reliability cannot be greater than 1, 

The next problem demonstfates the combined use of probability and computer 
simulation to determine the volume of an irregular solid. 
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PROBLEM 6. An internal fuel tank on a space vehicle has the shape of an ellipsoid tfuncated by 
three pla^ as shown in Fig. 5.4. Our problem is to detefmine the volume of 
this fuel Let us use, for an example, the ellipsoid whose equation is 

|2 + ^ + ~ - 1 and with the planes being .v - ±7, and z ^ =1,5, where the units 

are meter:^. 



Fig. 5.4 




a. If the tank is surrounded as tightly as possible by a rectangular prism with faces 
parallel to the planes formed by the coordinate axes, what inequalities must the 
coordinates of the ponts inside the prism satisfy? What is the volume of this prism? 

Solutions If Gr,>?,^) is inside the prism, x must satisfy -Kx <7 because of the truncating 
planes x ^ —1 and x ^l\y must satisfy =3 < < 3 because v = =3 and y = 3 are 
the planes tangent to the ellipsoid and parallel to the x-z plane; z must satisfy 
= L5 < z < 2, since z is bounded below by the truncating plane z - —1,5 and 
above by the plane z = 2 tangent to the ellipsoid and parallel to the x-y plane. This 
rectangular prism has dimension 14 m x 6 m x 3.5 m, and the resulting volume is 
294 m^ 

b. Let Vp be the volume of the prism and let be the volume of the tank, which 
we are seeking. If a point is randomly chosen inside the prism, express the proba- 
bility that it is also inside the tank, in terms of Vp and V^. 

Solution: This probability is equal to VJ V^, the ratio of the volume of the tank to that of the 
surrounding prism. 

If N points are chosen at random inside the prism and / of these points are also 
inside the tank, express V, in terms of M L and Kp, 

Solutions The probability that / points are in the tank out of the points chosen randomly 
inside the prism is approximated by UN. So we get UN - VJV giving 
K - V,(I/N). - 
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d. Write a computer program to perform this simulation, using a random number 
generator to get coordinates of points within the prism. 



Solutions 



10 REM MQLUME ilHULATION UiING 250 

PROBABILITY 

12 DIM K(15?I DIM VQLI25)i DIM P 250 

CT<25) ZBO 
iS M - I 

10 RtM INTERNAL FUEL TANK # TRUN 270 

CATEP ELLIPSOID 210 

□0 FIEM eOUNDARIES K^2/S*Z 4^ ¥-2 1000 

+ 1-^2/2*2 ^ I fNs -7#K=7 1010 



PRINT s PRINT 5 PRINT "TYPE 
ANY QTHIR LETTER IF FINIiHED 
INPUT ASjM^M .+ 1 



IF AS THEN K(H) 

1)1 GOTO iO . 

IF At "0« THEN 40 

END 

REM TEST POINT 
P ^ ' K # H / Y 



KCM - 



« Y / 8 + 







s;.- ■ . y - . - - ■ ^ 




2 # E / 4 -r ■ ■•: . ■ 




PRINT 


"THIS PROGRAM DDMPUTES 


1020 


IF P > m 1 THEN 1050 




THE 


MQLUME DF " 


1030 


I ^ . I 1 - ' . . 


34 


PRINT 


"A SPAGIERAFT"! INTERNA 


loao 


RETURN 




L FUEL TANK » 


1§Q0 


REM SUMMARY OF TRIALS^ 


36 


PRINT 


-WHDiE SHAPE IS A TRUNC -. 


1510 


PRINT s • PRINT " # • * * 1 




ATED 


ELLIPSOID- 




UMMARYT • # • »" 


37 


PRINT 


"MITH BOUNDARIEI 2^ 1 *g 


■1520 


PRINT l-^RINT ^'DD YOU WANT; 






-7»X§7i AND* 




TO SEE A SUMMARY OF THE";: V - 


31 


PRINT 


;'*K"2/S"2 PLUS + 2" 


1525 


. PR 1 NT " 1 1 MUL AT I ONI ; ID F ART « 




2/2' 


= 1" 




5 INPUT g« , ' - ' ' \ 


39 


PRINT 


— THE-USER WILL EHOOSE T 


. 1S30 


I F Bf Y ? THEN '1 BOO 



HE-NUMBER OF" r PRINT "POINTS 
■ ■ IN THE IIMULATION"! PRINT 
50 PRINT "HOW MANY POINTS DO YOU 

CHOOSE?" 
41 INPUT K(M> 
50 PRINT "COMPUTING. *." 
80 VP=la * 6 ♦ 3.S 
SO I - 0 ■ 
100 FDR N - 1 TO K (M) 
110 K s - 7 + la • RND M) 
120 Y^. - 3 :+ S ». RND ( I I . 
130 Z = - 1.5 + 3,5 * RND (IJ 
140 GOiUD iOOO : v 

150 NEHT N 

iBO ypL(M) = /yp * I / K(M)^ 

151 VVQL(M>iVOLtM) - 
1B5 PCTCM) = 100 • I / KIM) 
leS P% =- iCT(M) JPCT<M) P% 

170 PRINT II" OF THi "5K(M)5" PO 

INTS'! ' ^ ^ 
ISO PRINT rWERE IN THE TANK fFOR 

A 'VOLUME ESTIMATE"^ ' 
1S5 ' PRINT MOL(M)J** CUBic' UNITS, " 

ISO ^PRINT^ ■MHE .PIRGINTAGE OF PO 

I^4T1 IN THE TANK" 
200 PRINT "WAS "IPCT<M) 
202 GDSUB liOO. : ^ ' . . . 

205 PRINT i, PRINT » AND 

T HER TRIAL? 
110 PRINT i: PRINT i PRINT "IF YO 

;U WANT 'TO REDO THE S1MULATI0- - 

N» ^ 

220 PRINT "WITH THE S^ME^ NUMBER 

■ OF, POINTS iTYPE S". . 
230 PRINT : i PRINT WIT Hj A DIFFER 
ENT NUMBER •TYPED" 



iSaO PRINT ; PRINT s PRINT"- 

N PGT , VOL" 

15^5 PRINT . , ' 

15^0 FOR J ^ I TO M ' ' ^ ' 

1560 HTAg (5) J PRINT KCJJt 

1570 HTAB : ( IB) I PRINT PCTt J) 5 

lOao. HTAB '125)1 PRINT UOLiU) 

1530 NEKT J V - ; . 

IBOO RETURN 



]RUN ' ' ' „ 

THIS PROQRAM-CQMPUTES -TrtivVOLyHt > 0F^ 
A SPAOEGRAFT'S INTERNAL FUEL TANK^ ° 
WHPiE SHAPE IS; A TRUNCATEpUiLL I PiOIO 
WITH . BOUNDARIES] '^1 ,5 ,K*/-,7fXs7. ' ' j 
AND X?2yi*2 .Y-2/3*2;t'- 2?2>1*1'^1 ■ . 
THE UIER WILL CHOOSE, ^Hi 'NUMBlR iDf^- . 
POINTS IN THE SIMULflTJONt;r t'^ " 

* » * # #^suMMARY,?^^^L,#- 

DO "YOU WANT TO SeV 6^SUMHA8y*:b?vfH^^ 
SIMULATIONS ID, FAR?*, . i.'.V/--*i'Vv' "''i' i " " 
7 Y - 'V; ' - ^ " - ' .* '''J' 



N 

100 
100 

zoo 

200 
400 



per 

71 

Gl ' 

62 

64 

62 



208 
179 
'183 
189 
184. 
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PROBLEM 7. Sunspots were observed and recorded as long as two thousand years ago. The 
invention of the telescope around 1610 permitted the systematic observation of 
these solar features, their motion, and their frequency of occurrence, (Problem 6 
of Chapter 7 illustrates the use of trigonometry in analyzing sunspot motion,) It 
is relatively easy to observe sunspots by using a long-focus telescope to project an 
image of the Sun on a piece of white cardboard. 

Fig, 5,5 shows the record from 1610 to 1975 of what is now tommonly referred to 
as the "sunspot cycle. " The vertical scale represents the number of sunspots 
observed. The data since 1740 are considered reliable. 

Although sunspots are still not well understood, it has been established that they 
are regions in the solar atmosphere that contain enormous magnetic fields relative 
to their surroundings, along with cooler temperatures. Moreover, there appear 
to be connections between the level of sunspot activity and the occurrence of 
^'magnetic storms'' in Earth*s ionosphere, the density of Earth's upper atmo- 
sphere, and changes in Earth*s weather and climate. 

Since variations in upper atmosphere density can affect the orbital lifetimes of 
satellites, the prediction of sunspot activity is an important aspect of the plan- 
ning of some space missions. The mean cycle length as well as its variability must 
be taken into account, making statistical analysis vital to such predictions. 
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The following table sumniarizes sorne of the data of Fig, 5.5. The first step in the 
statistical analysis for the prediction of sunspot activity is to determine the mean 
and the standard deviation for each of the following measures: the rise time; the 
fall time; the period from minimum to minimum; the period from maximum to 
maximum. Compute these means and standard deviations. 



Table 5,1 





Year of 


Year of 




Year of 


Year of 


Cycle 


Minimum 


Maximum 


Cycle 


Minimum 


Maximum 


1 


1745 


1750.3 


12 


1867,2 


1870.6 


2 


1755.2 


1761,5 


13 


1878. 9 


1883,9 


3 


1766,5 


1769,7 


14 


1889,6 


1894.1 


4 


1775.5 


1778.4 


15 


1901.7 


1907.0 


5 


1784,7 


1788.1 


16 


1913.6 


1917.6 


6 


1798,3 


1805,2 


17 


1923.6 


1928.4 


7 


1810,6 


1816.4 


18 


1933,8 


1937.4 


S 


1823,3 


1829,9 


19 


1944,1 


1947.7 


9 


1833,9 


1837,2 


20 


1954.2 


1958.2 


10 


1843,5 


1848,1 


21 


1964.6 


1970.6 


11 


1856.0 


1860.1 









Solution' The computations were done by mierocomputar. The program Hsting and the results 
of the run are shown below and on page 82. 
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130 HTAB Sfli PRINT AAV? 
i3B FD^ J s 1 TO M 
laO ZiJ) ^ 5(J) 

NiKT J 
ilO SQ1U5 □000 
f iO BAU 3 ZAU : 
570 HTA§ 3as FRINT BAM 
iOO REh CQMFUTl AND PRINT STAND 

ARD DfMlATIBNi 
605 DIM D(H> 
SIO PRINT "liD* •* I 
620 FOR J = 1 TO H 
^□0 ZCJ) * R(J5 

iflO NlXT J : 

BSO ZAM ■ SAU 

BB0 OOlUB ^000 

170 RID « 2SD^ 

680 HTA5 Si PRINT RtD I 

eeo FDR J - 1 TO rtiZ(JJ « FfJ5i NIKT 
J 

700 lAU =1 FAV 
710 GOSUi 4000 

720 FSD ^ZBSs.HTAl iBi PRINT FID 



1 ) 



TO nil CJ3 SAC J) 4 NiXT 



730 FOR 

750 ZAU^ AAU 
750 GOBUB 4000 

7B0 ASP = 2SDr HTAB 24 J PRINT Ago 

3 : ^ . .. . . - . 

770 FDR i TO; Ni2(J) > B( J>i NEXT 

J 

7B0 lAV =■ BAy 

760 GOSUB 4000 . 

sop BSD sZiOr HTAB 34t PRINT BSD 
610 END 

REM- ViCHO INPUT. DATA : 
HOME 

PRINT ^CYCLE VR OF MIN 

.YR OF MAX" 
FOR g = 1 TO M - 
HTAB 3i PRINT J I 
HTAB Iti^PRINT X( J) 1 
HTAB 25s, PRINT Y(J>i = 
NEXT" J 
RETURN 

REM PRINT YEARLY RllE ?^^TlHr 
.FALLtlME iHlN-HlN PERIOD IMAX 
-MAX PERIOD 

PRINT 
RS 
PRINT 



1000 : 

1010 

1020 

1030 
1040 
1050 
1060- 
1070 
1060 
2000 



2010 



2020 



3050 Z = SUM / (M - 
3060 Z% ^ 100 * 2 
3070 1AM = Z% / iOO 
3060 RETURN . » 

4000; REM COMPUTER , STANDARD Df Ul" 

ATIONS TO 2 DECIMAL PLACES 
4020 Dili » lC4K-dAV^ - ^ 

4030 SUM ■ o<i> **'0'ii'r 'V-f/ " ; 

4040 FOR J ^ 2 tJ&'-'lH 1 } . 
405O DiJ) « Z(Ji i..,tAM' — -w',/- 
4010 SUM » SUM + l£0iJi.r4'DCJU_ 
4070 NEXT ' ^ 
4060 SO 

4030 SDZ « 100 ID I 

4100 iiD s tor / -loo 

4110 RETUl^N ^-r- 



J RUN 

CYCLi 
1 ' 
2 
3 
4 
5 
B 
7 
B 

a 

10 

11 

" 12 
.13 . 

14 

15 

1^; 

17 

18 
-13 

20 

21 

ALL TIMi 




1823.3 

1833|1 

1143.5- 

1858 

18G7.2 

1678*3 

1883.6- 

1313. 

132Q,6 

1933.8 

1344; 1' , 



1317. B f , 



MESl ARE^aN?^YlA^i^^^^r'V%^' , 



2030 



2040 

2050;^ 
2080 




FOR 
HTAB: 

HTAB Jj^FRINTt'RtJ) I 
2070„ HTAi,"lBs .PRINtJfijJ T 

20B0 /^HT^Pi^atip^RiNT^Acjri : 

2010 MNTAB*C.34i;4>RINT-Bf J) T ' 
21 00--' NEXTV4 ■ '-y-^'! < 'I A* . ^ 
1110 ^.RETURN.- 
□000 Rr ^ — * ^ 

3010 SUM 
3020 FOR 

3030 SUM 

.3040 NEXT 
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PROBLEM 8. Among the studies arising from Landsat observations are several concerning 
the evaluation of properties of snowpacks. In many areas of the world, water 
resources are heavily dependent on winter accumulations of snow. 

Computer models are being developed whereby potential water resources can be 
predicted from satellite measurements of microwave emission in snow-covered 
areas. Predictions from such models are tested and the models refined by making 
comparisons with'ground-based measurements of snow depth and temperature. 
Such measurements, when graphed, inevitably show a large amount of scatter, 
and it is the regression line for the datp. that is used as the standard for comparison. 

Fig. 5,6 shows such a comparison, where the horizontal scale is temperature in 
degrees Kelvin. (The Kelvin scale of temperature is obtained from the Celsius 
scale by adding a constant, 273.15, so that O^C ^ 273.15^K, and 
10^C-2S3.15*K0 

The data points of Fig. 5.6 are listed below. Find the parameters of the equation 
of the regression line. 



30 ^ 



25 



S- 15 




Model 

Gencrnted 

Curve 



Model generated curve ef the snow depth versus 
ihe actUiil dcpihs inken with trusk mouniod or sir^ 
er^ft senHors- 



210 



220 



230 
Tb(K) 



240 



250 



260 



Fig. 
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Chapter Five 



Solution: 



As is the previous problem, the computations were done by micrucomputer. T* 
program listing and the results follow. 



3LXiT 

10 REM SNOyPACK HODEL 
15 PRINT 

10 REH RlGREeSlON LINE 
,30 DATA 195»15»207iig*10S,l5il0 

S>17ill0>2S 
aO DATA 111 >20*111 .11 ,lia,21 ill 

i iii iiii »ia 

so DATA 118 » 16 #113 .la .117 ,20 *li 
. U 8j22#130»13 

SO DATA 231»13»231 ill »232 # 17 #23 
; - 3iiB*233>llfl33fBfl35f21 fl3S 
1^ *?'^^^»^»237fl2' 
. 70 DATA . 23Bf8 »l3i ,li ,l3g ,s ,250 * 

: e i240>a »2ap ts #2^1 izaz *5 #25 

3 V7 >ia3 ie»2a3#3 
80 DATA ia3»l»255^1 ,1^7, 2.157 »3 
: tZHB »a »15S >5 #250 il i2Sl .3 i2I2 
■ : ■ iZ »255»1 

eo DIM K(ae>£ dim vcai) 
. 100 - .PDR i = i TO 

110/ READ KC I) 2 MEAD ¥(1) 
120% NEKT I . 
130 GOSUB 1000 

140:KSUM> OiVSUM = OiSP ^ 0 i i© = 

SiO FOM i 3 i TO as 

leo HSUH ^ XeUH + KU )t¥SUM = YiU 

. m:+ Y(ii 
170 SP = S> + K(I) # VfDsSO g go 

+ HCn * ) 
liO NEXT J 

ISO XMIAN ^- XiUH / 4S;Yf1EAN ^ YiU 

/ M / as 

200 B = (IP al « KMEAN • VHEAN) 
/ (SO - fli # HMIAN » KMIAN) 

V 2i0 B% =. 1000 * BsB^ B% / 1000 
i2X> A% = 100 * YMEANiA = A% / 100 

;225 Ci;^ 10© • KMEANiC = C% / iOO 

130V PftlNT " THI REORESSION LINE 
, ■ ; HAS ILDPE " IB 
240 PMINT " AND A MEAN iNOW DBPT 

; , H 0P>« iA 
2Sq' PRINT "EORREiPONDe TD A HEA 
; N TEMPiRATURE"s PRINT " DF " 

„. . • '"iG-i. - : , ^. ■ 

leo END 



1000 RtH ECHO INPUT DATA 
1010;: PRINT " J«i"HCn" "Y(IP 
1020 FOR I = 1 TO as 
1030 . PRINT I ,HM J iVC n 
iOaO NEXT 1 

logo PRINT '*1ND OF DATA" 
lOiO RETURN 



z 


X ( I ) 


Y C I ) 


1 


105 


25 


2 


207 


19 


3 


206 


15 


a 


206 


17 


B 


210 


25 


e 


211 


20 


7 


211 


21 


i 


2ia 


21 


i 


217 


22 


10 


216 


16 


1 1 


216 


16 


12 


223 


14 


13 


217 


10 


la 


221 


22 


15 


ISO 


13 


16 


132 


13 


17 


231 


16 


16 


131 


17 


19 


133 


It 


20 


233 


11 


21 


233 


6 


22 


235 


21 


13 


236 


6 


2a 


237 


5 


25 


137 


12 


26 


13i 


s 


17 


13S 


IB 


it 


233 


6 


IS 


2ao 


i 


30 


2ao 


4 


31 


2ao 


8 


32 


2ai 


a 


33 


laz 


s 


3a 


las 


7 


35 


ia3 


6 


36 


ia3 


3 


37 


2a3 


1 


3B 


ia5 


1 


33 


247 


Z 


ao 


2a7 


3 


ai 


248 


4 


a2 


las 


4 


a3 


150 


Z 


aa 


151 


3 


a5 


252 


2 


46 


155 


1 


END OF DATA 







THE RiGRESilDN LINE HA6 iLQPE - , a^ 
AND A MIAN BNQM DEPTH OF 11.3 
GORRESPDNOe TO A MEAN TEMPERATURl 
OF 232. S3 
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APOLLO 1 7 EVA— Astronaut Eugene A. 
Cernan, Commander of the mission is photQ= 
graphed by Astronaut Schmitt whose photo is 
reflected in the gold visor. 
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Chapter Six 



The early work that led to our undefstanding of the planetary motions and gave 
us the description of the solar system we know today would have been virtu- 
ally impossible %vithout the use of logarithms to reduce the labor of the com- 
putations. Although computers and calculators have replaced logarithms as com- 
putational tools, logarithmic and exponential functions are still essential for the 
study of Earth's atmosphere and rocket propulsion, examples of which are cited 
in this chapter, 

PROBLEM 1. Experimentation and theory have shown that an approximate rule for atmospheric 
pressure at altitudes less than SO km is the following: Standard atmospheric pres- 
sure, 1035 grams per square centimeter, is halved for each 5.8 km of vertical 
ascent. 

a. Write a simple exponential equation to express this rule. 

Solution: Letting P denote atmospheric pressure at altitudes less than 80 km and h the altitude 
in km, we have 

P - 1035(i/2)^^^-^g/cm^ 

b. Compute the atmospheric pressure at an altitude of 40 km. 

Solution: From the equation of part (a), 

P ^ 1035 (l/2)-*"^-' ^g/cm= 
- 1035 (l/2)^^g/cm2 
^ 1035 (0.0084) g/cm= 

^ 8,7g/cmv 

c. Find the altitude at which the pressure is 20 percent of standard atmospheric 
pressure. 

Solution: Substituting in the equation of part (a) gives (0,20) (1035) ^ (1035) (1/2)'^^' ^ where 
h is in km, and so (0,2) ^ (1/2)^^^ ^ Now, taking logarithms. 



and 



log (0,2) log (0,5) 



log (0.2) 

= 5-8 , )n Js km - 5,8 (2,32) km - 13.5 km. 
log (0,5) 
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PROBLEM 2. The rule for the variation of atmospheric pressure with height which was given in the 
previous problem can also be written 

P - 1035 (lyf^!^^ 

Atmospheric scientists often use this rule in one of its equivalent forms where the 
base is 10 or e, the base of the natural logarithms, instead of 2. Find and so 
that P ^ 1035 (2)^" ^ 1035 (10)^*'^ ^ 1035 (e)^^^^ 

Solution: We need to find ki so that 2"^^ ^ 10^*. Taking logarithms, 0.17 log 2 ^ ky or 
&i = (0.17) (0.301) - 0,051. For it. we have 2« >^ - or = (0.17)log,2 
= (0.17) (0.693) ^ 0.12, 



Exponential and Logarithmfc Functions 



PROBLEM 3. Sometimes different bases are used together in the same application in atmospheric 
work. For example, atmospheric absorption of electromagnetic radiation from 
the Sun and other sources is dependent on the wavelengths of the incoming radi- 
ation. Instruments carried by rockets, balloons, and satellites have shown how 
far in the atmosphere such radiation penetrates before being reduced by a factor of 
1/e, the conventional measure used in this work. The results are given in Fig. 6,L 
Both the wavelength scale and the altitude scale are logarithmic, with the 
horizontal scale in base 10 and the vertical scale in base 2. (How much of this 
information could be displayed using linear scales even on a wall-sized chart?) 

Fig, 6,1 shows that visible light and radio waves penetrate the atmosphere com- 
pletely and reach Earth's surface. However, gases such as oxygen, ozone, nitro- 
gen, and water vapor absorb most of the infrared, ultraviolet. X-ray, and shorter 
wavelengths. At what altitude will solar infrared radiation of wavelength 10^^ m be 
reduced by a factor of 1 /e? 



10^ 



10-2 



1 0-^ Ur^" !U ^= Ur'^ Wavy lengm 
^ ' ' ' ' ' in Meiers 




hurth prontc 



I level 



Fig. 6.1 
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Solution: The equal intervals on the altitude scale have length log 2, The ordinate we are 
seeking, J, is 1/4 of the way between log 50 and log 100. This means that 

logj' ^ log 50 + |log2 ^ L699 + i(0,3010) 

= L699 + 0.U75 ^ 1,774. 
Theny ^ 10'^^^ ^ 59. 

If a calculator with ^y ' key Is available, we can solve this problem without 
actually finding logarithms, as follows: 

log>' - log 50 + I log 2 - log [(S0)(2)' - log [(50) (L 189)] 

= log 59, Soy ^ 59, 

In the foregoing problem, we saw how the use of logarithmic scales made it possi- 
ble to display information over an extremely large range of values. The next two 
problems show another use for logarithmic scales, that of fitting a mathematical 
function to experimental data. 
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PROBLEM 4* Very high energy p.^rticlcs (electrons and protons) are found in the radiation belts of 
some planets (e.g, , Earth, Jupiter. Saturn), and a plot of the number of particles 
found at different energies is called a spectrum. Often the spectrum has a shape 
that can be represented by an equation of the form N = KE'^ where N is the 
number of particles at a certain energy, E; ^ is a proportionality factor; and m is 
called the spectral index. 

When the spectrum has such a shape^ we call it a power-law spectrum, and the 
experimenter studying such a spectrum wants to know the values of m and A'. 
Table 6. 1 shows values of measured several £s during the flight of Pioneer JO 
past Jupiter, For these data, find the besi value of m and of (N is really the 
number of particles hitting a detector per unit time, or the counting rate, which is 
why the number can be u fraction.) 

Table 6, 1 



Energy, E 


Number^ 


, jV 




LU X 


10" 






H/ 


UM) 




10^ 


KO 




10^ 


L6 


1.0 X 


io^ 


4.5 20 


10,0 1 


20,0 


OA 





Solution: Using logarithms on the expression ^ A'^""' results in log A/ — log K + mlogE, or, 
to obtain the form of a linear equation y ~ nLv ^ b, log N ~ m log E + log K. 

We can find logarithms for the values of E and in the table (or we can use 
log-log graph paper and circumvent this step), plot the points, and draw the best 
straight line through this set of points. Then m will be the slope of the line, and 
K will be the value of N for which log £ - 0, (Note that this is the value that lies 
on the best straight line, and not necessarily any value in the data set.) 

We observe in Fig, 6.2 that the intercept on the log N scale is 3,5. Since log 

N ^ 3,5 when log E ^ 0, we have log K = 3,5, m K = 10^^ ' ^ 3200. The points 

(0,3,5) and (1.0,0) are on the best fit line, so m ^ ^zTi q ^ "3.5. So 
- (3200) £^'-' 
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Log N Log E 




Fig. 6.2 
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Solution: 
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Many of the control functions in a space vehicle system are automatic, handled 
by computers and electronic feedback devices. However, the adaptability and the 
decision-making ability of human monitors of these systems, whether crew mem- 
bers or ground controllers, greatly increase the overall effectiveness of mission 
control. Because of this, managers of projects in the space program have become 
interested in some results from psychological studies of human decision making. 

One such study measured the time it took to respond when faced with varying 
numbers of choices. Experimental results are given in the table below, where N 
is the number of choices presented and R is the reaction time in seconds. Graph 
these data on semilogarlthmic graph paper with N on the logarithmic scale (or 
graph R against logioA^ if semilog graph paper is not available) and find an empir- 
ical expression for reaction time as a function of the number of choices. 

^^123456789 10 
R 0.17 0,34 0,37 0.42 0.48 0.52 0,56 0.58 0.59 0.57 

The points are graphed and a ^^best fit" line drawn (see Fig. 6.3). Since the point 
(N,R) - (1,0.17) does lie on this line, we have R ^ 0.17 + m logi„M To find m, 
we can use the points (1,0,17) and (9,0.59), since both are on the "best fit'' line; 

0.59 ^ 0.17 + m login9 

so 0.42 ^ m (0.954), 

or m ^ (0.42)/(0.954) ^ 0.44 

The raquested relation is ^ ^ 0.17 + 0.44 loginM 

As we have seen in Chapter 4, solar cells, which convert solar energy Into elec- 
trical energy, can be used to supply power in space vehicles. Nuclear energy 
derived from radioactive isotopes is also used. Nuclear energy sources gradually 
lose power in a manner described by the exponential function. The next problem 
illustrates some computations of the available power and operational life of a 
satellite using a nuclear power source. ^ ^ 
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PROBLEM 6, A sateriite has a radioisotope power supply. The power output In watts is given by 
the equation 

where t is the time in days and e is the base of natural logarithms, 
a. How much power will be available at the end of one year? 
Soiution: Applying the given equation, we have 

P ^ 50e^^^^^'^" 

^ 50 X 0,232 
= 11,6 

Thus approximately 11,6 watts will be available at the end of one year, 

b* What is the half-life of the power supply? In other words, how long will it take 
for the power to drop to half its original strength? 

Solution: To find the half-life, we solve the equation 

25-50e^^^2sn 

for r and obtain 

- --0.693 

f ^250x0.693 

- 173. 

Thus the half-life of the power supply is approximately 173 days, (Note that In x 
is a shorter expression for log^ x,) 
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c. The equipment aboard the satellite requires 10 watts of power to operate 
properly. What is the operational life of the satellite? 

Solution! Solving the equation 
for f gives 

-In 0.2 
- -1.609 
f ^ 250 X IM9 
^402. 

Hence the operational life of the satellite is 402 days. 

PROBLEM 7. The theory of rocket flight shows that the velocity gained by a launch vehicle when 
its propellant is burned to depletion is expressed by the equation 

V ^ c In ^ 

where v is the velocity gained by the vehicle during launch; 
c is the exhaust velocity of the engine; 
In R Is logg/?. or the natural logarithm of Ri 

J n *i ■ f i i t_ taReoff Weight 

and R is the mass ratio of the spacecraft, defined by ^ ^ ~— . 

burnout weight 

a. The takeoff weight consists of propellant or fuel. F, structure, S, and payload, 

R At burnout, assuming all the fuel has been used, the remaining weight is S + 
^ + 5 + p 

so that R = p ' In general, the weight of fuel cannot be more than about 

10 times the weight of the structure in order for the vehicle to withstand the 
stresses of operation. Show that if F = lOS, then an upper limit for is IL 

Solutions If F ^ 1 05, then R = — ~ — ^ - 

5 + F 5 + F 

11(5 +F) " lOF 
5 +F 

So the largest possible value for R, is 11, but we see that in order to actually 
achieve this value, it is necessai^''^or F to be 0— in other words, the launch vehicle 
could carry no payload! 
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b. The minimum altitude for a stable orbit about Earth is about 160 km. At lower 
altitudes, air resistance slows the spacecraft and causes a rapid detarioration of 
the orbit. As will be shown in Problem 1 of Chapter 9, the spacecraft must attain a 
velocity of about 7.8 km per second to orbit at 160 km. However* in order to 
overcome the retarding effect of Earth's atmosphere while the spacecraft is 
ascending, the total velocity imparted by the launch vehicle must be at least 9.0 
km/s. What is the minimum exhaust velocity needed by the rocket engine if 
^^11? 

Solutions Substituting v=9.0km_/s 

and ^ 11 in the rocket equation^ 

9.0km/s^cln 11 

9.0 , , 9.0 ^ ^ , . 
c = : — ^ km/s ^ — = 3.8 km/s. 
In 11 2.4 

c. The propellants used for engines such as those of the Delta, Centaur, and 
Saturn launch vehicles could produce exhaust velocities averaging at most 3 
km/ s, which would not be sufficient to achieve orbit. The main engines of the 
Space Shuttle use a mixture of liquid hydrogen and liquid oxygen, which can 
produce exhaust velocities of 4.6 km/s. However, in order for the Shuttle to per- 
form its tasks and return to Earth with its crew, it has an /^-value of around 3,5. 
Could the Space Shuttle achieve orbit with its main engines? 

Solutions If c ^ 4.6 km/s, and 

R -3,5, then 

V ^ 4.6 In 3,5 km/s 

- (4,6) (L25) km/s 

^5.8 km/s, 

which is not sufficient for orbit. 
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PROBLEM 8. It is apparent from the rocket equation that the burnout velocity increases when the 
mass ratio increases. We can get a higher mass ratio by using a solid propellant 
because the stiff, rubberlike pfopellant mass serves as part of the structure. If no 
payload, or a very small payload, is included, a solid-propellant rocket could 
have a mass ratio of about 19. A typical average exhaust velocity for a solid pro- 
pellant might be about 2.4 km per second. Could this launch vehicle achieve a 
160 km Earth orbit? 

Solutions Using the rocket equation, 

V ^2.5 In 19km_/s 
-(2,4) (2.94) 
^7,1 km/s, 

which is much less than that needed for orbit. 

The solution to the problem pointed out in the preceding examples is to use stag- 
ing. That is, the launch vehicle is divided into two or more parts, or stages. As 
soon as the propellant has been burned in the first stage, there is a brief coast 
during which the heavy motors and structure in the first stage are jettisoned and 
permiited to fall into the ocean. Freed from this deadw^eight, the second-stage 
motors are much more effective; the same procedure is repeated for the remain- 
ing stages, 

PROBLEM 9. Let us.design o two-stage vehicle to place a payload into Earth orbit. We shall make 
some simplifying assumptions to make this problem easier while preserving the 
basic idea; (1) the structure weight of each stage is 10 percent of the fuel weight, 
the remaining weight being payload; (2) the gain in velocity is divided equally 
among the stages, each contributing 4,5 km/s to the required final velocity of 9.0 
km/s; (3) all stages use the same propellant with an exhaust velocity of 3.7 km/s. 
This third assumption is generally not true in practice^ — ^for example, the Space 
Shuttle uses solid rocket boosters in addition to the main engines— but our goal 
here is to see how staging works. For the sake of having a numerical example, we 
shall also assume that the total weight at liftoff is 5.0 x lO'* kg. For this numerical 
example, determine the weight of fuel to be carried by each stage, the structural 
weight of each stage, and the weight of the orbital payload. 
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Soiutioni Let Fu Su Pi represent fuel, structure, and payload weight, respectively, of the first 
stage, and F;, S3, and P2 those of the second stage. Since the '^payload" of the 
first stage includes the entire second stage and the orbital payload, 

Pi^ F. ^ S2 + Pi^ 

First stagei v ^cln 

4.5 ^3.7 In ^1 

In R,^~^l.22 

^1 -e'-=-3.4 
So F,+S, + F, _ 5.0X10^ 

Then Si + F| - - L5 x 10^ kg 

and Fi ^ (5,0 = IJ) x 10^ ^ 3.5 )c lO"* kg. 

By assumption 1. Si ^ 0.10 (3.5 x 10^) -3. 5 x 10^ kg. 
Then Pi - IJ x 10' = 3,5 x 10-^- 1,15 x 10' kg. 

Second stagei We again have, from the rocket equation, 

4.5 ^3.7 In ^2, 
so ^2^3,4, 

L15 X 10-* 



Also, ^2 



S. + P. 



Then S. + 1-15 x 10^* ^ 3^ ^ 

3.4 * 

Therefore, F.^ l,15x 10-*-3.4x 10'^ 8,1 x 10^ kg 

Sa^ 0.10 (8.1 X 103)^0.Sx 10-^ kg 

P. - (3.4 = 0,8) X 103 ^ 2,6 X 10^ kg. 
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Our design for the two-stage launch vehicle may be checked as follows: 



Weight of fuel: 
/^ i ........ . 

Total .... 



kg K 10^ 
35.0 
8T 
43T 



PROBLEM 10, 



Solutions 



PROBLEM IL 



Weight of structure: 

5i , 3,5 

Total 43 
Weight of orbital payload 2,6 
Total weight of vehicle 



50J^5.0x 10^ kg 



Thus, although the single-stage launch vehicle discussed in Pfoblem 7 could not 
place any payload into orbit, this two-stage vehicle can place nearly 5 percent of 
its weight into Earth orbit. 

Show that when all stages use the same propellant, the total mass ratio of a multiple- 
stage launch vehicle is equal to the product of he individual mass ratios. 

Indicate the burnout velocities and mass ratios of the first, second, third stages, and 
so on, by the subscripts 1, 2, 3, and so on. Then, using a three-stage vehicle as an 
example, 

Vi + v'2 + V3 ^ c In /?i + cln /?2 + c In R3 
V ^c(ln^i + ln^2 + ln/?3) 
y ^clog,(^i^2^3) 

(Note: Making the structure stronger so that it can support large payloads feduces 
the mass ratios. However, if we have several stages, the total mass ratio can 
become very high, producing much greater performance,) 

Using the equation derived in Problem 9, show that the launch vehicle constructed 
in Problem 8 can indeed orbit its payload. 



Solution! Given 



^i^3-(3,4) (3.4)-ll,56 
V -2,7 1ogalL56 
^3.7(2.45) 
^9.06km/s 

Th^ launch vehicle will impart sufficient velocity to overcome drag losses and 
insert the payload into a 160-km Earth orbit. Note that dividing the launch 
vehicle into stages increases the overall mass ratio to 11,56. 
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ngle measurements and the trigonome tric analysis of such measurements are 
m[M used extensively in space science. Among the examples we shall consider 
B she re are some involving tfansformations betsveen terrestrial (or celestial) and 
spacecraft coordinate systems, a variety of photogrammetric corrections, and the 
tracking of spacecraft from stations on Earth, 

PROBLEM 1* A conventional right-handed three-dimensional spacecraft coordinate system is 
shown in Fig, 7. 1 . The angular motions of the spacecraft with respect to the .v-, 
y-, and z-axes respectively are called roll, pitch, and yaw, shown in Fig, 7,1 by 
curved arrows. We shall develop the transformations between this coordinate 
system in a moving spacecraft and a reference coordinate system whose origin 
coincides with the one in the diagram but does not undergo rotation. Here, we 
shall consider a single rotation at a time. In Chapter 8, "'Matrix Algebra/' we shall 
investigate a series of such rotations. 

When the spacecraft performs a rotation, the reference system remains fixed, but 
the spacecraft coordinate system undergoes the same rotation as the spacecraft. 
If the point Q has coordinates (x, y. z) in the reference system, we need to find its 
coordinates in the spacecraft system after such a rotation takes place. Let us 
consider each of the motions roll, pitch, and yaw separately = 

a. Let the spacecraft coordinate system initially coincide with the reference sys- 
tem, and let the spacecraft undergo roll through angle R. Express the coordi- 
nates (x^, zr) of a point Q on the spacecraft in terms of (x, y, z) and R after this 
motion is performed. 
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Solutions Since the roll is around the jf-axis, the jc -coordinate of Q is the same in both systems: 
^ X, Now consider the plane parallel to the y-z plane, which contains Q, The 
roll moves Q to Q' as shown in Fig, 7»2, Let r ^ OQ - OQ' and let ^YOQ' ^ 9, 
Then Z^Y^^OQ' ^ B ^ R, Q' has coordinates (>% z) in the reference system, where 

y = r cos B and z — r sin &, 
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In the spacecraft system, Q' has Co ^rdinates ( Vr. zr) 

^vhere vr ^ rcos {G = R) and zr ^ r sin (0 - R). 

Expanding the sine and cosine of this difference resuits in 

Vr ^ r cos 6 cos R + r sin 0 sin ^ ^ y cos R + z sin R 
Sr = r sin 0 cos R = r cos 0 sin — z cos R — y sin R. 

b. Find the comparable transformations if the rotation is either a pitch through 
an angle P or a yaw through an angle K 

Solutions For a pitch rotation, this takes place around the y-axis, so if the coordinates in the 
spacecraft system are (.Vp, Jp, Zp), we have yp = y. We next consider a plane 
parallel to the x-z plane, and the analysis will be just as in part (a) with y replaced 
by z, z replaced by x, and Z.R replaced by resulting in Ap - a- cos P - z sin 

P; Zp ^ z cos P + X sin P. 

A yaw rotation takes place around the z-axis, so if the coordinates in the space= 
craft system are (xy, yy, Zy), we have zy - z; now we consider a plane parallel to 
the x-y plane, and this time the analysis is just as in part (a), with y r eplaced by x, 
z replaced by y, and replaced by The result is Xy ^ x cos Y + sin Y; 
yy - y cos Y - x sin Y. (We note that the right-handed system dictates that a 
positive angle of rotation take place so that the cyulic order x y z x is maintained.) 

c. An Earth-based computer monitoring the coordinates of Jupiter in Voyager^s 
reference frame recorded Jupiter at (2.03, -2.81, 0,336) (in units equivalent 

to 10^ km) at one point. If Voyager had performed a yaw rotation of 28^ just prior 
to this reading, what were Jupiter's coordinates in the spacecraft 
coordinate system? 

Solution! Using (a% z') for the spacecraft coordinate system: 

x' - A cos 28^ +>^sin28^ 

- 2.03 cos 28^ - 2.81 sin 28^ 

^ 0,473 

y' ^ - A sin 28^ ^ y cos 28^ 

= =2.03 sin 28^ - 2.81 cos 28* 

^ -3.43 

^ z ^ 0.336 

So the coordinates were (0.473, -3,43, 0.336) in these units. 

We next calculate the length of some of the latitude circles on Earth, 
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PROBLEM 2. Although Earth is not really a sphere, it can be treated as though it were spherical 
for many purposes. 

a. Show that the length of any parallel of latitude around Earth is equal to the 
equatorial distance around Earth times the cosine of the latitude angle (see 
Fig. 7,3), if we assume a spherical shape for Earth, 




Fig. 7,J 



102 



Solutions By the definition of the cosine function, cos 6 — r/R, or r = R cos 9. The length 
of the parallel of latitude is Cp, If Q denotes the equatorial circumference of 
Earth, then 

Cp - 2Tfr 

^ litR cos 6 
= Cg cos 0. 

b. Find the length of the 30^ parallel, north or south latitude. Use R = 6400 km. 
Solution: Applying the formula for the length of a parallul of latitude derived in part (a) gives 

Cp - (6400 km) (cos 30^) 
^ (6400 km) (0.866) 
^ 5500 km. 

Determine the length of the Arctic Circle (66^33' N). 
Solutioni Using the formula from part (a), the length is 

Cp = (6400 km) (cos 66^330 
^ (6400 km) (0,398) 
2500 km. 
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d. How far is it "around the world" along the parallel of 80- north latitude? 
Solutions Using the result of part (a), the distance is 

Cp - (6400 km) (cos 80^) 
^ (6400 km) (0,1737) 
^ 1100 km. 

PROBLEM 3* Two tracking stations s miles apart measure the elevation angle of a weather balloon 
to be Of and respectively (Fig, 7,4), Derive a formula for the altitude h of the 
balloon in terms of the angles a and jS, Ignore Earth's curvature. 




Fig, 7,4 



Solution! Writing an equation for the cotangent of each angle and solving for x gives 

cot a ^ - --- - - 
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X = h cot a — 5 



and 



X - ^ cot p. 

Now the two expressions for x are equated: 
/i cot a - s = A cot jS 



so 



and 



h (cot a - cox ^ s 



h ^ 
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cot a — cot p 
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PROBLEM 4* A satellite traveling in a circular orbit 1600 km above Earth is due to pass directly 
over a tracking station at noon. Assume that the satellite takes two hours to make 
an orbit and that the radius of Earth is 6400 km. 



a. If the tracking antenna is aimed 30- above the horizon, at what time will the 
satellite pass through the beam of the antenna? (See Fig, 7.5.) 




Solutions In the triangle formed by the station, the satellite, and the center of Earth, y ^ 120^ 
From the law of sines, 

sin a _ sin y 
6400 ^ 8000 

6400 sin 120 ^ ^ 

8000 "Q-^^^- 

Then 

a ^ 44^ 

and 

^ ^ 180^ - (120^ + 44^) ^ 16^ 

The time between 

p - 16^ and p - 0.0^ is ^ (120 min) 
— 5,3 min. 

This means that the satellite will pass through the beam of the antenna at 
12:00 - 5,3 minute^ o- 11:54.7 a\m. 
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b. Find the distance between the satellite and the tracking station at 12:03 p.m. 
Solutic m Computing angle p gives 

- 120 min 

By the law of cosines, 

x" - (6400)- + (SOOO)' - 2(6400) (8000) cos 9^ 

^ (40.96 + 64 " 10Li4) x 10" km^ 

^ 3,82 K 10^km= 

X ^ L96 K 10^' ^ 2.0 K 10^' km. 

We have found that the distance betw^een the satellite and the tracking station is 
2000 km (to two significant figures) at 12:03 p.m, 

c. At what angle above the horizon should the antenna be pointed so that its 
beam will intercept the sateUite at 12:03 p.m.? (See Fig. 7,6.) 



12.03 p.m.— ^ Noon 




Solutions Again, applying the law of sines, 



sin 9° ^ sin (y + 90°) 
2000 " 8000 



sin (r + 90") = sin 9= - 0.626 



cos y - 0,626 

y = 51°. 
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PROBLEM S. Two of NASA's tracking stations are located near the equator; one is in Ethiopia, at 
40- east longitude, another near Quito. Ecuador, at 78' west longitude. Assume 
both stations, fepresented by E and Q in Fig. 7,7, are on the equator and that the 
radius of Earth is 6380 km, A satellite in orbit over the equator is observed at the 
same instant from both tracking stations. The angles of elevation above the hori- 
zon are 5^ from Quito and 10^ from Ethiopia, Find the distance of the satellite 
from Earth at the instant of observation. 




Solutions In Fig. 7.7, OQ = OP^ OE = 6380 km: ^QOE = the longitude difference of the 
two stations, so ^QOE ^ 78^ - ^ 118^ Since 4QEO is isosceles, 

^OQE ^ ^OEQ ^ i (180^ - 118^) ^ 3r. 

Further, since the horizon is perpendicular to the radius, ^SQE ^ 5^ + (com- 
plement of ^EQO) ^ 5^ + 59^ ^ 64^ and Z.SEQ ^ 10' 4^ (complement of 
^QEO) - 10^ + - 69'. Also, ZQSE - 180' = (64' + 69") - 47'. These angles 
are all shown in Fig. 7.7. We are looking for the distance SP. If we can determine 
OS, then SP ^ OS ^ OP ^ OS - 6380 km. We note that OS is not an angle bisec- 
tor for either Z.QOE or ^QSE, so we must use an indirect method to find OS, 



We can evaluate QE from .^^^^ ^ 



then 
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sin 118" sin 31"' 

6380 sin 118' 



QE = 



sin3r 



^ 1,094 X 10^ km; 



SE _ QE 



sin 64^ sin 47" ' 

1.094 X 10-* sin 64' 



SE ^ 



sin 47" 



1.34 K 10^ km; 
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now 

OS - VTOE)- + (SE)2 ^ 2(OE) (SE) cos ^OES 

- V(6.38 X 10^')= + (L34 x lO"*)^ - (2) (6.38 x 10') (1.34 x 10^)cos 100^ 



^ 10'V40.70 + 119.68 + 24.24 
^ 10V249.62 ^ 1,58 X 10* km. 

So 

SP ^ 15800 - 6400 - 9400 km. 

Although the Sun is mora than a hundred times as large as Earth, as we noted in 
the first problem of Chapter 4, it subtends an angle of only about half a degree in 
the sky as viewed from Earth. In the next problem, we consider some aspects of 
the observation of sunspots. 

PROBLEM 6. a. Find the angular separation between two large sunspots when viewed from Earth 
(or Earth orbit) if they are separated by 30* In longitude along the Sun's equator. 
Consider two cases: 

1. A time when the midpoint between the spots is on the center of the visible disc 
of the Sun; 

2* A time about a week later when the Sun has rotated so that the leading spot is 
just about to go over the Sun's limb (edge). 

Recall that the Earth^Sun distance is 1.5 x 10^ km. The radius of the Sun is 
7,0 X 10^ km. In the first case it will help, and in the second it will be necessary, 
to make a suitable approximation (Fig. 7.8). 




Solutioni Case !• In the adge^on drawing shown above wa have: 

CD ^ Earth-Sun distance ^ 1,5 x 10^ km 
CA ^ CB ^ radius of Sun ^ 7.0 x 10^ km 
^ABC ^ 30^ and CD bisects ZACB 
Let AE be the perpendicular from A to CD and let h be its length. 
Then 

CE ^ h cot Z.ACE; 
ED ^ ^cot ZADE; 



and 



So 



h ^ CA sin Z.ACE. 

CD ^ CE + ED. 

^ h cot ^ACE + h cot ^ADE 

^ CA sin ^ACE cot ^ACE + CA sin ^ACE cot ^ADE. 



Then 



cot ^ADE ^ CD - CA sin ^ACE cot ^ACE ^ CD - CA cos Z.ACE 
CA sin ^ACE CA sin Z.ACE 

^ 1.5 X IQg - 7.0 X 10^ cos 15- 
7.0 X 10^ sin 15^ 

^ 1,5 X IQS 7.0 X IQg X 0.97 
7.0 X 10^ X 0.26 

^ L5 X IQg " 6.8 X 10^ _ T5 x 10» 



1.8 X 10^ 1.8 X 10^ 



^ 0.82 X 103. 



Z.ADE ^ arccot 0.82 x 10^ - 0.070®, so the angular separation between the 
sunspots ^ 2 ^ADE ^ 0,14®. A simpler solution can be found if we approximate 
AD by saying AD ^ CD. Now we can use the law of sines: 

CA AD 



sin ^ADC sin ^ACD ' 
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than 



and so 



CA PA 

sin ^ADC - ~ sin ^ACD - ~ sin ZACD 
AD CD 



7 0 X 10^ 
" 15 X IQg " "^'^ ^ ^ 

^ 1.2 X 10^3. 
Z.ADC ^ 0.07^ and ^ADB ^ 2 Z.ADC ^ 0.14-. 



Case 2, There is more than one way to solve this, but we present just one 
solution and use an approximation. In the drawing for Case 2 (see Fig. 7.8), 
construct the perpendicular AE from A to BC. For the approximation, we shall 

use tan ^ADB ^ ^ ^ In AAEC, cos ^ ACE ^ ^ , so CE ^ CA cos ^ACE ^ 

7.0 X 10^ cos 30^ ^ 7.0 x 10^ x 0.87 ^ 6.1 x 10^ km. 

Then EB ^ CB - CE ^ 7.0 x 10^ - 6.1 x 10- ^ 0.9 x 10^ ^ 9 x 10^ km. 

9 X 10^ 

Now, from our approximation, tan ^ADB ^ g - 6 x 10"^ giving the 

angular separation ^ADB - 0.036^. 

b. The unaided eye can distinguish a sunspot if it is 1.5 minutes of arc, or 0.025 
degrees, across. Sunspot sizes are usually measured in units of 0.001 of the Sun's 
area. What is the minimum size of sunspot that can be seen without a telescope? 

Solution: Sun*s area ^ 4w2 ^ 47r(7.0 X 10^- km^ ^ 1967r x lO^^ ]^^2 gince 1 sunspot 
unit ^ 10^^ of the Sun*s area, we have 1 sunspot unit - 196ir x lO'^ km^ 
Now if we assume that we have a sunspot that is approximately a disc sub- 
tending an angle of 0.025^ at Earth, we see from Fig. 7.9 that the disc has 

/n n25-\ 

radius = (Earth-Sun distance) x sin|^ j 

^ 1.5 X 10^ X 2.2 X 10-^ km ^ 3,3 x 10^ km. 



■Earth'Sun Lin€ 



Sunspot 
Dmmgten 



Fig. 7S 
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Such a disc has area 77(3 J X 10*)^ ktn^ 

/a w 1AR i ^ w 1 sunspot unit 
■ - 7T X 196 X 10' km^ 

_ lOg.9 X 10' 

^ 196 X ir'^"'P°' 



— 0.55 sunspot units. 



Historical note: Very faw sunspots exceed an angular diameter of 1.5 minutes of 
arc. Normally, the Sun is too dazzling to permit an observation of such a sunspot 
by the unaided eye; however, if the Sun is low on the horizon and shines through 
a thick haze, sunspots can be obsorved. Pretelescoplc sunspot observations have 
been recorded by Chinese and Japanese viewers. (Caut^uni Nevar look directly at 
the Sun.) 



The photographic scale factor for vertical aerial photographs was developed in 
Problem 7 of Chapter 4. We now consider the situation when the camera is tilted 
so that the film is not parallel to the ground. The result of such tilting is shown in 
Fig, 7.10, where the broken lines represent a square grid as it would appear in a 
vertical photograph and the solid lines show the actual image on a tilted photo- 
graph. (This is sometimes called the ^'keystone effect,*') In order to use the 
photograph to produce an undistorted picture, numerical relationships must be 
established between the actual shapes and their photographic images. 
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Scale variation on tilted photograph. 
Fig. 7.10 



Axil of tilt 



Fig* 7.11 
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PROBLEM 7* Fig. 7.11 shows the geometry of the configuration, where the camera is located at C, 
N is the nadir, V is the photographic nadir point, P is the image of ground point 
A, and / is the tilt angle of the film (the acute angle made by the film with the 
horizontal). If CT is the normal from the camera to the film, CT ^ /, the focal 
length of the camera. CN - H is the height of the camera above the ground, wnich 
we assume to be leveL It is customary in this work to use the film "positive" QW 
instead of the "negative'' PV, This is obtained by choosing W on CN and Q on CA 
so that CW ^ CV and CQ ^ CP. We let R be the point on the film positive so 
that CR is normal to the film. 
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Let 0 be the angle made by the line froni the caniera to the ground poitit A with 
respect to the vertical, where 0 > 0 If Q and R are on the same side of CN and 
# < 0 if they are on opposite sides of CN {$ ^ Z.NCQ). Express the ratio of the 
length of the image WQ to the length of NA in terms of $, t^fj and Hfor the case 
where B > t. 

Solutions From the geometry, wo ^ne that since CR ± WQ, Z.NCR - ^VCT ^ t and 

QVs QR + RW ^ / tan (g - t) + f tan t 
AN " an H tan 0 

/ (tan e - tan f) _^ . , , 

^ : + / tan t 

_ 1 + tan g tan t ^ 

H tan 9 

_ / tan g - / tan f + / tan f + / tan 6 tan^ t 
^ tan # (1 + tan 0 tan /) 

_ / tan fl (1 + tan^ t) 



H tan a (1 + tan & tan t) 



QW ^ /(I + tan' r) 
AN H(l + tan # tan r) ' 

b. Show that if f = 0 (untilted camera) or f ^ ^ (camera aimed at point A), then 
QW f 

- — — — ^ . (Recall from Problem 7 of Chapter 4 that this is the scale factor of a 
AN H 

vertical photograph.) 

Solution: For ^ — 0, tan f = 0 and the result fallows. For t ^ the second factor in the 
denominator becomes (1 + tan^ f), which cancels, and the result follows. 

c* Show that the result of part (a) still applies for the cases where Q is between 
W and R (Fig. 7.12) and where Q is on the side of CN that does not contain R 
(Fig, 7.13), taking into consideration the sign of 0, 
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Solutii nt h J is between W and R, then 



QW ^ RW - RQ ^ / tan t - / tan (t - $) 
AN AN H tan 6 



Solution: 



^ / (1 tang r) 
^ (1 + tan e tan /) ^ 

If Q is as shown in Fig. 7,13, then 6 is negative, and the positive value of the 
angle in the diagram is (-0). 



So 



QW „ QR - RW /tan ( (-g) + r ) - / tan f 
AN AN H 



/ tan U - B) -f tan t 
-H tan d 



f (1 + tan^O 



H {I + tan t tan B) 



Note that this implies (since e is negative in the last case) that points on the 
**down" side of the fi^m will have their images ^^stretched out/' whereas points 
on the ^'up'' side (at least those for which e > t) will have their images **shrunk." 

d. Fig. 7.10 showed a point I (called the isocenter) at which there is no distortion 
in the scale of the tilted photograph. Show that I is the point of intersection of 
the bisector of Z,NCR in Fig. 7.11 with the film positive QW by establishing that a 
vertical photograph taken with the camera in its position at C would contain ine 
point I. 

The bisector of ^NCR is shown in Fig, 7.14. along with a horizontal through I that 
intersects CN st Y. Since CN is vertical and lY is horizontal, ^CYI is a right 
angle, friangles CYI and CRI have corresponding angles equal and share side CI 
and are therefoie congruent. Since CY - CR - / a verticaf photographic posi^ 
tive and the actual photograph positive from the same camera position C both 
contain the point I. 
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Fig. 7.14^ 



FIm. 7.15 
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In Problem 8 of Chapter 4 wa developed a formula to correct an aerial photograph 
for distortioii due to Earth's curvature. The distortion occurs because the camera 
cannot distinguish how far away an object is— it has no "depth perception." For 
aerial photography^ the picture is interpreted as though everything is in the plane 
tangent to Earth at the nadir; in satellite photography, as we shall see in Problem 9, 
pictures will be interpreted (unless corrected) as though everything is in the 
horizon plane sensed by tha satellite, 

PROBLEM 8. Depth perception in humans has two aspects, called monoscopic and stereoscopic. 

Monosconic judgments of distance use only one eye and are based on an inter- ^ 
pretation of relative sizes of objects, shadows, hidden portions of abjects, and 
other attributes of this type; such judgments are very rough and frequently fail. 
Stereoscopic judgments of distance use both eyes and are quite accurate in most 
people. Stereoscopic judgment depends on the physical separation of the eyes, 
which causes an object to be viewed at a different angle by each eye, as shown in 
Fig. 7.15= The angle subtended by the ''eye base'* (the distance LR where L is 
the le; : eye ar : R the right) at the object O is called the parallactic anglei it is 
evident that the closer the object, the larger the parallactic angle. 

The smallest parallactic angle discernible by human eyes is about 0.025^, and the 
average adult eyes are spaced about 6.5 cm apart. What is the largest distance at 
which the average adult can judge depth? 

Solutions Let d be the distance of O from LR in Fig. 7 15. We present two methods of solution. 
The first uses the fact that 

tan (I zLOr) ^ (I Lr) /d,sod^ (0.0325)/tan(0,0125^) m 

- 150 m. 

For another approach, we may approximate LR as an are of a circle with radius d 
vv^here LR subtends an angle 6 = 0.025^, If $ is in radians, then LR - 0 ^ So 

e ^ 0.025'* ^ ^ 0.00044 rad and therefore 

. _ 0.025 „ ^ isn™ 



Satellites such as the Landsats, Seasat, and the Synchronous Meteorological Satellites 
(SMS-1 and -2) have made it possible to study Earth and its oceans, resources, 
and weather patterns as never before. They have returned observations and data 
that are being used by botanists, geologists, oceanogr aphers, and meteorologists, 
among others, in numerous projects. To ate just two examples, Landsat observa- 
tions have been used in the assessment of soil moisture in agricultural fields, and 
SMS observations have been useful in predicting severe storms. 
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PROBLEM 9. 



A ^acecraft at a distance h from Earth in synchronous orbit can sae only a portion 
of Earth s surface, as illustrated in Fig. 7.16. The circle that is the boundary of 
this spherical ''cap- will be called the horizon circle, and the spacecraft has^ensors 
that can recognize this horizon. 

Although every spacecraft uses its horizon sensors to find its angular direction 
with respect to Earth's center, those satellites whose purpose is to observe Earth 
can also use this angle measurement to determine the size of the spherical cap 
that can be observed. - 

In Fig. 7.17, S is the position of the spacecraft, C Is the canter of Earth H is a 
point on the horizon circle seen by the spacecraft, P Is the subsatallita point 
on Earth (the intersection of Earth^s surface with the Una from Earth's center to 
the satellite), and Q is the center of the horizon circle. We have SP - A and 
CH^ CP - r, the radius of E^Tth. p is the angular separation of tha horizon seen 
by the spacecraft from Earth's canter, and A is the angle subtended at Earth's 
canter by the radius of the horizon circle. 



Fig* 7.16 




SolutJoni 
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a. Find the relationships among p. A, h, and r. 
Since 4SHC has a right angle at H, sin p - cos A ^ 



r + h 



b. Listed below are soma Earth-observing satellites and their parigea or apog( 
distances from Earth. For each, find the angular radius (A) of *he horizon cird 
seen by the spacecraft. (Earth's radius is 6378 km.) 



ea 
e 



Landsat 2 
Seasat 



916 km (apogee) OGO-1 
790 km (apogee) OGO-1 
36 000 km (apogee) 



260 km (perigee) 
150 000 km (apogee) 



(OGO is the Orbiting Geophysical Omervatory.) 
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Solutioni Using A ^ cos' 



for Landsatl we have 



cos 



Similarly^ we get angular radii of 27^, 16^, and %T for the remaining 
casaSj respectively. 

If a satellite sees a horizon circle of angular radius 30% what is its distance 
from Earth? 



Solutioni 



cos 30^ ^ ^ -r - 0.8660 



6378 + h 

^ 6378(1 - 0.8660) 
0.8660 



^ 987 km (to the nearest 



In observing Earth froni space using spacecraft sensors^ distortions are introduced 
because of Earth's spherical shape. For example, suppose a thick black line is 
painted along the equator, the 10^ parallel of latitude, and the 50^ and 90^ west 
meridians of longitude as shown in Fig. 7.18(a). Uncorrected obser. jtiorisofthis 
"rectangle*' would appear as shown in Fig. 7.18(b). The diagram in Fig. 7.19 illus- 
trates how this distortion comes about. Although spacecraft sensors can ineasure 
the angle at which point R on Earth is observed , they cannot measure the (listance 
to R — all observations are interpreted as though lying in the same plane, so the 
image of R is treated as though it were at R '5 in the plane of the horizon circle, 
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The data can be corrected by the spacecraft's computers so that the information 
relayed to Earth is distortion-free. The actual computer program that does the 
correction depends aUso on the particular hardware of the sensors, but the first step 
in the correction is to express the relationship among the angle of observation of 
R the angular Alleviation of R from the line joining Earth's center to the 
sateUite (£.0), and th« angle of observation of the horizon (£p\ Since a and p 
can be measured, the computer can then find 0 for the proper mapping of R. 

d. Show that the relatzionship linking a, 0, and pis given by 

sin p iin 0 

tan a = — ^ _ 

1 ^ sin p cos 0 

Solutioni In Fig, 7.19, if T is thte foot of the perpendicular from R to CS then 

Un^a^Bl^ RT ^ rsinfl 

TS CS - CT (r + /i) - rcos^ 

sm p 



^ r _ _ _ sin p sin 0 

1 — sin p cos 0 




Fig. 7*19 



If the spacecraft sensors measure Zp as 30^ and a point R is observed at an 
angle of 25^ from the su satellite point, what is the actual angular displacement of 
R from the subsatelHte _tJOint with respect to Earth*s center? 

Solutioni We have z.p ^ 30% ^ 25% and we are seeking £0, 

From the last equation, tan 25^ * sinSO^smff 

1 - sin 30* cos 0 ' 



so 



0.446 ^ 



_ (0.5) sin 0 
1 - (0.5) cos^ 



_ sin B . 
2 - cos ^ * 
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then 5j-,2 ^ 

(0.466)^ - — 



4 - 4 cos ^ H- cos^ ^ ' 

0.217^ 2. ' 

4 4cos /3 + cos )3 



Clearing fractions gives 0.868 — 0.868 cos^ ^ 0.217 cos^ p ^ 1 — cos^ and col= 
iecting terms gives 1.217 cos^ /3 - 0.868 cos p — 0.132 ^ 0. 



Then 

^ 0.868 ± V(--OJ68)^ — 4(1.217) (-0.132) 
^ ^ ^ 2(L217) 

^ 0.868 ± VL396 
2.438 ■ 



Sine© we know that | ^ | < 90^, we discard the negative root, and so 



p ^ 32.8* ^ 33^ 

We have already seen in Fig, 2.2 of Chapter 2 that the celestial coordinate systefti 
uses a^ngies of dedination and right asceniion in a manner analogous to the 
latitude and longitude angles of the coordinate system of Earth. We now compare 
the three-dimensional spherical coordinate system commonly used in mathe- 
matics with the one generally used in astronomy and space science. 



PROBLEM 10» Texts in analytic geometry or calculus with analytic geometry usually define a 
spherical coordinate system so that if forP(p5 6^ ^) we let Q be the foot of the 
perpendicular from P to the Jf-y plane (Fig, 7*20)5 then 

p ^ the distance OP, p a 0 

9 — the angle made by OQ with the positive x-axis, the positive 
angular direction being a rotation from OX toward OY* 
0 m e Kl'n 

^ ^ the angle made by OP with the z-axis, with the positive angular 
direction being away from OZ. 0 ^ ^ ^ 




EKLC 



Chapter Seven 



ERIC 



In this syitem, as the reader may verify, the transformation between (x, y, z) and 
(p, 0, ^) are as follows: v _f , / - 

y ^ p sin ^ sin 0 a - (arctan (y /^) 

z — pcos 0 ^ ^ arccos (^/(x- + y^ + 

In the spherical coordinate system used by astronomers and space scientists if P 
has coordmates (r, 8, a) and Q is the foot of the perpendicular from P to the'x-v 
plane (Fig. 7.21), then ■ 

r — the distance OP, r a 0 

5 ^ the angle made by OP with OQ, the positive angular direction 
being from OQ toward the positive z-axis, 

a ^ the angle made by OQ with the positive jc-a?iis, the positive 
angular direction being a rotation from OX toward OY 

0 s a < 2^. 

Develop the transformations from (r, 8, a) to (x, y, z). 
Solutloni From the definitions, it is evident that 



p ^ d ^ a, ^ ^ 2 ^ a. 



So we have 



X ^ p sin ^ cos ^ ^ r sin - 5j cos a ^ r cos 8 cos a 
y ^ psin 0sin 0 ^ rsin - ij sin a ^ rcos fisin a 

^ ^ p cos ^ ^ r cos - a] ^ r sin 8 

(Recall that 8 is the daclination and a the right ascension in the celestial 
coordinate system.) 
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PROBLEM 11, On March 5, 1979, the spacecraft Voyager 1 passed close to the Jovian mOQn lo. This 
close encounter took place just after Voyager's closest radial appmach to Jupiter, 
which occurred at about noon on that day. If we set up a Cartesian coordinate 
system centered at Jupiter with the x-y plane as lo's orbital plane and the 
Jupiter-to-Sun vector as the positives-axis (see Fig, 7,22), then Voyager's 
spherical coordinates at fn ^ 13 hours were r ^ 5,0, S = ^5.0°. and a = 127°, 
(We measure lengths in units of Jovian radii, ^j, where 1 = 70 000 km. 
The spherical coordinate system used here is the one defined in the 
previous problem,) 



X 




Y 

Fig. 7.22 



a. What was the Voyager's radial distance from Jupiter in km at /n — 13 hours? 
Solution: r = 5.0 (^j) ^ 5.0 x 70 000 km ^ 350 000 km 

b. What were its Cartesian (x, y, z) coordinates in the system defined above? 
Solutioni X ^ 5,0 cos (=5,0^) cos (127^) ^ -3,0 

y - 5,0 cos ("5,0^) sin (127^) = 4,0 Rj 
z ^ 5,0 sin (-5,0^) ^ =0,44 
Three hours latef , at /i = 16 liours, its coordinates were 
r ^ 6,5, 8 - -^1,9% a ^ 166^ 
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c. What were V1o>?ager's Cartesian coordinates at /i ^ 16 hours? 

A' ^ 6.5 cos (=-1.9°) cos (166^) - -63 Rj 

y ^ 6.5 cos (-Lr) sin (166') ^ 1,6 

z ^ 6,5 sin (-1.9^) ^ -0.22 

As you can see, in the interval Voyager hts moved away from Jupiter in the anti- 
Sun direction (itsx-coordinate has become mora negative), toward the Sun- 
Jupiter line (its j-coordihate has decreased), and it has moved toward lo's orbital 
plane (its z -coordinate has decreased in abolute value). 

If we assume that Voyager's Cartesian coordinates change linearly with time 
between /o and fj, this means that we assume that Fojagerhas constant velocity 
components in the a% y, and z directions. 

d. Under this assumption, what are Voyager's velocity components in the x, y, 
and z directions? 



Solutions 



Solutionp 



Soluti 



xi - xo _ -6.3 - (-3,0) 



= 1,1 RJh 



ti - to 16 = 13 

LI X 70 000km/h ^ 77 000 km/h 
yi - yo _ 1,6 - (4.0) 



ti = to 

Zl ^ Zq 



_ 23 ^ "0.8 Rj/h ^ 56 000 km/h 
0,22 - (0.44) 



16 - 13 



0.073 ;?j/h - 5 100 km/h 



a. Under the assumption that Voyager's Cartesian coordinates vary linearly with 
time, find expressions for ^(f), >>(/), and z(r). 

Since ^ ^ Xa + V^ (t - fa), we have;c ^ -3,0 - 1.1 (/ = to). Similarly, 
y ^ 4,0 - 0,8 (t - to) and z - =0,44 + 0,073 (/ - fo). While Voyager was 
moving, lo had been progressing in its orbit. Consider lo's orbit to be a circle of 
radius r - 5.9 and recall that in this coordinate system, lo's 8 equals 0 at all 
times. At tn ^ 13 hours on 5 March, lo's phase n igle a was 139^. 

This phase angle is a linear function of time. Knowing that lo's orbital period is 
42.5 hours (i,e., it takes lo 42.5 hours to move 360^ in a), derive an expression 
fora(/). ^ 

360 

lo moves through |— — degrees per hour, and a = 139° at /„, so 

360 (/ - fo) , 
"^'^ " 42.5 + degrees = 8.47 (f - h) + 139 degrees. 
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g. Find lo's rectangular coordinates as functions of time. 
Solution: X = rcos 5 cos a = 5.9 cos (8.47 (i — to) + 139), since cos 5 = 1 

>' ^ r cos 5 sin a - 5,9 sin (8,47 (f - to) + 139) 
z ^ r sin 6 ^ 0 

hp Derive an expression for the separation distance A between Voyager and lo as a 
function of time. Use J for (/ = to). 

Solution: = (Xvay - Xin)' + C>Voy - yuO^ + C^Voy - ^In)" 

Froni parts (e) and (g), 

A- - [-3,0 = LU " 5,9 cos (8.47? + 139)p 
+ [4,0 - 0.8£ - 5.9 sin (8.47? + 139)]^ 
+ [=0,44 + 0,073£]^ 

1, Use a calculator and evaluate A, for several valuer Df in the interval 0 < { s 3. 
Plot the results, and use the resulting graph to find when Voyager's closest 
approach to lo occurs and at what distance. 

Solution: 



5 A 

0 LSI 

0.5 ■ 1.23 

1.0 0.94 

1.5 0,62 

2.0 0.33 

2.1 0,28 

2.2 0.26 

2.3 0.25 

2.4 0,27 

2.5 0.31 
2.75 0.46 
3.0 0.65 



L5h 



1,0 



A 



0.5 



J _L 



Fig. 7.23 



The graph is shown in Fig. 7.23. 

We see that the closest approach occurred at f = 13 + 2.3 hours ^ 15.3 hours on 
March 5 at r distance of 0*25 Rj, or about 17 500 km. 
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Chapter Eight 



latrices are an invaluable tool in space science, making it posiifale to organize, 
handle, and uianipulate, with the aid of computeis, large qunn lities of 
I data. Most of the actual examples involving matrix algebra are^ too long and 
coinplex for inclusion here; however, by considering simplified examples, w can 
get some sense of the role of matrix algebra in this context. 




PROBLEM L In Chapter 5 we considered some simple error-detecting binary codes rfor telenietry. 

A more complex system, the Hamming Code, will not only detect the presence 
of an error in a received message but will identify the erroneous bit in cases vvhera 
a sixigle error has occurred. If two bits are wrong, this fact wlU be detracted but 
the locations of the errors will not be known. We use a very simple exrample to 
illustrate the method. 

Suppose we have a -^message'' in the form of a four-bit binary strings chat is, the 
message is in the form abed where each of a, c, d is 0 or 1. The Hazmming matrix 
for a message of this type is the 4x8 matrix H: 



H 



0 0 0 0 1 1 1 1 

0 0 1 1 0 0 1 1 

0 10 10 10 1 

.1111111 1 



The structure of the matrix is as follows: For a message containing 4 bmts, we 
need 2 - = 8 columns and 4 rows. The binary numerals for 0 through?^ , (Nviitten in 
3-di jit form as 000, 001 , 010, . , = ) are used, in order, as the first three * efitries in 
each column; the bottom entry is always 1. A Hamming matrix for a5— bit message 
would need 2"* = 16 columns and 5 rows in order to represent the binao-y numer- 
als for 0 through 15 (0000, 0001, nil) followed by 1 in the columxjii. 



If the message we wish to send is abed, we need to use four additional iparity bits, 
Pu pi, pi,, and p4, and form a message row vector M — [p\ pj a p^b ccG], The 

^o' 



parity bits must be assigned so that the product H - 



in modSarithrnetic. 



a. Find the conditions that p, , p^, p^, and p^ must satisfy so that H • ^ 
mod 2 arithmetic. 



m 
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P2 
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Solutions 



So the conditions are Pa + b + c + d = 0', pi^ + a ^ c ^ d = Oi 
p2 + + d + ^ 0;/7i + p2 + + G + p4 + d + c + d - 0, 

b. Find the message row vector if the actual message is 0 1 10, 

We have a ^ 0, b - l,c = l.d - 0. Substituting these values in the preceding 
equations in part (a) and solving in mod 2 gives p^ — 0, pi ^ l,p2 = 1 , pj ^ 0, 
The message row vector is then M — [0 1 10 0 1 1 0]\ 

c. The matrix H ^ is a column vector called the syndrome vecior S\ In the set- 
ting we are using, 5 will have four components. When a message is received, the 
syndrome vector is formed. If none of the bits of M was in error, the components 
of 5 will all be 0, If we find that 54 = 1 , we know that an error has occurred in 
transmission, and the binary number JiS^s^ gives the number of the component of 
M which is wrong, where the components are numbered from the left, beginning 
with 0, If 54 ^ 0 and one or more of Su is 1, then two bits of M are incorfect, 
but we do not know which two^ — ^the error is detectable but uncorrectable. If 
there are more than two errors, it is possible that they will be "corrected'' incor- 
rectly or not get detected. 

Suppose the message [0 0 1 1 1 0 0 0] is received. Compute 5 and, if appropri- 
ate, correct the message. 
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PROBLEM 2. 



Since s-4 ^ 1, there is an error; si 52^3 = lOlbaso 2 = 5, so the error is in position #5 
(recall that the first position is #0) and the corrected message is [0 0 1 1 1 1 0 0], 

In the last chapter we developed transformations from a spacecraft coordinate 
system to a reference system with the same origin when the spacecraft has per- 
formed a roil or a pitch or a yaw rotation. Matrix algebra is the natural tool to use 
to find the transformation in cases where the spacecraft performs a series of such 
rotations. This is developed in the next problem. 

Recall that in Problem 1 of Chapter 7, we showed that 



>'r — y cos R -f z sin ^ 
Zh — z COS R ~ y sin R 



xp = X cos P — z sin P 

yp ^ y 

zp ^ z COS P + ^ sin P 



xy — X cos Y + y s\n Y 
yy ^ >^ cos K = ;r sin y 



EKLC 



where the uppercase R, P, Kare the angles of roll, pitch, and yaw respectively, the 
coordinates (x,y,z) are those of the reference system, and the subscripted coordi- 
nates are those of the spacecraft coordinate system after performance of the rota- 
tion designated by the subscript. 
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X 






- M - 


y 


where the sub- 






2 





a* Express these transformations in matrix form 



script sc designates the spacecraft coordinate system, by finding Mr, Mp, My, the 
matrices of roll, pitch, and yaw, respectively. 



Solution! Expressing each set of transformations above in matrix form, 



0 0 

cos R sin R 
^sin R cos R 



cos P 0 =sin P 

0 1 0 
sin P 0 cos P_ 



Mv - 



cos Y sin Y 0 
-sin Y cos y 0 
0 0 1 



If the spacecraft and reference systems are initially concurrent and the space- 
craft performs in sequence a roll through angle a pitch through angle P, and a 
yaw through angle F, then the transformation from reference system coordinates 
to spacecraft coordinates svill be given by 



where M ^ My ^ Mp A/r, 









X 






- M ^ 


y 










Find M if ^ - 30% 


P ^ 


45% Y 





60% 



Solution: 





"10 0 " 




1 


0 


0 




Mr - 


0 cos 30^ sin 30^ 




0 ^ 


/3 
2 


1 
2 






0 -sin 30^ cos 30% 




0 = 


1 

2 


3 

2 






°cos 45^ 0 -sin 45^^ 




2 


0 


2 


Mp ^ 


0 1 0 




0 


1 




0 




.sin 45^ 0 cos 45% 






0 


V2 








2 




2 . 








1 






0 




cos 60* sin 60^ O" 




2 




2 


Mv- 


-sin 60^ cos 60* 0 




2 




1 

2 


0 




0 0 L 




0 




0 


1. 
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M - My • Mp = Mr ^ 



1 

2 

2 
0 

1 

2 

2 
0 

4 

vs 

^'4 



0.35 
=0.61 
0.71 



2 

1 

2 

0 

2 
1 

7 



2 
0 

. 2 

2 
0 



0 -.^^ 



0 
\/2 



^2 
4" 

2 

'4 



"4 
1 

4 "J 



VI 3 

8 4 

V6 ^ ^ 
8 4 



0.93 0.13 
0.13 0.78 
=0,35 0.61 



8 4 



3 V2 
4 



1 
4 



0 


0 


v3 


1 


2 


2 


1 


V3 


1 


2 



Solutions 



c. The matrix M can be used to find the orientation of the spacecraft coordinate 
axes with respect to those of the reference system in terms of direction cosines. 

^1^ 

is a column matrix whose elements are the direction cosines of the 



X ^ M 



spacecraft ^-axis with respect to the x-,y-, and z-axes of the reference system. 



Similarly, Y - M ■ 



and Z ^ M 



produce column matrices whose ele- 



ments are the direction cosines of the spacecraft and z-axes, respectively, with 
respect to the reference system. For the motion of part (d), find X, Y, and Z, and, 
from these, the angles made by each of the spacecraft coordinate system axes 
with those of the reference system. 



^ - M 



The angles between the spacecraft x-axis and the;e-, j-, and z-axes of the refer- 
ence system are about 70'', -52*, and 45^, respectively. 



'l 




035 




cos 70^ 


0 




-0,61 




cos (-52*) 


.0. 




. 0.71. 




.cos 45^ 



f ^ M 







0.93^ 




cos 22^ 


I 




0.13 




cos 83^ 


.0. 




.=-0.35. 




.cos (-70*). 
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The angles between the spacecraft y-axis and the .r-, y-, and z-axes of the refer- 
ence system are about 22®, 83", and =70®, respectively. 



"o" 




"0.13^ 




^cos 83^^ 


0 




0.78 




cos 39® 






.0.61. 




,cos 52^ 



The angles between the sppcecraft z-axis and the x-,y-^ and z-axes of the refer- 
ence system are about 83^, 39®, and 52®, respectively. 

Analyzing the !?ght emitted from sources in space is a very important part of the 
astronomer's or space scientist's task* Some of these sources^ such as the stars, 
are too far away for their shapes to be discernible; but others are close enough for 
the emitting volume to be made out— that is, light can be seen to come from 
separate parts of the volume— and such sources are said to be "spatially resolved/' 
Among such sources are the solar Btmosphere, glowing at temperatures ranging 
from 250Q®C to well above a million degrees Celsius, depending on the particular 
location, and comet tails fluorescing under the Sun's radiation. 

If such a source is transparent to its own radiation— that is, light emitted at any 
point within it can escape from the source volume without being scattered or 
reabsorbed — then an observer looking at a particular area of the surface of the 
source will see the sum of all the light emitted behind that area, in the "line of 
sight/' The actual distribution of emitting intensity within the source, which in 
practice is always an u nknown function of position, is not directly available to an 
outside observer. 

However, when the source geometry is of an especially regular or simple shape, 
such as spherical or cylindrical, mathematical methods are available to '"invert" 
the observed intensity data, thereby "reconstructing" the source. 

In the next problem, we illustrate the basic idea with a very simple but concrete 
example in two dimensions, 

PROBLEM 3, Consider a small checkerboard, three squares on a side, on which a few lighted 
candles have been placed in some squares at random, as shown in Fig. 8. 1. If one 
looks down any row, the combined light of all the candles in that row will be seen; 
this combined light is simply the arithmetic sum of the separate candles in the 
row. Referring to the figure, if we look along row 1, the light of three candles will 
be seen; whereas, looking along row 3, we see the light of two candles. In similar 
fashion, one can look along a column or even along a diagonal. 
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Fig. S.l 
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To a ''two-dimensionar' observer in the plane of the checkerboard, this is, in fact, 
the only direct information available^ — the actual distribution of the candles on 
the checkerboard is unknown. This observer can, however, designate the number 
of candles in square (ij) as an unknown variable jCjj and proceed to set up a 
system of equations for these nine unknowns: 



Row equations 

X\\ "h Xi2 4- Xij, — 3 
^21 + ^22 -^23 ^ -3 
X31 + X32 + ^33 = 2 



Column^quatioris 
X\i ^ Xji 4" JC31 ^ 2 

^12 + ^22 + X^2 = 3 

xn + X27i + A'33 — 3 



Since we need nine equations to solve for nine unknowns, we may look along three 
of the diagonals to get 

^21 + X32 — 3 

■^11 -^22 ~^ Jf33 ^ 1 

-^12 + X23 = 3 

In this simple case, it is not difficult to solve the system by eh'mination; however, 
it is easy to see that a more general method is usually necessary, 

a. Write a matrix equation for this system of linear equations. 



Solution* 
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0 




^33 




3 



Use elementary row operations to find the source distribution for a case that 
produces the following matrix equation: 
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Although it is certainly possible to solve part (b) manually, it is no doubt obvious 
that a computer solution is more desirable even in this vastly simplified context. 
Any of the uommercially available programs to solve such matrix equations could 
be employed to produce the solution to part (b) or to discover that the solution 
to part (a) is not unique. 

In practice, the physical radiating sources encountered are more complex in 
several ways: (a) they are continuous distributions rather than discrete ones, as 
in the example just treated; (b) they are three-dimensional sources; (c) they do not 
have simple geometric shapes; and (d) distant (astronomical) sources cannot 
usually be observed from a sufficient number of directions to obtain a complete set 
of emission data. What this means is that each observation must be modeled as 
an integral rather than a simple sum and the integrals are generally complicated 
expressions that are difficult to ^invert/' However, such inversions can be car- 
ried out for certain types of local radiating sources. 

One recent example of this same technique in the medical field is Computer-Aided 
Tomography, or CAT scanning, in which X-ray radiation through a section of the 
human body is used to mathematically reconstruct a three-dimensional image of 
the section. For example, one kind of scanner measures the X-ray intensity that 
penetrates the portion of the body being imaged (such as the brain or the abdom- 
inal cavity). This scanner records the received radiation at 160 different positions 
in each scan direction; the entire unit is rotated one degree at a time around the 
head or abdomen, in a complete semicircle, to obtain 180 x 160, or 28 800, 
*'sums/' The computer then processes this information to produce a -'picture ' of a 
cross section of the organ by reconstructing the X-ray absorption in each square 
(or **pixer') of a 160 k 160 grid. The complexities enumerated in the foregoing 
paragraph also apply in this context, requiring the use of additional sophisticated 
mathematical techniques. However, the basic idea of the checkerboard model 
underlies this useful application. 
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The ApoHo 1 1 lunar module photographed 
from the command and service module dur- 
ing rendezvous In lunar orbit with Earth rising 
above the lunar horizon. 
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The conic sections play a fundamental role in space science. As shown in the 
AppendiXj any body under the influence of an inverse square law force (i.e,, 
where force is invefsely proprrtional ^o the square of distance) must have a 
trajectory that is one of the conic sections. In celestial mechanics the forces are 
gravitational; however, it is also of interest that the forces of attraction or 
repulsion between electrically charged particles obey an inverse square law, and 
^ such particles also have paths that are conic sections. 

Telescopes with mirrors that are conic sections are also important in space tech* 
nology because of their reflective properties. We shall close this chapter by con- 
sidering the design of an X-ray telescope that requires two reflections in sequence 
from surfaces whose cross sections are conies. 



In the analysis of orbits, where a celestial body, such as a planet, comet, meteor, 
star, or artificial satellite moves under gravitational attraction to a primary celes- 
tial body, the center of mass of the primary body is at one focus of the conic 
section along which the satellite moves. Because the simplest nontrivial conic sec- 
tion is the circle, we shall begin with a consideration of circular orbits. (The word 
"nontrivial" is included because a conic section could be a point or a pair of 
intersecting straight lines, if the sectioning plane passes through the cone's ver- 
tex.) Most of us understand from experience Newton's first law of motion, which 
states that an object in motion continues in a straight line unless it is acted on by 
some force. If we wish to make an object move in a circular path rather than in a 
straight line, we must give It a constant push toward the center. Thus a central, 
or centripetal, force is required. For example, when we tie a string to an object and 
whirl it in a circle, the pull of the string is the force that keeps the object in the 

2 

circular path. If we represent the centripetal force by Fu then Fi ^ , where m 

is the mass of the object, v is its speed or velocity, and r is the radius of 
the circle. 



When a spacecraft is moving in a cifcular orbit about any primary body, the force 
toward the center is supplied by the force of gravity Fj. According to Newton's 

law of universal gravitation, Fi — . In this equation, G is the constant of 

universal gravitation, assumed to be constant throughout the universe: M and m 
are the masses of any two bodies; and r is the distance between their centers of 
gravity. The physical situation, if the forces Fi and Fi are equal, is represented 
in Fig. 9.L 

The arrow toward the center represents the lorce of gravity, the dashed arrow 
represents the tangential velocity of the spacecraft, and the curved arrow indi- 
cates the circular path. (In rigorous use, velocity is a vector quantity, because it has 
both magnitude and direction, whereas speed, having magnitude only, is a scalar 
quantity. We will be using the symbol v for speed, the magnitude of the velocity 
vector.) Thus the force of gravity holds the body in the circular orbit, 

If we set Fi = Fi, we obtain = GMm §Q|^jj,g for y giyQs us 

r OB 
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This simple equation enables us to find circular orbital velocities about any pri- 
mary body, if M is the mass of the body and r is the radius of the orbit measured 
from the center of mass of the body. Because the varue of GM is constant for any 
primary body, it is convenient to substitute its numerical value rather than to 
compute the value of the product for each individual problem. If the primary body 
IS Earth, then CM — 3,99 x 10"^ m^/s-. Thus for bodies in circular orbits 
around Earth, 



^3.99 X 



m/s 



where, of course, the distance r is expressed in meters, 

PROBLEM L Most manned spacecraft in Earth orbit have been placed at altitudes of about 

160 km or more because atmospheric drag at altitudes below this causes a rather 
rapid deterioration of the orbit. Find the velocity needed for a body to stay in 
Earth orbit at an altitude of 160 km. 

Solution: Using the given equation. 



PROBLEM 2, 



/ 3 ,99 K 10'^ 
VEorth- V(63S0+ 160) X 10^""^' 



^7.81 X 10^ m/s, Of 2.81 x 10^ km/h. 

The formula for circular orbital velocity is quite general and can be applied to 
circular orbits about any primary body. C is a universal constant. We need only 
to change the value of M when we are concerned with another primary of different 

mass. 

a. The mass of the Moon is approximately 0,012 times the mass M of Earth, 
Write a formula for finding circular orbital velocities about the Moon. 



Solulioni Multiplying the numerator in the previous equation by 0,012, 

^ ^0,012 X fl 



Vmo 



.99 X 10*' 



m/s 



■4- 



4.8 X 10*^ 



m/s. 
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b. During the ApoHo flights the parking orbit for the command and service mod- 
ule about the Moon had an altitude of 110 km. The radius of the Moon is about 
1740 km. Find the velocity in this orbit. 



Solution: 



PROBLEM 3. 



Solution: 



4,8 X 10^^ 



m/s 



(1740 + 110) X 10^ 
^ lO^vTim/s ^ 1600 m/s, 
Of 5800 km/h. 

A synchronous Earth satellite is one that is placed in a west-to-east orbit over the 
equator at such an altitude that its period of revolution about Earth is 24 hours, 
the time for one rotation of Earth on its axis. Thus the orbital motion of the 
satellite is synchronized v^ith Earth's rotation, and the satellite appears, from 
Earth, to remain stationary over a point on Earth's surface below. Such commu- 
nication satellites as Syncom, Early Bird, Intelsat, and ATS are in synchronous 
orbits. Find the altitude and the velocity for a synchronous Earth satellite. 

The velocity can be found from the equation for circular orbital velocity. It can also 
be found by dividing the distance around the orbit by the time required; that is, 

^ ^ ' Because the two velocities are equal. 



l-nr [cm 



GM 



GMt" 

47T^ 



It is apparent that t = 24 hours = 86 400 seconds. Substituting the other values 
yields 



,/ 3.99x iO'''x (86 400)' , 
1 4 X (3.14)' = 10 

= 4.22xlO'ni, or 42 200 km 

Altitude = 42 200 - 6400 km = 35 800 km 

, „ 2x3.14x42 200 . 



24 



■= l.lOx 10*km/hr 
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!_ 




Fig. 9.3 



! 



To understand orbits, we must know something of the nature and pfopefties of the 
conic sections. They get their name, of course, from the fact that they can be 
formed by cutting or sectioning a complete right circular cone (of two nappes) with 
a plane. Any plane perpendicular to the axis of the cone cuts a section that is a 
circle. Incline the plane e bit, and the section formed is an ellipse. Tilt the plane 
still more until it is parallel to a ruling of the cone and the section is a parabola. 
Continue tilting until the plane is parallel to the axis and cuts both nappes, and the 
section is a hyperbola, a curve with two branches. It is apparent that closed 
Orbits are circles or ellipses. Open or escape orbits are parabolas or hyperbolas 
(see Fig. 9.2). 

Another way of classifying the conic sections is by means of their eccentricity. 
Let F be a fixed point (focus) and d a fixed line (directrix). For nonzero values of 
eccentricity e, a conic section may be defined as the locus of points such that the 
ratio of the distance PF to the distance from P to d is the constant e. The use of 
polar coordinates permits a unified treatment of the conic sections, and it is the 
polar coordinate equations of these curves that are used in celestial mechanics. 

PROBLEM 4. Use the eccentricity definition above to show that the equation of a conic section in 

polar coordinates can be stated as r = . ^ , where p is the distance 

*^ 1 ecos 6 ' ^ 

between F and d, and the polar axis is perpendicular to, and pointing away from d, 
with the pole at F as shown in Fig, 9,3, 

Solutions If Q is the foot of the perpendicular from P to d, and P has coordinates (r, &), then, 
by definition. 



PF 
PQ 



p + r cos 0 



r^ep + er cos 0 
r ~ er cos 0 = ep 
r(l = e cos 9) — ep 



e cos d ' 
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Fig. 9.4 shows a family of conies, each of which has directrix d and focus Fo, for 
different values of e. If a Cartesian coordinate system has origin Fn and ^-axis 
along the polar axis, the Cartesian equations of these conies have this form: 

ellipse: l£^ + r«i 

parabola: ^ q (x = h) 

hyperbola! — — ~- + ^ i 




Fig. 9.4 
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PROBLEM B. Shovv that the polar aquation of Problem 4 can be transformed into the Cartesian 
equation of an ellipse if 0 < e < 1 ; a parabola if e = 1 ; and a hyperbola if ^ > L 
Express the parameters h, a, b, or q, as appropriate, in terms of e and p. 



Solutions 



1 = e cos 6 
r ^ er cos 6 + ep 
If € — 1 J then r = er cos 0 + ep becomes r = r cos 6 + p. Since 

r cos ^ = ^ and r = Vx^ + , S/x^ + y^ ^ x + p. 

Squaring, 

^2 ^ y2 ^ ^3 2xp + p' 



125 



Conic Sections 



so 

/i = — ^ and q ^ 2p. 
If e # 1, then r ^ ercos 6 + ep becomes 

Dividing by (1 = e^) and completing the square, we get the following: 

X-- f + • 



\- e'l 1 - e- (1 - e-)- 



CP f e-p' 
1 - 1 _ 



If 0 < e < 1, the denominator of the term is positive, and we have an ellipse 

pe" ^P t _ 
with h = =^ ^ , a ^ — = , o ^ " 



1 ^ e^'" " 1 = e^- Vr^' 

If e > 1 , the denominator of the term is negative, so we may rewrite the 
equation as 



1, 



e- - 1 

and we have a hyperbola with 



I pg' gp . „ 



1 - e" - 1 ' " 
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PROBLEM r 



Recall that for an ellipse a= 6= - c- and for a hyperbola a' + ^ where in 
both cases, c is the distance between the center of the conic and a focus and a is ' 
the length of the semimajof axis. Show that the results of the preceding problem 
nre consistent with this and that in both cases e ^ da. 



Solution: For the ellipse. 



For the hyperbola. 



1 - 



c" ^ + b- 



o + ~ 



c = 



We see that for an ellipse, 2 ^ c is the distance between the foci. Since e ^ da \f 
c ^ 0 we have 0^0; but if c ^ 0, the two foci coincide with the center and we 
have a circle rather than an ellipse. A circle can therefore be considered the conic 
section with eccentricity 0. 

It is shown in the Appendi:t that the total energy £ of a two-body gravitational 
system and the eccentricity e of the orbit of the less massive body (mass m) with 
respect to the more massive body (mass M) are related by 

^ ^ GMmje^ - 1) 
2ep 

Since it is virtually impossible in the real world for the total energy to have a 
value that would result in e - 0 of e - 1 exactly, orbits that are exactly circles or 
exactly parabolas do not occur in nature. However, such orbits are of interest as 
limiting cases of actual trajectories. The energy equation of the Appendix 



1 



GMm _ ^ _ GMm(e^ ^1) 
r ^ 2ep 



provides the means to determine the velocity of an orbiting 
orbit, - 



at any point in its 



PROBLEM 7. 



Solutloni 
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in am orbit 
problem 5), for each 



Solve the energy equation for v, and then express the velocity at any point 
m terms of G, M, r and a, if needed (where a is defined as in problem 5\ 
type of orbit, 

1 _ ^ GMm GMm (e- ^ 1) 
- m V = — — + ^ * 

2 r lep 



v^^GM 
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For an ellipse. 



For a circle, r = constant — a 

(Recall that this was shown at the beginning of this chapter in the preliminary 
discussion of circular orbits.) 

For a parabola, e = 1, and 

For a hyperbola, 

' - , .0 that 

€p a 

The minimum escape v elocity of a rocket-borne space probe is the paraboilc 
velocity Vp = \^2GM I r , Velocities greater than this produce a hyperbolic orbit, 
and lesser velocities produce an elliptical orbit (or no orbit if too small). 

Elliptical orbits are frequently analyzed in terms of orbit parameters, such as 
apogee and perigee distances. These distances are indicated in Fig, 9.5 by the let- 
ters A and P respectively. Before we discuss elliptical orbits, It will be necessary 
for us to avoid ambiguity by clarifying our terminology and mathematical notation. 
Most of us know from our reading of space events that in NASA news reports the 
point in an orbit nearest the surface of Earth is called the perigee, whereas the 
farthest point from the surface is called the apogee. These points are indicated by 
C and D, respectively, in Fig, 9.5. In common usage the word is used to refer to 
either the position of the point or the distance to the point. 





Chapter Nine 



PROBLEM 8. 



However, usage is not uniform; some references state that the distances are 
measured, not from the surface of Earth, but from the center. In this article we 
shall use distances measured from the center. The distances from the center'to 
the perigee and the apogee will be indicated by P and A, raspectively. In most 
discussions, the context will make this clear. If in any situation confusion could 
result, then distances from the surface, if used, will be called perigee altitude or 
apogee altitude, whereas distances from the center will be called p^n^ee radius or 
apogee radius. Incidentally, the mathematics is simpler when distances are mea- 
sured from the center. 

a. Express the distances A and P in terms of the semimajor axis a and the eccen- 
tricity e of an ellipse. 



Solutions From Fig. 9.5, 

A = a + c— fl +ea = a (1 -¥ e) 
P ^ a — c ^ a — ea = a (1 — e ). 
b. Express the eccentricity of an elliptical orbit in terms of A and 
Solutions The following relationships are apparent from Fig. 9,5: 

a^i(A ~h P), 

c ^ a ^ P ^ ^(A ^ P) ^ P ^ i(A - P), 

and 

a 1{A P) 

A - P 

e ^ ^ ^ 

A ^ P 

This formula is a quick and easy way of finding the eccentricity of an elliptical 
orbit . As a check, we note by inspection that e ^ Q when A ^ P, which is the 
condition for a circular orbit. 

PROBLEM 9. Derive formulas for and v^, the velocities at apogee and perigee, in terms of A or 
P, respectively, and e of the elliptical trajectory. 

Solutioni From Problem 7, the velocity of a body in an elliptical orbit at a distance r from the 
focus is 



V = 
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If r ^ w4 ^ a (1 + e), we can substitute l/o = (1 + e)/A and r = A to get 
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If r — F = a (1 — e)i by a similar substitution. 




(1 + ^), 



PROBLEM lO, 
Solutions 



PROBLEM IL 
Solutioni 



These gquations can be written in other ways as well, because numerous ways of 
expressing relationships among e, a. A, and P are possible. The parlicular form 
for the formulas reflects personal preference. 

Show that the velocities at apogee and perigee are inversely proporlional to tha 
distancei from the center. 

If we divide the equation for Va by the equation for vp (see Problem 9), we obtain 



Thus the velocity at perigee is inversely proportional to P^ and so on, That is, when 
the orbital distance from the center of the primary body is smallj the velocity at 
that point is large; when the distance is large, the orbital velocity is small, This 
result agrees with Kepler's second law of planetary motion, which states that a 
planet moves about the Sun in such a way that the radius vector from Sun to planet 
sweeps out equal areas in equal timas. 

Derive a formula for the period of an elliptical orbitj given that the period of an 
elliptical orbit with semimajor axis a is the same as that for a circle with radius 



Following the method used in Problem 3, we express the velocity in terms of the 
distance around the orbit and the time p required to make one transit of the orbit 





Also 




Then 




jlitrf ^ GM 
P'~ ~r- 



^ GM 
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Because the period is the same when r - a, we may write 



= 27^ 
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PROBLEW12, An Eajrth satellite is placed in an elliptical orbit with perigee altitude of 160 km and 
apogee altitude of 16 000 km. Use 6380 km for the radius of Earth, 

a. If injection is at perigee, what must be the injection velocity? 
Solution! We firet find the eccentricity as follows: 

P - 6380 + 160 - 6540 km or 6,54 x 10^ m 

A - 6380 + 16 000 - 22 380 km or 2,24 x 10^ m 



By Problem 8, 



By Protilem 9, 



22 380-6540 ^ 15840 
22 380 + 6540 28 920 



- 9.72 X 103 n^/g^ 3 ^ ^^j^ 
b* Fin^ the spped at apogee, 

SolutiQii! By Problem 9, 

- V 2 24 x l0^ <^ ^ 0-55) - VL78x10^x< 



0,45 

- lO-'VOSm/i - 2,83 X 10^ m/s, or 1.02 x 10^ km/h. 
c. Find the period in this orbit, 
Solutioii! From Problem S, 

_ 22380 + 6540 ^^^^^^ 

^^14460 km, or 1.446 x 10^ m, 

and, frorri Problem 11, 



^ /(L446X 10^) 
V 3.99 xin'^ 



3 

s 



= (2it)(10') VTils 

= 17.2xl0's, or 4.77 h. 
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PROBLEM 13. During the Apollo flights, the Apollo spacecraft and the third stage (SIVB) of the 
Saturn V launch vehicle were placed in a parking orbit 190 km above Earth, Find 
the velocity and period in this orbit. 



Solution* Because r = 6380 + 190 - 6570 km, or 6,57 x 10^ m, we find from Problem 7, 

IO^VO.6073 m/s - 7.79 x lO^ m/s, or 2 J x 10^ km/h. 



PROBLEM 14. 



PROBLEM 15, 



/3,99 



X 10^^ 



X 10- 



From Problem 11, 



p ^ 2tt 



(6,57 X lOT 



3.99 X 10^ 



5300s-L47h, 



During the flight of Apollo 11, the SIVB stage was reignited and burned long enough 
to place the Apollo spacecraft on a trajectory to the Moon. At the end of the 
burn, the spacecraft had a velocity of about 3,90 x 10^ km_ per hour at an altitude of 
336 km. Was the Apollo spacecraft given escape velocity? 



Solution! Using the results of Problem 7, the escape velocity e 



llGM _ I 2(3.99x10^^) 
^^^i^ ^ V(6380 + 336)xl0'"^^^ 



^ lO^vTTi m/s ^ 1,09 X 10^ m/s, or 3,92 x 10^ fim/h. 

Thus the velocity imparted was about 200 km per hour less than escape velocity, 
thereby assuring a free return trajectory. That is, if the major propulsion systems 
failed, the spacecraft would be going slowly enough to be pulled around and ori- 
ented back toward Earth by lunar gravity, the attitude-control system being ade- 
quate to make needed course corrections. 

A spacecraft, as illustrated in Fig, 9.6, is in a circular orbit 800 km above Earth , The 
spacecraft must be transferred to a lower circular orbit 160 km above Earth, 
Computs the velocity changes needed at A and P to achieve this transfer. 



A 
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Solutioni We first find the eccemricity of the transfer orbit, which is. of course an elllnse 
with ^ -7180 km and/' = 6540 km. ' ' 

7180-6540 ^ 640 , ^ 
7180 + 6540 13720 

We then compare the velocities at A in the circular orbit and the elliptical orbit 
to fmd what chanies must be made. Since CM has units m'/s=. we express A and P 
in meters. From Problem 7, 



3.99 X 10" , , , , 

7.18 X 10* '"'5 = 7.45 X 10- m/s, or 2,68 x lO'' km/h. 



and from Problem 9, 



/3 99 X lO''' 
= y ^ (1 - 0,047) = 7,28 X IQJ m/s. 



or 2.62 X 10^ km/h. 



Therefore a propulsion engine on board the spacecraft must be fired long 
enough so that a retrothrust (opposUe to the direction of motion) will slow down 
the spacecraft by 600 km per hour. The spacecraft will then leave the 800-km 
circular orbit and follow the elliptical transfer orbit, remaining in It indefinitely 
unless additional changes in velocity are made. 

When *he spacecraft reaches the point P, however, we want it to move from the 
elliptical orbit into the 160 km circular orbit. Therefore we must use the results of 
Problems 7 and 9 to investigate velocity changes at P, 



/ 3.99X lO ''' , „, , 
V 6.54 X 10« *" ■ ^ °^ 2.81x10^ km/h 



/3.99 X 
' V 6.54 X 



10''* 

—y(l,047) = 7.99x 10' m/s, or 2.88 x 10^ km/h. 



10* 

That is, a retrothrust must reduce velocity again, this time by about 700 km/h. 

This method of transferring a spacecraft from one orbit to another is known as a 
Hohmann transfer, named after Walter Hohmann, city engineer of Essen Ger- 
many, who published the method in 1925. There are many paths that could be used 
to move the spacecraft from the 800 km to the 160 km orbit. But the Hohmann- 
transfer ellipse, requiring only two short burns, is the most economical, taking the 
mLjiimum amount of energy. Therefore this method is called a minimum-energy 
transfer. It has many applications. 
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PROBLEM 16^ A satellite is placed into a synchronous orbit by a technique involving a Hohmann* 
transfer ellipse. We computed in Probleni 3 that the altitude of such a satellite is 
about 35 800 km and its orbital speed is about 9370 km per hour. Fig. 9.7 suggests 
the details. 



Soluttoni 




Fii. 9.7 



We shall assume that injection is at the perigee point, which we shall place 
160 km above Earth. Then obviously 

P ^ 6380 + 160 ^ 6540 km, and 

A ^ 6380 + 35800 ^ 421S0 km. 

We wish to find the velocity change needed at A. 



e ^ 



42180 ° 6540 ^ 35640 
42180 + 6540 48720 



^ 0.732 



'4 



g^^^ (1.732) - 1.028 X 10^ m/s, 



or 3.7 X lO^km/h 



' ■ ■ - % (1 - 0.732) - 1.59 X lO^m/s, 



4.22 X 10^ 



or 5.73 X lO^km/h. 



But the tangential velocity needed at point A is 9370 km per hour. Therefore the 
velocity of the satallite must be increased in the direction of Earth's rotation by 
9370 5730 — 3640 km per hour. This extra push or kick would be provided by 
the firing of a motor on board the satellite, and the thrust and firing time must be 
such as to give the desired increment in velocity. Such a motor to be fired at 
apogee is called an apogee motor, and the thrust it provides is calli^d an 
apogee kick. 

The relative efficiency of using this method is easy to understand. Placing a heavy 
final stage of the launch vehicle at the synchronous altitude and then having a 
burn to give the entire assembly circular orbital velocity would take much fueL 
Instead we send up to the synchronous altitude only a Telative] ii^ht satellite 
and a small apogee motor. The numerical values used jw ihi^ pi.>uiem are merely 
illustrative. If th^^p4p|.ee altitude is higher or lower than the one we have 
assumed, all the other numbers are changed. 
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PROBLEM 17, 



One more maneuver is needed to make the satellite synchronous. It now has a 
period equal to the t me of Earth^s rotation. However, the satellite will appear to 
be stationary over a jiven point only if it is in equatorial orbit. Unless corrections 
were made during la.mch, the plane of the orbit will be inclined to the plane of 
the equator. One method of solving this problem is to fire a motor at the precise 
instant when the satellite crosses the equator, adjusting the burn time and direc- 
tion of thrust so that the vector sum of the burn velocity and the orbital velocity 
make the angle of inclination equal to zero. 

The first step in lunar orbit injection in the Apollo 11 flight was to place the 
spacecraft in an elliptical orbit of 110 by 313 km, the low poxnt^ot periUme 
(corresponding to perigee for Earth)— being on the back side of the Moon. 



a. Compute the ve 
this orbit. 



locity needed at perilune to inject the Apoilo spacecraft into 



Solution: Using the data developed for lunar orbits in Problem 2, 

P - 1740 + 110 - 1850 km 
A - 1740 + 314 - 2054 km 
2054 - 1850 204 



2054 + 1850 3904 



^ 0.052 



/4 8 X 10^^ 

V L85 K 10^ (1 052) IO^VTtI - 1,65 x 10^ m/s, 

or 5.95 X lO^km/h. 



b. Find the period in this orbit. 
Solution: Evidently a ^ j (2054 + 1850) km, or 1.95 x W m 



PROBLEM 18. 



Solutions 
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and 



(1.95 K 10^)-' 

^4.8 X 10^^ ^ ^ 2tr(1.24) x 10^ s ^ 782 s, or 130 m 

The lunar module descent orbit insertion during the Apollo 11 mission began with a 
Hohmann transfer. The command and service (CSM) and lunar modules were in 
a circular orbit 110 km above the Moon. The lunar module was detached and its 
descent engine was fired to reduce velocity so that it would enter a 110-by-15^km 
lunar orbit. Find the reduction in velocity needed to achieve this orbit. The CSM 
remained in the 110 km parking orbit. 

In this case, the change to the elliptical transfer orbit was made at apolune fcorre^ 
sponding to apogee for Earth). 



A - 1740 + 110 - 1850 km or 1.85 x 10^ 
P - 1750 + 15 ^ 1755 km or 1.755 x 10^ 



m 



m 
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1850 - 1755 
1850 + 1755 



3605 



95 



^ 0,026 




4.8 X 10'^ 
L85 X 10' 



(1 - 0,026) - ]0V2.526 ^ L59 k 10- m/s, 

or5700km/h. 




no km 



IS km 



Fig, 9,8 



We found in Problem 2 that the circular velocity in the llO-km orbit was 5800 km 
per hour. Thus the reduction in the velocity needed, achieved by a retroburn of 
the lunar module descent engine, was 100 km per hour. At perilune altitude of 
15 km, several retroburns and attitude changes were made—both automatically, 
and manually by the pilot— causing the spacecraft to descend to the surface. If for 
any reason the descent from the 15-kni perilune could not be made, the lunar 
module could have remained indefinitely in the elliptical transfer orbit until a 
rendezvous and docking with the CSM could be made. Thus this maneuver, 
which seemed so tricky and dangerous as we watched before our television sets, 
was actually a routine Hohmann transfer. The tricky maneuver, requiring some 
manual control, came when the powered descent to the lunar surface was made 
from the 15-km altitude. 

We will conclude our discussion of orbits by considering the classic analysis 
known as Kepler's Problem, which in modern times makes use of high-speed com- 
puters to produce final results. It is the task of determining the exact position of 
a body in an elliptical orbit at any given time, Kepler, of course, was interested in 
establishing the nature of the planetary orbits around the Sun, but today the 
same analysis is used to predict the location of artificial satellites in their orbits 
around Earth, 

We shall make use of a number of the relationships involving elliptical orbits 
already established. For the orbit illustrated in Figure 9.9, 



1 ^ e cos B 



1 = e cos B 



FP ^ 



P 



CP ^ CA = a 



CF ^ 
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whera F is one focus of the ellipse and the location of the primary body in the 
gravitational system; A, P, C are apogee, perigee, and center of the ellipse, 
respectively: ^nd s 15 the eccentricity of the ellipse, 

Kepler's Problem is stated in terms of theangle ij (called the true ano/naiy of the 
rlhpse) between the Earth-perigee ray (FP) and the radius vector (FS), rather 
than the angle e, as shown in Fig. 9,9, Since vh the supplement of the ellipse 
equation may be written in terms of as 

r = . — 

1 + ^ cos V ' 

Recall also that the rate at which the radius vector traces out the ellipse Is not 
constant, but is in accordance with Kepler's first law: The radius vector sweeps out 
equal areas in equal time. This makes the task of expressing rand indirectly in 
terms of time extremely difficult. Kepler circumvented the problem by considering 
the projection of the ellipse on an ^'auxiliary circle" having the same center and 
passing through P and A as shown in Fig, 9.10. If a satellite is at S on the ellipse 
and Q is the foot of the perpendicular from S to AP, then is the intersection of 
OS with the circle. Kepler defined three new quantities: the foreshortening factor 
^ T SQ/S^Q; the eccentric anomaly E ^ ^S^CF, and the mean anomaly M, 
which IS a fictitious angle through which an object would move at a uniform angu- 
lar speed with respect to F, That is, M - (4f /T) « 2it radians, where Tis the time 
for one complete orbit and At is the time of interest. He then established the 
following relationships: 

(1) k^VT^P 

(2) r ^ a(l - 6 cos E) 

(4) Af ^ E = e sin E 

These will be derived in Problem 20. 




Fii' 9,9 Fig. 9.10 
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In ofder to determine a satellite's position at any time, we must be able to compute 
r and u Since in general, e and r(and therefore M) are known for an orbit, if 
equation (4) can be solved for then (2) and (3) will provide r and v. But aqua- 
tion (4) is transcendental in E, so that no analytical solution is possible, This 
difficulty has been the core of many computational schemes generated by astrono- 
mers, mathematicians, and physicists. 

High-speed computers now make a numerical, iterative solution both possible 
and feasible. The iterations would proceed as follows: 

El ^ 0 

Ej — M + ^ sin El 
E3 ^ M + e sin E2 



Ek i- 1 — M + e sin 

The iteration continues as long as necessary to compute E to a desired accuracy 
(say 10""^ or 10""-); in other words, when jEi^ + i E^] < 10""^, we can use 
E ^ Ej, + I, if this is our desired level of accuracy. Since Ek + 1 = Ek ^ 
e (sin Ejj — sin Ei^ ^ 1), it can be shown that |Ek + 1 ^ Ekl m e^^ so that the 
sequence Ej^ converges to the new value of E, since e < 1, This process is highly 
efficient for small values of and after a few iterations, it is usually found that 
the difference is within tolerance. 



PROBLEM 19. a. Write a computer progream to perform the iteration outlined above, and then to 
use the value of E so found to compute r and 1^, where E found to an nr*^"^ .cy 
of 10^^^ 




b* Use this program to find the position of the satellite discussed in Problem 12 
one hour after it passes the perigee point in its orbit. 



Solutions We had e - 0.55, a ^ + p) ^ 1(22 380 + 6540) km ^ L45 x 10^ km. and 

T = AJl h. Running the program of part (a) with these values produces the 
following results: 

3RUN 

WHAT li 'THE ECCENmiCITV OF :TH6^ ORBIT? 
to * ' > - ^ ^ ^ ' 

WHAT li THi SEF1IHAJ0R AXIS LENeTHV^ 
IN KILOMETE^i? v j.. . 

-^la^oQ ^ ' 

WHAT IS THE PERIOD OF THE ORiITt - 
IN HOURS'' 

HDW MANy HOURS AFTER - PERIGEf^ IS THE 
POSITIQN ,T0 BE OfTfRMlNED? t " . . ' - v 
? I ' , ^ : 

12 ITERATIONS Unnm RIQUIRED TO 
ACHIEMI DESIRED ACCURAeV. , . 

THE satellite' Ii^les7e'KM^'^^^^/'^'^^^^^■ ' ■ 
DISTANT FROH lARTHf : ' .... 

AND ITS ANOnALY IS 2i3? RAOlANp*^ . : 



PROBLEM 20. The four relationships of Kepler's Problem can be established using the geometry 
and trigonometry of Fig. 9.10. 

a. Show that ^ vT^7 

Solution: k was defined as /c = SQ/S'Q. From Fig. 9, 10, we have 

SQ ^ r sin 

and 



S'Q - V(S'C)- = {CQf - Va^ = (CF = QF)^ 
= = (ae + rcos y)" . 



Then 



^ rsin V / r"(l — cos^yj 

Vfl' - (ae + rcos py y — (ae + rcos i^)^ ' 

(1 - e^) 

Substituting r ^ i + ^ ^os y ""^ simplifying produces, after some labor, the result 

k - Vl - e', 

b* If a rectangular coordinate system is placed on Fig. 9. 10 with origin at F and 
positive X-axis along the polar axis, express the rectangular coordinates of S in 
terms of E, a, and e. Then use the fact that S can also be given by (— r cos i>, r sin v) 
to show that r ^ a (1 — e cos E). 
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Solutions Since CF = aeandS'C = a, the. ^ ^coordinate of S is OF ^ CP = CQ ^ as - acosE. 
The y-coordinate is 

SQ = *S'Q 

= \ I ^ (a sin E). 
So the coordinates of S are (ae - a cos E, Vl - a sin E), 
Now ^ r^cos^i^ + r"sin"y — (=rcos v)^ 4- (rsin v)^ 

^ (ae - a cos E)= + (VT^7 s sin £)' 

= ah^ - 2a^^ cos E + ,g"cos^£ a^sin^E. - a^^a si^^^ 

— a^e^ = 2s3e cos E + — a^e^(l ~ coi^E) 

= " 2a cos E + a^e~ cos^E 

= fl^ (1 ^ 2^ cos E + cos'E) 

^ (1 ^ e cos E)", 
So r=afl = f cos E). 

c. Use (b) and the identity 

, 6 ll - cos 6 

tan = = \ = ~- 

2 V 1 + cos d 

to show that 



^ - 1^ _ IjtA * E 

Solutioni Since r ^ a (I — e cos E) and =r cos = a cos E), 

cos E ^ e 



cos — 



1 — e cos E 



Then 



and 



Then 



1 _ 1 — e cos E — cos E ^ e (1 + (T = cos E) 

1 ^ cos V ^ — - — — ^ : — ^ — - — 

1 ^ e cos E 1 ^ e cos E 



1 + cos y = 1 " e cos E 4- cos E - g ^ (1 - g)(l cos E) 
1 ^ e cos E X ^ e cos E 



V __ /l - COS y _ /l 4- g /l — cos E /l + e E 
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d. As the satellite moves in its orbit so that the radius vector sweeps out equal 
areas of the ellipse in equal times, FS' sweeps out equal areas of the auxiliary 
circle in equal times. If the area enclosed by FP, FS% and the arc S'P is swept out in 
time At, and T is the tims for the satellite to traverse the entire ellipse, then this 
area is given by Ai/T(.Ta% Recall that M is {At/T) (liv), and use the geometry of 
Fig, 9, 10 to show that M ^ E ^ e sin £, 

Solution: Area FPS' =- area of sector CSP' - area of CFS' 

At _ n _ 1 ^ ^ 1, , 




^ e sin £■ 



M ^ E ^ e sin E 

Among the first telescopes used to explore the heavens were those based on the 
reflective properties of paraboloidal mirrors. It is the fact that all light striking 
such a mirror in the direction parallel to the axis of the paraboloid is reflected to 
the focus that provides the light-gathering capacity of the telescope. The reflec- 
tive pfopeftie^ of ellipsoidal and hyperboloidal surfaces are also important. In 
both cases, light striking the surface in a direction toward or away from one focus 
is reflected in a direction either away from or toward the other focus. These prop- 
erties are illustrated in Fig. 
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PROBLEM 21, 



Solutions 



The technique of ray tracing is used in the design of optical instruments. One class of 
such instruments owes its focusing properties to the law of reflection. This law 
states that the angle between the incident ray and the reflecting surface must eo" 1 
the angle between the emergent ray and the reflecting surface. In ray tracing, 
equations are written for the lines containing incident and emergent rays of el " 
tromagnetic radiation. Use this technique to prove the reflective property of th 
parabola, given that the slope of the tangent at the point (xp, yg) on the parabo^c. 

— 4px is 2p /yQ, (The slope of the tangent at any point on a conic section gra^^h 
wiii be derived in Chapter 10.) 

We must show that an incident ray parallel to the axis Is reflected through the focus. 
The geometry of the reflection is shown in Fig. 9.12; since the lines TI and PF 
are parallel, and since the angle of incidence equals the angle of reflection, we 
have that triangle FTP is isosceles and so ^ ^ 20, The equation of the line con- 
taining the incident ray isy - yu. The equation of the line containing the reflected 
ray is (y — yo) ^ tan ^(x - Xq) - tan 26 (x - jco). 




Fig. 9A2 



Since 



tan 26 - 



2 tang ^ 2p 

~7^and tan 0 ^ — , 

1 = tan^e yo 



yo^ ^ yi^ yu-4p'- 



But 



yS ^ 4pXo 



so 



tan 26 - 



4pxo - 4p- 
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Substituting the slope into the equation of the line. 



Itx ^ p,y — 0, so this line passes through the focus. 

Special instruments have been designed to study the eleetromagnetic radiation of 
stars and other astronomical sources in wavelengths outside the visible region. If 
X-fays are to be reflected, the incoming rays must form a very small angle (gfazing 
angle) with respect to the reflecting sufface; otherwise the X-rays are simply 
absorbed. However, with grazing angle incidence, incoming rays that are not par- 
allel to the axis are not focused at all (making it impossible to form an image of a 
source that is not a point) unless an even number of reflections is used. The X-ray 
telescope on the High Energy Astronomy Observatory (HEAO) satellite was 
therefore designed to use two reflections from conic section surfaces, Fig. 9,13 
shows some of the possibilities that were considered. Notice that in each, the 
focus of the paraboloid coincides with one focus of the other conic. 



Foflo cling surfa^ES 



paraboloid = 




hyperboIoid= 



paroboloid 




rGfleeiing surfaces 



paraboloid 




Fig. 9.13 
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PROBLEM 22. In designing an X-ray telescope that uses hypefboloidal and pafaboloidal reflectinf 
surfaces, the most effective placement of the x- and ^^-axes is such that the x-axis 
coincides with the axis of the paraboloid and the^-axis passes through the inter- 
section of the two surfaces. This is illustrated in Fig. 9.14 in cross section. (Note 
that figures 9.14 thru 9.17 are distorted: grazing angles are much smaller than shown 
in these diagraipi.| 
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paraxiol 
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a* If p is the distanca between the vertex and focus of the parabola in Fig, 9,14; c, 
the distance between the center of the hyperbola and each of its foci; the 
distance between the center of the hyperbola and each of its vertices; and k, the 
distance between the center of the hyperbola and the origin, find the equations 
of the two conic sections in this coordinate system. 

Solutions For the parabola, the vertex is at (—(p + c + A), 0) and the focus-vertex distance is 
p, so the equation is 

= 4p (x -Hp + c + ^), 

For the hyperbola, the center is at {—k, 0); the role of the param eter a in the 
standard equation is taken by h and that of the parameter b by Vc^ — , so the 
equation is 

b. Fig, 9,15 shows the ray paths, which form angles a and ^ with respect to the 
^-axis, and the tangents to the parabola and hyperbola, which form angles $ and ^ 
respectively, with respect to the of-axis. Experience in this field has shown that 
when successive reflections take place, surface reflection efficiency is maximum 
when an incoming ray parallel to the axis strikes each reflecting surface at about 
the same angle. Show that this condition, together with the fact that the angle of 
incidence equals the angle of reflection, means that /3 ^ 20, ^ ^ 3d, a ^ 46. 
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Solutioni We number the angles 1, 2, 3, 4, 5 in Fig. 9,15, as shown in Fig, 9,16. 




Fig, 9.16 



We see that 

Zl = ^2 and Z3 — ^5 of incidence ^ Z of reflection); 

Z3 = Z4 (vertical angles); 

Zl = 9 (incoming ray is parallel to axis); 

Z2 = Z3 (for maximum reflection efficiency). 

So 

Zl ^ Z2 ^ Z3 ^ Z4 ^ Z5 ^ e. 

Now, since an exterior angle of a triangle is equal to the sum of the nonadjacent 
interior angles, 

P -20; 

$ ^ p -¥ Q ^ 3$; 
a ^ ^ -¥ e ^ 40. 

c. In designing the surface of the X-ray telescope, the designer must be able to 
express the parameters and h of part (a) in terms of two initial design 

parameters F and >o of the instrument, where F is the distance along the x -axis 
between the origin and the focus of the hyperbola, and yo is the distance along 
the ^ -axis between the origin and the point of intersection to the parabola and 
hyperbola. (These are shown in Fig. 9.15.) 

Recall that under grazing-angle inddence, 0 is a very small angle (this is defi- 
nitely not shown in the figure; the angle occurs where these lines finally intersect). 
It is also true in this situation that p is small compared to F. Show that the 
parameters p, c, and h can be given, at least approximately, in terms of F and yn 
by the following 
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Solution! 



Let (^,, y^) be the point on the hyperbola where the second reflection takes place, and 
iet (0, y.) be the yMntercept of the line containing the ray after the second reflec-' 
tion. (See Figure 9,17,) Then tan a ^yz/Fmd tan ^ ^ yi/(c + k) ^ yj(2c + F), 
since F^k -c. Because 0 is such a small angle, a and p are also small. This 
suggests the following approximations: 



yj^yo^yu tan « ^ a; tan ^ * /3, 



(0. yj) 




Fig. 9.17 



(Fig, 9.17 does not show this because the angles are not small enough. The reader is 
encouraged to imagine how the figure would change if angles 6, and a shrink.) 



Since a ^ 2p, we get 



and if ^2 ^yi, then 



So 



yi ^ 2yi 
F 2c + F ' 



2F ^ 2c + F, or F - 2c, 



1 - % 

c i -F, and F + c ^ |F 



The parabola has equation ^ 4p (x + p + c ^- ^), Since y ^ y^ when ^ ^ 0, and 
since c -¥ k ^ 2F, 



Kp %' X, 



yo- - 4p(p +2F)^ 4p2 + SFp. 
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If p is small compared to F\ the term in may be neglected giving p ^ y,r/(8f ): 
h can be expressed exactly in terms of p. c\ and k by observing that (0, y„) is on 
both the parabola and the hypefbola. This means thot 



4p (p + c + k ) 



and 



/r 



c 




Then 




so 



or 



h-{c" + + 4p (p + c + A')) + c'A-^ ^ 0. 



Using the quadratic formula, 

, _ + fe- + 4p (p + c + ± V(c^+ + 4p (p + c + ^))^ - 4cH- 



The positive sign before the radical is discarded, since it produces a physically 
unraalistic value of h larger than c, so 

2/r - c= + + 4p (p + c + it) = \/(P^TT^TlJ'^~TTTT)P"^lc^ , 

This analysis produces an initial set of parameters p, k, and /j. A ray-tracing 
computer program, based on the principles discussed in Problem 21, is then used 
to check the actual focusing capabilities of a hypothetical instrument with these 
specifications. The use of the computer with such a program makes it possible to 
refine the values of the parameters for best focus under desired conditions. Actual 
preparation of the reflecting surface is also controlled by computer-driven 
machinery once the optimal values of the parameters are established. 
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Although calculus is used extensively in space science and technology, we shall 
consider in this chapter just a few problenis, most of which extend or 
aniplify Ideas discussed in previous chapters. Calculus is also used in the 
Appendix, 

PROBLEM 1* Until recently it was accepted that there were three possible states in which matter 
could exist: solid, liquid, and gas. Under conditions that normally prevail on 
Earth, these are the only states in which matter is found. However, it is now known 
that if the temperature is very high or the density is very low, a fourth state of 
matter can exist; it is called plasma. A plasma consists of electrons and positively 
charged ions rather than neutral atoms, and so it has both electric and magnetic 
fields. (An ion is an atom that has lost one or more of its electrons,) On Earth, 
plasmas exist, at least temporarily, in lightning, electrical sparks, fluorescent 
lamps, and in the ionosphere. 



In addition to the electromagnetic radiation we sense as heat and light, it is now 
known that the Sun emits partiela radiation having a wide range of energies. The 
particles (or plasma) appear to come from specific regions on the Sun, some as 
highly energetic particles which move radially outward into interplanetary space. 
Some of these highly energetic particles that reach Earth's ionosphere produce 
auroral displays (the northern lights) and affect shortwave radio transmission by 
modifying the ionospheric structure. 

A lower energy component of the particles is emitted from the Sun on a continuous 
basis, and these lower energy particles also move away from the Sun in a straight 
line (radially). The study of this interplanetary plasma, which has been called the 
solar wind, is of great concern to astronomers and other scientists for several 
reasons. One is that the Sun is the only star we are close to, and the emission of 
plasma means that it is very gradually losing matter, an important factor in stellar 
evolution. Another is that the plasma state of matter is difficult to study on Earth 
because it is hard to reproduce in the laboratory the conditions of high tempera- 
ture and low density that exist naturally in the solar atmosphere and in interplane- 
tary space, 

A number of space probes and satellites have been used to investigate the prop- 
erties of the interplanetary plasma. The Interplanetary Monitoring Platform (IMP) 
series of probes from 1963 to the present, the Orbiting Ceophysical Observatory 
(OGO) series from 1964 to 1974, the International Sun-Earth Explorer (ISEE) 
satellites from 1977 to the present, and the Mariner^ Pioneer^ and Voyager deep» 
space probes have all carried experiments resulting in a series of measurements of 
flow direction, density, velocity, and electric and magnetic fields of the 
solar wind, 



It has been postulated, on theoretical grounds, that the magnetic field lines of the 
solar wind coincide with the locus of particles emitted from the Sun, and the 
experimental findings to date seem to support this hypothesis. 

a. Determine the shape of this locus, given that the solar atmosphere from which 
emission takes place rotates at a constant angular velocity and that particles 
move outward with constant velocity in the radial direction. Assume the direction 
of rotation is clockwise* 
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SDlution: 



An intuitive solution using **time lapse'' polar graphing is displayed in Fig. 10.1, and 
this shows that the locus is an Archimedean spiral. This can be verified analyii- 
cally using calculus. We are seeking r - f(0) such that dr /dr - (particles emitted 
with constant radial velocity), and dB/dt ^ C (emissive origin is rotating with 
constant angular velocity). From the chain rule. 



dr 
dr 



dr 66 
de dt 



and substituting from above, we have 



de 



dr 

de 



V 

c 



a constant we may call k. 

Integrating and choosing the coordinate system so that f(0) ^ 0, we have r ^ kS, 
which is the equation of an Archimedean spiral. 




Fig. 10.1 
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b. It has bean observed that the equatorial region of the solar atmosphere rotates 
at a rate of C " 2,94 x 10^^ radians per second with respect to the distant stars. 
This is known as the sidereal rotation rate and is equivalent to a sidereal rotation 
period of 



2lT 



2S4 X 10 



ns ^ 2.14 X 10' s ^ 24.7 days. 



Spacecraft measurements of the solar wind velocity show time variations, with 
velocity peaks at approximately 25-day intervals as welL For example. Mariner 2 
measured velocities varying from 400 km/s to about 750 km/s at some peaks. 
Determine k and plot the graphs of the Archimedean spirals for velocities of 
400 km/sand 750 km/s. 



Solution! 



For V 400 km/s, 

K^V^ 400 
C 2.94 x 10^^ 



^ 1.4 X lO^km/rad. 



For V ^ 750 km/s, 



k ^ 



750 



2.94 X 10^^ 



^ 2.6 X lO^km/rad. 



The graphs of r - kB are shown in Fig. 10,2, (Note that it is the practice to use 
km/rad as th& unit for k, but this is equivalent to km, since the radian is 
dimensionless,) 




I Sir 



Scale for r: 2 units = 10** km 




(a): V^400 
km/sec 



(b): V-7S0 
km/sec 



Fig. 10.2 



In Chapters 4 and 7, we considered some of the corrections needed to produce 
undistorted pictures of spacecraft observations. Here is another such correction. 
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PROBLEM 2p Most satellite photogfaphy makes use of scanning techniques. This is illustrated in 
Fig. 10.3(a), where the scanning is done in the difaction orthogonal to the flight 
path. In performing the scan, a system of mirrors and lenses rotates around an axis 
parallel to the flight path. Although the scanning system rotates at a constant 
rate, we can see from Fig. 10, 3(b) that the rate at which the scanning beam moves 
along the ground depends on the angle it makes with the vertical, 




(a) 

Fig. 10.3 



If wa imagine that a square on the ground has the pattern shown in Fig, 10.4(a), 
the result of this variable Earth-scan rate will be the distorted pattern shown in 
Fig. 10.4(b). In order to produce an undistorted picture, the actual recording of 
the images must be done at the Earth-scan rate rather than the rotation rate. 
This panoramic distortion correction requires the ability to express the scan rate 
along the ground, dx fdt in terms of the satellite height, h, the angle and the 
rotation rate, d0/dr, of the scanning system, where fl, A, and x are as defined in 
Fig. 10.3(b). Find such an expression. 
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Solution: (This is a straightforward related rate problem,) From Fig. 10 J(fa), we see that 
tan ^ ^ = , so that x = h tan a Differentiating with respect to time, 

f = /isec^^f , 
dt dt 

PROBLEM 3, In Problem 9a of Chapter 4, we estimated the surface area of an antenna '*dish'' by 
tfeating it as though it were a flat circle. In actual practice, such an antenna is a 
paraboloidal "cap" whose depth is from 10 percent to 20 percent of its radius. Let 
us see how good this estimate was. 



a. Find an expression for the surface area of a paraboloidal cap that is bounded 
by a circle of radius r and has depth a. 

Solution: We can consider the paraboloidal cap to be generated by revolving the illustrated 
portion of the parabolic curve (Fig, 10.5) around the ^-axis. We must first deter- 
mine the function y = f(x) for this curve. Then the surface area will be given 
by the integral 



5 - 2w 



2wJ^ f(^)V^ 



1 + [f(.t)]-dv 



Since this curve is a parabola with axis horizontal, op ening t o the left, and 
with vertex at (fl,0), its equation has the form >^ ^ bVa - x, wher e b mu st be 
determined so that (0,r) satisfies the equation. This means b\ ^a = 0. or 
b = r l\fa. The function we need, then, is t(x) = (r /V5)Va = x = rVl ^ x /a. 
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Fig, 10 J 

Then 



fix) ~ 



(4 



2 ^ 1 



S « 2'Trf% s/l - - Jl 
Jo V a V 



..'4 



r s i - i VI + , . ' 
^0 V a V 4a~ = 4ax 



^ lit \ Vr^ + 4a^ — 4ax dx 
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This integral may be evaluated by the substitution « - + 4a^ = 4ax to produce 



OS' 

b. Recall that the paraboloidal cap had a radius of 10 meters. If its depth was 1 m, 
find its surface area (assuming exact numbers) and then find the relative error of 
the estimate made in Chapter 4. 

Solution: For r ^ 10 m and s ^ 1 m, 

5 - ^[(104)-^^-^ = 10^'] * 317 m^ 

The estimate* approximating the paraboloidal cap as a circle, was 
S ^ irn ^ 100^ = 314 m% 

f,e, ^ ^^^g^^^^^ ^ 0,0098, or about 1 percent, 

c. Find the relative error for a depth of 2 m. 

Solutiom S ^ ^[(116)^^- 10^] ^ 326 m^ 

error = |326 — 314| m" = 12 m^ 
r.e. = ^ 3-8 percent 

PROBLEM 4. In Problem 1 1 of Chapter 7, we observed that a spacecraft at a distance h from Earth 
can observe only a portion of Earth*s surface. 

a* Derive a formula for finding the fraction of the observable area as a function 
of height above Earth*s surface. 

Solution: The portion of Earth's surface visible from the spacecraft is shown shaded in 
Fig, 10,6. Let ^4^ be the area of the zone with altitude BE. If we set up a 
rectangular coordinate system with origin at A, and if the coordinates of a point on 
the arc EC are (g(y), y)* then this surface area is found by evaluating the integral 



A, = liT f ' g(y)Vl + [g'iy)f dy 



where g(>^) = x = VA^ - >»' and andyE are the y-coordinates, respectively, 
of the points B and E, 
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Fig. 10.6 



To determine tha>'-coordinate of B,}'B,we observe thut triangles ABC and ACD 
are similar, so that 

AD, AC 
AC AD 

or 



giving 



We have 



R ~ R + h 



R + h 



and 



so the integral is 



R ^ h 

If we let represent the area of Earth'ssurface, then - 4fTR ^ so that 

_ _ 2^R^l\ h 
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b. In ^pfil 1983, two members of the Spaca Shuttle Challenger crew, Story Hub- 
grave and Donald Peterson, performed an extravehicular activity (an activity 
outsid: ejhe spacecraft) while Challenger was at an altitude of 280 km. What frac- 
tionojf Earth did they see? (Use 6380 km for the Earth's radiusO 

Solution; For h = 280 and R = 6380, 

% " 2(638?+ 280) " or 2.1 percent, 

c, Discuss the manner in which the fraction AjA^ varies with the altitude h. 

Solutions Intuition suggests that as h increases, the value of ^^M^ should vary from mto to 
1/2, O XI the surface of Earth, the fraction is zero. As h increases, so does the 
fractio ri, and yet it must always be less than 1/2 ; that is, one cannot hope to view 
more tlian a hemisphere at any one time, A little algebra bears this out 



A^ ^ h 

2(r+/i) 

is certainly zero when h = 0, Observe that 




k^h inereases, the denominator of the right-hand side decreases, which forces 
the entire fraction AJA^ to increase. Furthermore, as h rlh — * 0, and con- 
sequeritly ^^/y4g approaches 1/[2(1 +0)]= 1/2. 

L At what altitude will an astronaut see one-fourth of Earth's surface? 
Solution: We must find h such that 

1 ^ h 

4 2(6380 + /i) 

4^/1^ 2 (6380) + 2 v/i 

2- /I ^ 2(6380) 



/2 ^ 6380 km. 



e, The first astronauts to travel that far from Earth were Apollo 8 crew 
(Anders ^ Borman, and Lovell), v^ho orbited the Moon on Christmas Day, 1968. 
What percent of Earth's surface could these astronauts see u they passed the 
Moon, a distance of 3.76 x 10^ km from Earth? 

Solution. ^ ^ 3-76 x 10^ ^ 3.76x10^ 

2(3.76 X 10^ + 6.38 x 10^ 2 (3.82 x 10^) 

« 0,492, or 49.2 percent 
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PROBLEM 5* The reflective properties of the conic sections were discussed in the final two prob- 
lems of Chaper 9, where the formula for the slope of the tangent to a parabola 
was used. Use differentiation to find such a formula for each of the following 
conies at a point (ao, jn), 

a. The parabola = 4px 



Solutions Differentiation produces 



dx y 



So the slope of the tangent at (^o, jo) is ~ , or, if stated in terms of Xn, 

slope 



b. ihe ellipse — + - 1 



Solutions Differentiating this equation. 



a' b'dx - 



The slope of the tangent at (^n, >^n) is then 



a-y^ 



c* The hyperbola ^ ^ = 1 

Solution: This is exactly the same as (b) with one sign change , so that the slope of the tangent 
at (xo,>^o) is^ — , 



PROBLEM 6. Our final example in this chapter is a result of the spectacular discovery by the 
Voyager space probes that one of Jupiter's moons, lo, is the site of active volca- 
noes, the first known instance of volcanoes other than those here on Earth, Images 
returned by the spacecraft have provided measurements which scientists are 
using to develop and evaluate models by which both the behavior and possible 
causes of this volcanic activity may be understood. 
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The starting point is the familiar projectile problem with the valid condition (in 
contrast to such pfoblams applied to Earth) that air resistance is neglected, since 
So has no atmosphere. It is also convenient to begin the modeling process by 
assuming that the gas and solid particles ejected from the volcano's opening do not 
affect each other's motion upon ejection, that the opening is circular (roughly), 
and that all par tides start from the same point below the surface with the same 
initial velocity at all possible escape angles in any direction. These assumptions 
together have been called the balltsiic model. 

a. Photographs of the plumes of some of lo's active volcanoes show that any verti- 
cal cross section through the volcano vent's canter has the shape diagrammed in 
Fig. 10.7, where we have placed the r-axis along lo's surface and the z-axis per- 
pendicular to the surface. If (r, z) is a repfesentative point in such a cross sec- 
tion, find an expression for the escape angle / (or a trigonometfic function of i) 
which will cause an ejected particle to pass through this point (i is measured from 
the 2 -axis), Let fp be the radius of the circular opening, Vo the initial velocity of the 
particle, and d the distance below the surface of the point where the 
particle originates. 



Solb^'on* Iff is the value of lo's gravity and /o the time at which the particle is ejected, then our 
model is given mathematically by the following conditions: 




Fig. IQJ 



dfz 



dz 
dr 



Vocos /■ z (/fl) ^ -d 



d/= 



dr 
dr 



Vosin /; r (to) - 0 



Routine integration and application of the initial conditions result in 



(A) 



^gt^ + gtat + vof COS i - d - /J - i'n/o cos / 



(B) 



r(0 



(voSinj)(f ~ to) 
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Chapter Ten 



To get an expression for i\ we eliminate t by solving equation (B) for f and substi- 
tuting this expression into equation (A), which gives 



(C) 



2 = —^coseri + rcot I — d. 
2v& 



Since equation (C) contains cosec^i and cot i, we may use the identity 
cosec'/ ^ 1 + cot^i to obtain an equation in cot i, namely, 

— jCOt^i ^ r cot I + 1^ + d H- z ^ 0. 



Applying the quadratic formula* 



cot i ^ r ^ 



and simplifying, 
(D) 



cot I 



We see that there are two possible ejection angles that will bring a particle through 
a particular point, one on the way up and the other on the way down. This is 
illustrated in Fig. 10=7, at the point marked (r, z). 

If Iq is the largest possible escape angle for the vent, express the height of the 
lower boundary of the portion of the plume that contains both upward- and 
downward^moving particles (identified in Fig, 10.7) as a function of r, and then 
find and of Fig. 10,7. 

Solutions From equation (C), with i ^ 



- =^QQ%ec% + rcot lo - d; 

and rp are the values of r in this expression for which - 0. Setting ^ 0 and 
multiplying by 



2vi sin-iQ 
8 



produces the equation 



^2 _ 2y|sin/oCGsiQ ^ ^ 2v| d sin% _ ^ 
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Using the positive sign in the quadratic formula to get the larger r. 



. ''m ^ -r- sin lo 

' * a 



cos 



. ^ I , . 2gd 
h + cor to - ^ 
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and using the negative sign to get the smaller value. 



tp - ~ sin !q 
8 



cos 



I 2gd 



c* Express the height Zm of the upper boundary of the plume of Fig. 10.7 as a 
function of r and find an expression for the maximum height of the plume. 



Solutions The upper boundary is the set of points for which the two solutions given by 
equation (D) of part (a) coalesce; in other words* where the radical vanishes, or 



\VqV Vq 

SO 

The maximum height hp occurs when r — 0: 
(E) ''^ " i - 

d* Express the coordinates (n , zi) of the point to the right of the z-axis, at which 
the upper and lower boundaries meet, in terms of Iq, Vg, and d. 

Solutions Since this point is on the upper boundary, we have the radical of equation (D) equal 
to 0, and since it is on the lower boundary, we have / = iq. From equation (D), 

cot la ^ — , so n ^ = tan iq. 



From our solution to part (c). 



2g 2vi 



or 



^1 - 



= ^ (1 " tan%) ^ d. 
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e, Mea surements of images obtained from the Voyager J imaging system have pro- 
vided v-^Iues for r^, r^, and hp. Use the foregoing results to obtain expressions so 
that In, ^'ilg, and d may be calculated from these measurdments. 



Solution; From equation (E), d ^ 
known ^ 

From part (b). 



(vi}/2g) =- hp. This will give d if we have v,T/g, since /ip is 



~^ rp ^ sm 1(1 cos li) 



or 



i^m + ^p) cosec(2/n). 



This wii3 give ^f/g if we have /,h since and are known. 
Also fr^ in part (b), 



Sin ln\/C0S'in = ^ 



Since ^ = (v^/2g) - hp. 



- — sm in\/C0S-/n ^ 1 + 



2vS 



Sin^in 



cos h 



Squaring to ehlminate the radical and substituting * 2 sm /„ cosji 



rm + rp 



2ghp 



, 2 ^/ 2sininCOs4 \ , . 



COS^In 



4/ip , . 
~— - tan lu ~ tan-Ji), 



COS"/(i 



which we can solve for tan in using the quadratic formula: 
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Since /n appears on the images to be smaller than 45-, the minus sign is 
chosen, giving 



tan In 



f. Voyager 1 delected eight volcanic plumes, of which plume 1 and plume 3 were 
closest in shape to our diagram in Fig. 10.7. Table 10.1 gives the observed mea- 
surements of Tm, ^p, and /ip for these plumes. For each, calculate /n, vUg^ and d 
from these data, then use the results of part d. to predict n and Zx for these plumes. 



Table 10.1 

MeesiirEnients of lo 'i Plumes 



Plume 








No, 






hp (measured in km) 


1 


500 


17,5 


280 


3 


125 


7,5 


70 



Solution: For plume 1 , 

jan /n 
.\/n ^ 28,05 



f 



(517.5) (cosec (56, P)) ^ 623 km 

623 



d-^^-h,^ 2 



280 ^ 32 km 



tan /o ^ (623) (0,533) ^ 332 km 



- ^ (1 - tan%) ^ fip[l = (0,533)= 



8 

For plume 3, the results are as follows: 

tan lo ^ 0,4821, so io - 25.74* 



- 32 ^ 191 km. 



— ^ 169 km; d ^ 15 km; n ^ 82 km; ^ 50 km 
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g, Fofplume 3, it was also possible taio measure This measurement was 125 km. 
Compute the relative error of our ca^lculated value for n as a percent. 



This error is large enough to deniaji^ refinement of the modeL The assumption 
leaitliWy to hold is that ejected ga^ and particles do not affect each other's 
motiooiit is more probable that the m combination of particle sizes and rate of gas 
flow i3iiich that the particles are cprs-ried by the gas into the central portion of the 
top cfthe plume and released Into ballistic trajectories only on descent 

Althoyih we shaU not consider thes^ modifications here, the reader may be 
mteieitedin knowing that the result^ of further refinements of the model are 
consiM with the theoretical pfopo osal that lo's volcanoes are due to tidal 
effeemnits surface generated pfima^rily by another Jovian moon, Europa. 



Solutions 




Relative errOr 



43 
125 



0.34 ^ 34 percent. 
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GRAVITATIONAL FORCES 
AND THE CONIC 
SECTION TRAJECTORIES 
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ha elliptical shape of planetary orbits was first assefted by Johannas Kepler 
on the basis of painbtaking observations made by him and by his prede- 
cessor, Tycho Brahe. It was Isaac Newton's great achievement to establish 
mathematically that the inverse square law force of gravitation must produce a 
trajectory that is one of the conic sections. We present this analysis as an appendix 
rather than in the calculus chapter (Chapter 10), since the manipulations needed 
in the development include some complexities that may be unfamiliar to the 
intended audience. 

Before beginning the main problem, we need to establish some properties of fifst 
and second derivatives of the vectors, The unit vectors 7, Ur, Ub are shown in 
Fig. A.l, We have 

Mf — / cos 6 + j sin $ 

and 

- T (=sin 6) + y cos 0. 
Differentiating with respect to time, 





In polar coordinates, the position vector is (r m j, and so 
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acce 



= ^ d?T dj d7"« ^ ''iH "» ^ 'IFTT 

"id,' 'id?ii«'+l2dfl7 + ''df5i"»* 

Wa now consider the statamant of the law of gravitation: The force on the orbiting 
body (mass m) is pfoportional to the product of the maiias of the two bodias 
involvad and inversely proportional to the square of the distance between thein. 
In symbols, 

d^(rSj _ CMm - 



6!^ 



whara the minus sign expresses the fact that the force of gravitatipn acts in the 
direction toward the mass M. Using the result above for acceleration, and cancel- 
ing the m, we get 



"dr 



"'IdFiJ^^^^L^drdF^^^J 



Since the ug term is missing on the right, we have 

^drde^ d20 ^ 

Kepler's first law, that the radius vector sweeps out equal areas in equal times, 

1 T 

can be stated mathematically ^ = A, a constant. If we differentiate this 
expression with respect to time, we get 



^ drde ^ ,dV 



^ 0 



or 



1 

2^ 



» dr d| dfl 
dtdt ■ df i 

idB 



^ 0 



So we may let the constant A ^ , and this is equivalent to the fact that the i4b 
coefficient vanishes, or Kepler's first law. 

Equating the u, terms, g _ , = = ^ ^ 
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Using the substitution B - CM, as well us .4 ^ r^— , proclucuH the differential 
equu 

(1) ^^dl^M^o. 

dr r 

We are seeking r as a function of d\ it turns out lo be easier to find 1 1 r as a function 
of 0. This can be done by letting w — l/r. or equivalently, r — 1/ h'. Then 

dr ^ _^dw ^ 1 dw dO 
di di do dt 

. dw 1 dd . do . 

do di dt 



Differentiating again. 



dfr ^ _ ^ ^ fdw \ ^ dfu; dO 
dr dt \dB ) ~ dO^ dt 



de- 



We can now rewrite (1) as 



d^' 



or 



d^iv , B 

- — ^ + IV ^ — ; , 

de- A- 

which has a solution, w — C cos (0 ^ ^i?) + . 

A- 

But this means that 

, ^1 1 A^B ^ 

' IV Ccos(0 - 6,0 + (BM^) ' 1 + {CAyB)cQsi6 - B^' 

¥ov Arl B - ep,C ^ = 1 /P. $i ^ 0. we get r ^ epi {\ — e cos 0, the equation of 
the conic section in Chapter 9, Problern 4. 

We next derive the ''vis-viva," or energy integral, and show that the eccentricity of 
a conic section trajectory is physically determined by the total energy of the 
gravitational system. Again, it helps to first establish some properties of the vec- 
tors involved » In this context, we will need the square of the velocity vector, v\ 
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so that 



= dv Id,.. 

The gravitational equation may be written in the form 

dv GMm 

Premultiplying by y, using the dot product, produces 

- dv GMm - - 



Substituting 



on the right. 



- dr^ ^ d0^ 



on the left, and recalling that u, ^ u, - I while Ur ^ Uy ^ 0, we get 

fl d . ,J GMm dr 



di 



Integrating, we get the energy equation 

/n\ 1 ^ GMm = 

U) jm y- - ^ E, a constant. 

In physics, the first term on the left, (1/2) m V", is the kinetic energy of the system; 

the second term, - GMm ^ gravitational potential energy; E, the constant 

of integration, is the total energy. We may evaluate E by considering a particular 
point in the orbit. Since r - epi (1 - e cos B), r attains its minimum value for 

e ^ ^: 



^ . = 

* mm 



1 + e 
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1 



so 



ep 

Now from the vector expression for velocity, 

^■=(fr-(4f)" 

But when 6 = ir, the velocity is entirely in the 5^ direction, so for w — Wm^x^ 
= ('^^f ) equation (2) becomes 

^mAh^^ — B m w — E — 0 , 

^ do 

recalling that A ^ ^" ^ ^ ^ ^^^^ Now, using the quadratic formula with 
the positive sign, 

Bm + V'iBmf ^ImA^E B ( K 2A^\ 



mA^ A^\" ' B-m 

Equating the two expressions for w^^^, and recalling that ep = A^/B'. 



and tharefore 
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Since m > 0, v c see that the nature of the trajectory depends on the total energy 
E: 

If E - 0, then e ^ 1 and the trajectory is a parabola; 
if £ < 0, then e < 1 and the trajectory is an ellipse; 
if E > 0, then e > I and the trajectory is a hyperbola. 
It is sometirnes useful to solve (3) for E, giving 

- B-m (e' - 1) 

^ GMm (e^ - 1) 
2 ep 
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Apogee 

Attitude of a spacecraft 

Celestial equator 
Celestial sphere 
Center of mass 

Conjunction of planets 

Cosmology 
Declination 
Direction cosines 

Ecliptic 

Electromagnetic spectrum 

Ephemeris 

Geocentric 
Iteration 

Jovian 
Julian day 

Orbital period 
Perigee 

Photon 



The most distant point from Earth reached by a body in an elliptical orbit with Earth 
at the primary focus. 

The orientation of the spacecraft in space, with respect to some chosen coordinate 
system. 

The projection of the equatorial plane of Earth on the celestial sphere. 

An imaginary sphere of infinite radius on which celestial objects appear projected. 

The point within a body at which all the mass could be located without changing its 
dynamical behavior. 

The position of the planets when they are on the same right ascension circle on the 
celestial sphere (in other words, when they appear closest together in the sky). 

The study of the evolution of the cosmos or universe. 

The analog, on the celestial sphere, of latitude circles on Earth. 

The cosines of the angles made by a vector in space with each of the three positive 
coordinate axes. 

The path described by the center of the Sun on the geocentric celestial sphere during 
the course of a year. 

Radiation of various wavelengths emitted in the form of washes carrying rapidly vary- 
ing electric and magnetic fields (light is an example of a portion of the spectrum). 

A list of the successive positions of a celestial object on the geocentric celestial 
sphere for a series of equally spaced times. 

Concentric with Earth, 

A repetitive mathematical procedure, on an initially chosen trial value, which can 
produce improved values of a desired quantity. 

Relating to the planet Jupiter. 

The number of days, and fraction of a day, measured from noon on 1 January of the 
year 4713 B.C, 

The time it takes for an object to complete one orbit. 

The closest point to Earth reached by a body in an elliptical orbit with Earth at the 
primary focus. 

The smallest unit (or ^'particle'') of electromagnetic radiation, carrying one quantum 
of energy. 
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Pitch 



An angulaf rotation of an aircraft or spacecraft around an axis through the wings 
(which has the effect of moving the nose up or down). 



Right ascension 

Right-handed three- 
dimensional coordinate 
system 



Roll 

Watt 

Yaw 



The analog, on the celestial sphere, of the longitude circles on Earth. 

A convention of three basis vectors most simply represented by the thumb, index, 
and middle fingers of the right hand when extended at right angles to each other; the 
X-axis is along the thumb, the y-axis along the index finger, and the z-axis along the 
middle finger. 

An angular rotation of an aircraft or spacecraft around an axis along its length (which 
has the effect of tipping its wings). 

A unit of power, (In the mks system of units, 1 watt — 1 joule per second,) 

An angular rotation of an aircraft or spacecraft around an axis perpendicular to its 
body (which has the effect of moving the nose from side to side). 
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